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Voorwoord 


Door het toenemende gebruik van microcomputers maken de grafische toepas- 
singen een sterke ontwikkeling door. Het beeldscherm is tegenwoordig, veel 
meer dan de printer, het afdrukapparaat bij uitstek. Grafische toepassingen 
(graphies) vinden we bij computerspelletjes, bij computerkunst en in toene- 
mende mate op de werkplek van ingenieurs, ontwerpers, architecten en 
anderen die zich bezighouden met Computer Ondersteund Ontwerpen (Com- 
puter Aided Design). 


De geïnteresseerde leek is vaak verrukt van de mooie 'plaatjes', die bij 
demonstraties getoond worden, en denkt dat hiervoor hele ingewikkelde pro- 
gramma's nodig zijn. lemand, die iets van wiskunde afweet, weet hoe een 
functie eruitziet en kan doorgaans zonder veel moeite grafische programma's 
maken. 


Ik heb voor het Zwitserse Computertijdschrift 'MIKRO+KLEINCOMPUTER' een 
cursus ‘programmeren met hoge resolutie! geschreven, die in een aantal afle- 
veringen van het blad is verschenen. Al direct na het verschijnen van het 
eerste artikel werd ik werkelijk overstelpt met enthousiaste reacties van de 
lezers. Dit heeft ons aangezet tot het maken van dit boek. 


Veel van de grafische toepassingsprogramma's worden geschreven voor een 
bepaald merk en type mierocomputer of voor een bepaalde plotter. Wie een 
dergelijk programma wil herschrijven voor een ander merk microcomputer zal 
veel moeilijkheden ondervinden en vaak zelfs zijn pogingen staken. 


De veertig programma’s in dit boek zijn geschreven in standaard BASIC en 
gebruiken geen POKE- en PEEK-opdrachten. Bovendien worden slechts twee 
grafische opdrachten gebruikt, te weten het inschakelen van de hoge-resolu- 
tiestand en het verbinden van twee punten door een rechte lijn. De illustra- 
ties die als voorbeeld dienen van hetgeen de programma's tekenen zijn 
gemaakt op de miniplotter van de Sharp PC-1500. 





Wij hopen dat de lezers veel plezier aan deze programma's zullen beleven en 
dat zij voor velen een aanmoediging zullen zijn om het terrein van de compu- 
tergraphies verder te verkennen. Mijn dank gaat uit naar de heer H.J. 
Ottenbacher van Micro*Kleincomputer; zonder hem zou dit boek nooit het 
daglicht hebben gezien. 


Marcel Sutter 
voorjaar 1984 


Voorwoord bij de Nederlandse 
uitgave 


De auteur, Marcel Sutter, heeft dit boek oorspronkelijk geschreven 
in standaard Microsoft BASIC. In zijn boek is een aantal appendices 
opgenomen met programmalistings voor een aantal microcomputers, 
zoals de Commodore 64, Apple II. Wij hebben besloten dit boek voor 
een zestal microcomputers afzonderlijk uit te geven. Het boek dat u 
nu leest is de ZX Spectrum versie. Alle programma's uit dit boek 
zijn getest en hebben gedraaid op een ZX Spectrum. De programma's 
zijn vanuit het Spectrum-geheugen op een printer afgedrukt en zo 
in het boek opgenomen, Tik ze dan ook precies zo in als ze in het 
boek zijn afgedrukt. 


Wij raden u aan al bij de eerste programma's te experimenteren met 
een aantal grafische mogelijkheden van de ZX Spectrum. De pro- 
gramma's gebruiken allemaal een zwart potlood (INK) op wit papier 
(PAPER), waarmee de zwart-op-wit tekeningen gemaakt worden. 
Probeer gerust andere combinaties door de opdrachten INK en 
PAPER te gebruiken. Met FLASH en BRIGHT kunnen ook leuke 
effecten bereikt worden. In de appendix vindt u een BASIC- 
programma voor het inlezen van een machinetaalroutine, waarmee 
vlakken gekleurd kunnen worden. 


Wilt u op uw printer laten afdrukken wat een programma getekend 
heeft, zet dan vlak voor de STOP-opdracht de opdracht COPY in 
het programma. Hierdoor wordt het schermbeeld op de printer afge- 
drukt. 


Het aantal beeldpunten (pixels) dat de ZX Spectrum gebruikt is 
256 horizontaal (genummerd 0 t/m 255) en 176 (genummerd 0 t/m 
175) vertikaal, Veel illustraties in dit boek zijn gemaakt door de 
auteur op een plotter met 220 beeldpunten horizontaal en vertikaal. 
Zij geven dus wel een goede indruk van wat u op het scherm kunt 
verwachten, maar zijn niet precies gelijk aan hetgeen u op uw ZX 
Spectrum scherm zult zien. 


Alle programma's zijn kort. Verfraai ze, verbeter ze en breid ze uit, 
De mogelijkheden zijn enorm. Omdat de programma's opzettelijk klein 








gehouden zijn, wordt er bijvoorbeeld niet getest of bepaalde waar- 
den, die u moet intoetsen, het programma door een of andere fout 
laten afbreken. Een verbetering zou zijn om deze tests wel op te 
nemen, waardoor een programma alleen met de juiste invoer aan 
het werk gaat. De appendix bevat een aantal suggesties voor het 
verfraaien en uitbreiden van de programma's. 


Het is eigenlijk een boek voor elke ZX Spectrum bezitter. Vindt u 
de meestal wat wiskundige uitleg bij de programma's teveel van het 
goede, tik de programma's dan gewoon ìn en draai ze; u zult ver- 
rukt zijn van het resultaät. Kunt u zich nog iets van de sinus en 
cosinus herinneren of zit u nog op school dan zult u weinig moeite 
hebben met de theorie in dit boek. Voor leraren en leerlingen zijn 
wellicht de programma's voor het tekenen van functies interessant. 
Het boek behandelt ook het principe van driedimensionaal tekenen 
met verborgen lijnen (hidden lines). Erg leuk zijn de vijf program- 
ma's waarmee getekend kan worden zoals dat in de taal LOGO kan. 
Dit geeft zeer fraaie tekeningen, vooral de Turtle-graphics. Ook 
bevat het boek een vijftal educatieve toepassingsprogramma's, waar 
in graphies de hoofdrol spelen, 


Tot slot willen wij Nico Huffels bedanken voor de assistentie bij 
het bewerken van de programma's in dit boek. 


Nok van Veen 
augustus 1985 
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1 ZX Spectrum Graphics 


Het beeldscherm van de ZX Spectrum is als grafisch scherm te 
gebruiken, waarbij horizontaal 256 en vertikaal 176 puntjes (pixels) 
aan of uit gezet kunnen worden. Op dit grafische scherm kunnen 
we zogenaamde hoge resolutie graphics tekenen met behulp van de 
POINT-, PLOT-, DRAW- en CIRCLE-opdracht. 


In alle programma's gaan we uit van de standaardinstelling voor de 
achtergrondkleur (PAPER) en de afdrukkleur (INK), de papier- en 
schrijfkleur. De standaardinstelling is wit voor PAPER en zwart 
voor INK, We tekenen dus met een zwart potlood op wit papier. U 
kunt echter bij het tekenen gebruik maken van acht kleuren (zie 
pagina 22). Experimenteer in uw programma's met PAPER en INK en 
kijk wat het effect ervan is. Ook kunt u met BRIGHT en FLASH 
beelden helderder respectievelijk flikkerend maken. In de appendix 
staat een programma met een PAINT-routine voor het kleuren (vul- 
len) van vlakken. Dit kunt u bij de programma's in dit boek gebrui- 
ken. 


De PLOT-opdracht wordt gebruikt om aan te geven waar het tekenen 
van een lijn moet beginnen. Met de DRAW-opdracht wordt vervol- 
gens aangegeven naar welk punt de lijn getrokken moet worden. De 
oorsprong van het grafische beeldschermcoördinatenstelsel ligt bij 
de ZX Spectrum in de linkerbenedenhoek. De x-as wijst horizontaal 
naar rechts en de y-as wijst vertikaal naar boven. De 256 beeldpun- 
ten op de x-as worden genummerd van 0 tot en met 255 en de 176 
puntjes op de y-as van 0 tot en met 175. Een punt op het beeld- 
scherm (een pixel) kunnen we dan aangeven met een x-coördinaat 
en een y-coördinaat, bijvoorbeeld het punt (128,88), dat in het 
midden van het scherm ligt (zie p.2), 


Om de ZX Spectrum een lijn te laten trekken tussen de punten P, 
en P‚, met respectievelijk coördinaten (x1,y1) en (x2,y2) geven we 
de volgende opdrachten: 


PLOT x1,yl 
DRAW x2-x1,y2-y1 
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N 
*255 


Als (x1,y1) het punt (55,70) is en (x2,y2) het punt (150,135), 
dan worden deze opdrachten: 


PLOT 55,70 
DRAW 95,65 


De lijn wordt dan vanuit P, naar Ps getrokken. Willen we vanuit 
P; naar Pj een lijn trekken. dan nemen we 


PLOT x2,y2 
DRAW x1-x2,y1-y2 


wat voor ons voorbeeld neerkomt op: 


PLOT 150,135 
DRAW -95,-65 


We zouden de lijn ook kunnen tekenen door een heleboel puntjes 
tussen P, en Pz, die precies op de verbindingslijn van P, en P 
liggen, te laten oplichten. Deze methode staat bekend als ‘punt- 
graphics'. Wij gebruiken echter in dit boek de 'lijngraphics'- 
methode, waarbij we steeds twee punten, die natuurlijk best vlak 
naast elkaar kunnen liggen, door een rechte lijn (of lijntje) met 
elkaar verbinden. Hierdoor hoeven we alleen de coördinaten van 
de uiteinden van de te tekenen lijn te berekenen en bovendien zien 
de tekeningen er beter uit dan met de puntgraphies-methode. 


Bewust hebben we ons tot het gebruik van PLOT en DRAW beperkt 
en ook tot het gebruik van 'zwart op wit'. Mede hierdoor zijn de 
programma's kort en overzichtelijk, zodat de tekentechniek goed 
tot uitdrukking komt. Wij raden u echter aan de programma's uit 
te breiden met meer grafische opdrachten en kleur. In de appendix 
geven we hiervan enkele voorbeelden. Ook geven we in de appen- 
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dix een programma voor het inlezen van een machinetaalroutine, 
waarmee een heel vlak gekleurd kan worden. Zo hier en daar in het 
boek vindt u suggesties voor het aanbrengen van kleur in de teke- 
ningen, 


Structuur van een grafisch programma 
In alle programma's gebruiken we steeds dezelfde variabelen en 
dezelfde programmastructuur. De verschillende variabelen beteke- 


nen: 


xl,‚yl coördinaten van het beginpunt van een lijn 
Xx2,y2 coördinaten van het eindpunt van een lijn 


u,v oorsprong van het wiskundige coördinatenstelsel; 
dikwijls is dit het midden van het beeldscherm (128,88) 
h de waarde 0,5 voor het afronden op helen 
k de vermenigvuldigingsfactor voor functiewaarden 
w‚wl _ hoeken bij trigonometrische functies 
rd het getal 1/180 voor het omrekenen van graden in radialen 


Een grof structuurdiagram voor de programma's is hieronder weer- 
gegeven. 





L begin programma 





beeldscherm schoonmaken 





variabelen u‚v‚h‚k,rd,enz, vastleggen 





punt Pl(x1,y1) vastleggen 





doe zolang nodig 





nieuw punt P2(x2,y2) berekenen 





verbindt Pl met P2 





punt P2 wordt punt Pl (xi=x2 
yl=y?) 














einde programma 





4 ZX Spectrum graphics 


In alle programma's gebruiken we twee belangrijke transformatie 
formules. 


Voor het grafische coördinatenstelsel van de ZX Spectrum is de 
linkerbenedenhoek van het scherm de oorsprong. Wijzelf gebruiken 
echter bijna altijd een coördinatenstelsel met de oorsprong in het 
midden van het scherm. Om de beeldschermcoördinaten x',y' (zoals 
de ZX Spectrum ze gebruikt) te berekenen uit de coördinaten x,y 
zoals wij ze gebruiken hanteren we de transformatieformules: 


x'=INT(utxth) : y' = INT(v+y+h) 


Bekijk de figuur hieronder eens. Duidelijk is dat elk punt p(x,y) in 
ons coördinatenstelsel (gestippeld assenkruis) door bovenstaande 
formules wordt getransformeerd (overgebracht) naar een punt 
p(x',y') in het coördinatenstelsel van de ZX Spectrum. Hierbij hoeft 
onze oorsprong (u,v) natuurlijk niet altijd precies in het midden 
van het beeld te liggen. 








Genoeg theorie, nu de programma’s. De ideeën voor de programma's 
hebben wij gekregen bij het doorbladeren van Amerikaanse, Duitse 
en Franse computertijdschriften. Alle programma's zijn echter eigen 
creaties of bewerkingen en vereenvoudigingen van reeds gepubli- 
ceerde programma's. 
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tekent een 'diagonaalweb'. Elk van de acht bovenste 


punten wordt door een rechte lijn verbonden met elk van de onder- 


ste acht punten. 





Programma 1 


S 
Kd 
a 
28 Ê 
we 
Nd a 
hal: AT Lo 
Hij 
aw z 
bdd uw 
Lob Ë 
WO u 
er 
Le aad z 
Ore LJ 
| es hd 
© sao > 
tou ww 
sHas ex 
>e OEL Z 
EJA XD 
n-KOAOWK O 
„WOL  Ziiuk 
Ok Zw 
SEESSSESGES 
AND HNO ONS 


added 
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Programma 2 is een fraai voorbeeld van het Moiree-effect. Geniet 
ervan, 


189 REM prooranna 2 
oiree effect 


118 25 
120 FOR j=0 TO 175 STEP 7 

158 PLoY 40,175 

Ìäg DRAW },-178 

153 PLOT 48,175 

188 DRAW 178,j=175 

170 NEXT J 

1589 FOR j=0 TO 175 STEP 7 

199 PLoï 215,0 

289 DRAW 178, J 

219 PLOT 215, 

229 DRAW -17849,175 

230 NEXT j 

240 IF INKEY&a"" THEN GO TO 240 
250 SToP 
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Programma 3 tekent een reeks driehoeken in perspectief. Het 'cen- 
trum van vermenigvuldiging' heeft de coördinaten 5(0,69). De drie 
hoekpunten van de driehoek die vermenigvuldigd wordt zijn 
P1(6,81), P2(20,78) en P3(12,63). Alle coördinaten zijn beeld- 
schermcoördinaten (oorsprong links onderaan). De drie richtingen 
van vermenigvuldiging (richtingsvectoren) zijn SP;, SP, en SPs. 
Deze worden in het programma als a(1);b(1), a(2);b(2) en 
a(3);b(3) vastgelegd. De vermenigvuldigingsfactor k is de lus- 
variabele van de buitenste FOR-NEXT-lus. k loopt van 0 tot 9 in 
stappen van 0,4. Hieronder zien we de uitgangssituatie, het resul- 
taat van het programma en het programma zelf. 











y 
7 
fe 
# 
1 
/ 
4 „7 
sid PÁ ze 
É, p 
184 p2 
Le 
69de” 
SI = 
63d is 
Ps, PN 
Ee Á) Wi 
__ 

(0,0) 6 12 20 X 
188 REM programma 3 driehoeken 
in Perspectief 

112 CLS 

128 DIM x(3) : DIM y (3) 

138 DIM ai3) : DIM b (3) 

148 FOR jel TO 3 

152 READ aij),bij) 

169 NEXT j 

172 LET xö0 ‚ LET 48-39 

188 FOR ka@& TO 9 STEP d,4 
198 FOR jel TO 3 

202 LET xj sxdekaa lj) 

21e LET gij eyBekeb ij) 

228 NEXT 

238 PLOT XI) ,y (1) 

248 DRAW XE -x (1) ,y () =y (1) 
259 PLOT x{i2),y ie) 

288 DRAW XI -x (2) ‚U LI) =y (2) 
2/8 PLOT X(3),y (3) 

28e pan, xii (3) ‚y (1) -y (3) 


E90 NEXT Kk 
328 IF INKEY$z"" THEN GO TO 909 
328 DATA 6,12,22,9,12,-6 
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Programma 4 tekent een reeks ingeschreven regelmatige zeshoeken. 
Behalve bij de eerste zeshoek, liggen alle zes de hoeken van elke 
zeshoek precies in het midden van een zijde van de omgeschreven 
zeshoek. De coördinaten van de hoekpunten van de eerste (buiten- 
ste) zeshoek worden met trigoniometrische functies (sinus- en 
eosinusfunctie) berekend. De hoek die hierbij als argument van de 
functies gebruikt wordt doorloopt de waarden 60°, 120°, 180°, 2409, 
300°, 360° en 420°. De middelpunten van de zijden van een zeshoek 
worden aangegeven met a(1);b(1) t/m a(6);b(6). Dit worden de 
hoekpunten van de volgende zeshoek. 


198 REM Programma 4 
ing@schreven zeshoeken 

119 CLS 

180 OEM x(2) : DIM u {23 

138 OIM at7) : Dim bl? 

180 LET u=ize : LET voos 

150 CET r=88 : LET h=0.5 

189 CET w=b0xp 17160 

18 FÖR jel 10 5 

188 LET wisjau 

18e LET xt sinT werscos wan 

+ 

200 LET VLDEINT Wwerssin Wa) 

Ed 

B10 NEXT J 

20 FOR n=1 TO 20 

230 FOR jz1 TO 6 

248 Plot xijd,8 u 

Pir hd AU Ox ijtii XG ‚U OSI ey 

260 NEXT 3 

279 FOR ket TOG 

280 LET ak) INT ((X (KJ +X (Kk +1 

1) Z8+h) 

292 LET BIKJZINT (Cy ik) +y (k+1l 

)) 28th) 

308 NEXT k 

210 LET atzsatu : LET oi) eb 

329 FOR jat TO? 

„2e LET kijlmaln s LET vise 

3Â0 NEXT J 

350 NEXT n 

360 IF INKEY&s"" THEN GO TO 960 

5% SToP 





Zie de appendix voor een speciaal effect met dit programma, waarin 
gebruik gemaakt wordt van een machinetaalprogramma voor het 
kleuren van vlakken. 
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Programma 5 tekent alle diagonalen in een n-hoek. De n hoekpun- 
ten liggen op de omtrek van een ellips of van een cirkel afhankelijk 
van de waarden die we in het begin van het programma voor de 
halve lange as a en halve korte as b kiezen. Voor de berekening 
van de coördinaten van de punten op de omtrek van de ellips ge- 
bruiken we niet de (cartesische) vergelijking (x?/a?) + (y?/b?) = 1, 
maar de parametervorm x = a.cospen y = b.sinw. 


Als de parameter p loopt van 0° tot 360°, dan beschrijft het daar- 
aan toegevoegde punt P(x,y), met x =acospen y = bsinp, precies 
de omtrek van een ellips. Hieronder zien we een bepaald punt 
P(x1,y1) op de omtrek van de ellips met de daarbij horende waarde 
van de parameter wj. Voor dat punt geldt: x1 = acosp: en 

Yi = bsinp:. 


y 


NR (x,y) 





P- 





hi 
ij 
ï 
XI x 








Als xj = aCOSPj en yi =bsinp, dan geldt ook 
x? =a?cos?pi en y? = b?sin? wr en ook 
2 2 
x ar 
ar =eos?i en 5 = sin?p;, dus dan is 
2 2 
xi Dn B 
az 5 = cos? + sin? 


2 
en omdat cos?p, + sin? pi gelijk is aan 1 geldt 55 + Er =1 


en hebben we de cartesische relatie tussen xj en ya afgeleid uit de 
parametervorm. 


Voor het berekenen van de coördinaten van het j-de hoekpunt 
(YH) verdelen we de maximale middelpuntshoek van 360° in 
n gelijke hoeken van elk 360/n graden. 
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Als a de lengte van de halve lange as is èn 
b de lengte van de halve korteas èn 
n het aantal hoekpunten van de n-hoek èn 
w gelijk is aan 360/n èn 


wl de hoek die hoort bij het j-de hoekpunt èn 

x(j) de x-coördinaat van het j-de hoekpunt t.o.v het 
middelpunt van de ellips èn 

y{) de y-coördinaat van het j-de hoekpunt t.o.v, het 
middelpunt van de ellips 


dan zou je in een BASIC-programma de coördinaten x(j) en y(j) als 
volgt kunnen berekenen: 


w_ =360/n : wl = (j-1)tw 
K()  = INT (axCO8(w1)) 
YÚ) _=INT(b*SIN(wi)) 


Als we dit doen maken we twee fouten. De eerste fout is dat de 
hoek wl, als argument van de COSinus- en de SINusfunctie, in 
graden is uitgedrukt, terwijl BASIC vereist dat het argument van 
deze functies in radialen wordt uitgedrukt. De tweede fout is dat 
het grafische coördinatenstelsel van de ZX Spectrum de oorsprong 
van het coördinatenstelsel voor het scherm in de linkerbenedenhoek 
van het scherm verwacht en niet in het middelpunt van de ellips. 
Wat we dus nog moeten doen is: 

a. bereken wl in radialen èn 

b. transformeer x(j) en y(j) naar beeldschermcoördinaten. 


We weten dat een hoek van 360° overeenkomt met 21 radialen, 
waarin n = 3,14159..... . Dit betekent dat een hoek van 1 graad 
overeenkomt met een hoek van 21/360 of n/180 radialen. De trans- 
formatieformules voor de transformatie van 'onze' coördinaten naar 
beeldschermcoördinaten hebben we op p.4 gegeven. De drie BASIC- 
opdrachten voor het berekenen van de beeldschermcoördinaten van 
het j-de hoekpunt van de n-hoek worden nu: 


w _= (360/n)* PI/180 : wl = (j-1)4w 
x(j) = INT (u + axCOS(w1) + h) 
YO) =INT(v + bxSIN (wl) + h) 


We gebruiken de BASIC-functie PI, 


In het navolgende (en ook in het voorgaande) programma komen we 
bovenstaande opdrachten tegen. Kiezen we voor a en b dezelfde 
waarde dan krijgen we een regelmatige n-hoek met al zijn diagona- 
len. Bovendien is de ellips dan een cirkel geworden. 


ZX Spectrum graphics 












SA 
SER 


183 REM programma 5 
diagonalen in een n-hoek 
118 CLS 
Je INPUT "HALVE GROTE AS A{ala 


ia 
Aer “HALVE KLEINE AS Bee8 


re 
& 
6 


140 INPUT “HOEVEEL HOEKPUNTEN N 
< zin 
158 CLS 
168 DIM Xin) : DIM yin) 
172 LET vel28 : LET v=88 
182 LET had,5 
192 LET w=i362/n) #PI/180 
288 FOR j=i TO n 
218 LET wiaij-l) w 
ae LET X{jJ=INT (U+asC0S (wid 
+ 
Pal LET yIjJzZINT (vebexSIN (wi) 
242 NEXT j 
258 FOR izl TO n= 
Ze FOR dik Ton 
28 PLOT XLid yi) 
288 DRAW xj lid, =y Cid 
290 NEXT j 
382 NEXT i 
318 IF INKEY&e"" THEN GO TO 318 
328 STOP 
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Programma 6 is wat veeleisender in die zin dat het wat voorberei 
ding nodig heeft. 


We willen een reeks ingeschreven vierkanten zo tekenen dat de 
hoekpunten van een ingeschreven vierkant uit de reeks op de zij- 
den liggen van zijn voorganger. We zien deze situatie voor het eer- 
ste, tweede en derde vierkant van de reeks in de onderstaande 
tekening. De afstand tussen een hoekpunt van een vierkant en een 
hoekpunt van het volgende vierkant (S: en Sz) is steeds een vast 
deel van de zijde waarop de hoekpunten liggen. Zo is S, = a;/k en 
S2 = ar/k. In het algemeen is Sn = an/k, waarbij an de zijde is van 
het n-de vierkant in de reeks en k de verkleiningsfactor. De waar 
de voor k kunnen we zelf kiezen. k=16 geeft een fraaie tekening. 


XLI, YY) 


X(4),Y4) a, | XY) 








XY) a, X(2), Y(2} 


XD, YY) 


CLC); (DY (2)D); (X(3),Y(B)) en (X(4),Y(4)) zijn de 
coördinaten van de hoekpunten van een bepaald vierkant uit de 
reeks. Hoe berekenen we nu de coördinaten van de hoekpunten van 
het volgende ingeschreven vierkant ((XX(1),YY(1)); (XX(2),...)? 


Voor het eerste hoekpunt van het eerste ingeschreven vierkant geldt: 


KX =x EDEN an 


vv =va + OD 


Ga na dat dit klopt in bovenstaande tekening. 
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Voor het tweede hoekpunt van het ingeschreven vierkant geldt: 


vvo + OO 
In het algemeen geldt: 
xx = xD èn 
vra = va « DO 
voor J = 1,2,3,4, waarbij X(5) =X(1) en 
Y(5) = Ya). 


Als we het tweede vierkant getekend hebben, veranderen we 
XX(I) in X(I), YY(I) in Y(1), XX(2) in X(2). YY(2) in Y(2), 
enzovoorts, en tenslotte X(5) in X(I) en Y(5) in Y(1) en we gaan 
met dezelfde formules opnieuw XX(1),YY(1),... enz. berekenen, 
maar dan voor de hoekpunten van het volgende vierkant. We zien 
dit gebeuren in de regels 290 t/m 340 van het volgende programma. 
De regels 210 t/m 240 tekenen het vierkant, terwijl de regels 250 
t/m 280 de hoekpunten van het volgende vierkant berekenen, Een 
structuurdiagram voor het programma ziet er zo uit: 





begin programma ingeschreven vierkanten 





scherm schoonmaken 





lees waarde voor K in 





bepaal coördinaten voor het eerste vierkant 





voor het eerste t/m 40-ste vierkant doe 





| teken dit vierkant 








bepaal coördinaten voor het volgende vierkant 





einde programma 
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Hier komt het programma: 


100 REM programma & 
ingeschreven vierkanten 


LS 
„128 INPUT “TOETS K IN 1cKcR@ 


138 LET hz0,5 
DIM x(8) : 


180 FOR jal TO 
17® READ Kij), 
NEXT j 


189 

198 CLS 

290 FÖR nel TO 40 

819 FOR jat TO 4 

280 Plot X(j),175-y(j) 

290 DAAU xt ek CI CH Wd 
+ 

Bad NEXT J 

250 FOR jz TO 4 

280 LET atjaxij+INT CixEj+L 
NX LI) ZK +h] 

278 LBT bij) sytj)+INT (Wy tJel 
ey O1 ZK th) 

288 NEXT j 

390 FOR ja TO 4 

309 CET xijmalj 

319 LET V(j abi 

389 NEXT | 

330 LET xi5) ex (1) 

340 LET yB) ay (1) 

352 NEXT n 

360 IF INKEY&u"" THEN GO TO 360 
379 STOP 

388 DATA 42,175,215,175,215,0 
389 DATA 40,5,40,178 











U kunt dit programma als uitgangspunt voor een uitgebreider pro- 
gramma gebruiken. Verdeel het beeldscherm in een aantal vierkan- 
ten. In elk van deze vierkanten tekent u, met bovenstaand pro- 
gramma, een reeks van ingeschreven vierkanten. Teken steeds 
identieke reeksen, of kies steeds een andere waarde voor k, of 
probeer ze ten opzichte van elkaar te laten draaien. 





2 Grafieken van functies in 
cartesische vorm 


In het eerste hoofdstuk hebben we alleen rechtlijnige patronen gete- 
kend, bijvoorbeeld het programma voor het tekenen van alle diago- 
nalen in een regelmatige n-hoek. 


In dit en het volgende hoofdstuk zullen we ons uitvoerig bezig hou- 
den met het tekenen van de grafiek van een aantal continue en niet- 
continue functies. In de wiskunde noemen we de grafiek van een 
niet-lineaire functie een kromme. De grafiek van een lineaire func- 
tie (bijvoorbeeld de functie y = 4x + 3) noemen we een rechte. We 
kunnen de vergelijking van zo'n kromme op drie manieren formule- 
ren: 





A ; met cartesische coördinaten y=f(x) 
B ; met poolcoördinaten r=f(@) 
C ; of in parametervorm :xe=f(t), y =g(t) 


Niet-wiskundigen kennen vaak alleen de gewone (cartesische) vorm 
y = f(x). Een voorbeeld van een niet-lineaire functie in cartesische 
vorm is de functie y = 2x? - 3x + 4. De grafiek van deze functie 

is een parabool. 


De functie r = 110.eos(4p) stelt een kromme in poolcoördinaten 
voor. Als we de hoek w laten lopen van 1 tot 360 graden en we 
tekenen bij elke hoek p onder die hoek een punt op afstand r van 
de oorsprong, dan krijgen we een bloemfiguur als op p.33 staat. 
Deze vorm r = 110.cos(4p) heet de poolcoördinatenvorm. Een punt 
in het platte vlak wordt nu niet gekenmerkt door de afstand van 
dat punt tot de x- en y-as (cartesisch), maar door de afstand tot 
de oorsprong (r) en de hoek wp die de x-as maakt met de lijn die dat 
punt met de oorsprong verbindt. Poolcoördinaten komen in hoofd- 
stuk 3 aan de orde (zie tekening op p.16). 


Er is nog een manier om de vergelijking van een kromme weer te 
geven en dat is de parametervorm. Hierbij worden de x- en y- 
coördinaat van een punt op de kromme afhankelijk gemaakt van een 
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Pe Pi ,ya) P(r,o) 





+ 

| 

| 

Î 

Î 

i 
ed 
oorsprong xr 


X 





Hetzelfde punt P in een 
cartesisch coördinatenstelsel en een poolcoördinatenstelsel 


bepaalde parameter. Een voorbeeld zijn de vergelijkingen 
Xx =a.cos(t) en y =a.sin(t). 


Als we t laten lopen van 0 tot 360°, beschrijven de daarbij horende 

punten (x,y) precies de omtrek van een cirkel met straal a. De car— 
tesische vorm van deze cirkelvergelijking is x? + y? =a?, die velen 

direct zullen herkennen als de 'cirkelvergelijking'. 


In dit hoofdstuk houden we ons dus alleen bezig met het tekenen 
van grafieken van functies en relaties die door middel van cartesi- 
sche coördinaten beschreven worden. De andere twee vormen komen 
in hoofdstuk 3 aan de orde. 


Continue functies 


Een continue functie is een functie waarvan de grafiek in één vloei- 
ende beweging van ons 'potlood', dat wil zeggen zonder het potlood 
van het papier te hoeven halen, getekend kan worden. Er mogen 
‘knikken! in voorkomen maar geen onderbrekingen (zie tekening op 
p.11). 


De linkergrafiek kunnen we in één beweging, zonder het potlood 
van het papier te halen, tekenen; bij de rechter grafiek kan dat 
niet. De linker grafiek is de grafiek van een continue functie; de 
rechter van een niet-continue functie. Als we alleen links of alleen 
rechts van het punt xo kijken dan is de functie op die intervallen 
natuurlijk wel continu! De functie is alleen niet-continu in het punt 
X0. 
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Y y 


XX 





_ 
X, EN %o 


1 
ij 
Ï 
Ô 
| 
1 
| 
\ 
| 
ï 


Knikpunten 


’ 
Functie is niet-continu 
in het punt Xa 


We willen nu een algemeen programma ontwikkelen dat, gegeven de 
grenzen a en b van een bepaald interval (asx=b) de grafiek tekent 
van een willekeurige continue functie y = f(x). Mochten de x-as en 
de y-as in het gebied liggen waarin de grafiek van de functie wordt 
getekend, dan willen we dat ook beide assen door het programma 
getekend worden. Om het programma kort, maar toch algemeen, te 
houden brengen we de volgende vereenvoudigingen aan. 


1. De vergelijkingen van de functie, die getekend moet worden, 
wordt in een subroutine, vanaf regel 1000, opgenomen. De 
functiewaarde y = f(x), bij een bepaalde x-waarde, wordt 
berekend door op dat moment met de gewenste x-waarde de 
subroutine aan te roepen (GO SUB 1000). Er wordt dus niet 
gebruik gemaakt van de methode waarin de functie als tekst 
(INPUT-opdracht met stringvariabele) wordt ingelezen waarna 
het programma deze tekst zelf omzet in een goede BASIC 
functiedefinitie. Een andere mogelijkheid zou zijn om de functie 
met een DEF FN-opdracht in het programma te definiëren. 


2. Op de coördinaatassen wordt niet automatisch een schaalverde- 
ling aangebracht en ook wordt er geen tekst bij afgedrukt. 
Veel micro's met HRG*-mogelijkheden kunnen niet tegelijkertijd 
tekenen en 'schrijven'. Wilt u toch graag weten wat bij een 
bepaalde waarde van x de functiewaarde (y) is, dan kunt u 
bijvoorbeeld na het tekenen van de grafiek de computer (op 
de printer) een lijstje met x- en y-waarden laten afdrukken. 


Terug naar het programma-ontwerp. In het eerste deel van het pro- 
gramma berekent de computer, na het inlezen van de waarden van 
de intervalgrenzen a en b, de grootste en de kleinste functiewaarde 
(y-waarde) in het opgegeven interval [a,b]. De waarde van hp 


*HRG = High Resolution Graphics 
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(hoogste punt) en Ip (laagste punt) worden op het beeldscherm 
afgedrukt. We weten zo in welk gebied de computer gaat tekenen. 
Hierna vraagt het programma of we hp nog groter en of we Ip nog 
kleiner willen maken om daarmee een fraaiere tekening te krijgen. 
Zo niet, dan toetsen we voor hp en Ip dezelfde waarden in als het 
programma heeft berekend; anders toetsen we de door ons gewenste 
maximale en minimale functiewaarden in. 


Het volgende programmadeel (regels 290-380 op p.20) bevat de lus 
(FOR x=a TO b STEP dx) voor het tekenen van de grafiek. De 
coördinaten x,y van een punt op de grafiek moeten met de juiste 
transformatieformules omgezet worden in beeldschermcoördinaten. 


Om het interval [a,b] voor te kiezen x-waarden (asx<b) om te 
zetten in het HRG-interval [0,255] (Osx's255) gebruiken we de 
volgende evenredigheid: 


Xa 


—a) : (b-a) =X!: Ge 
(x-a) : (b-a) =x' : 255 ofwel Da 


Xx! 
“255 





in BASIC: 
XX =INT(kx*(x-a)th), met kx = 255/(b-a). 
We doen dit ook voor het afbeelden van het functiebereik Ipsyshp 
op het grafische beeldbereik Osy's175: 
EERE Le -N 
(y-1p) : Chp-lp) =y' : 175 ofwel bop = 155 
Dit geeft in BASIC: 
yy = INT (ky#(y-Ip)+h) met ky = 175/(hp-Ip). 
We zien deze berekeningen van xx en yy in de regels 310 en 320. 


In de volgende regels worden twee, op deze wijze berekende, pun- 
ten (x1,yl) en (x2,y2) door een recht lijntje met elkaar verbonden, 


Vervolgens (regels 390-420) wordt de ligging van de x- en y-as 
bepaald. Kiezen we in bovenstaande transformatievergelijkingen 
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voor x en y de waarde nul, dan vinden we respectievelijk de ligging 
van de y-as en van de x-as (regels 390 en 420). 


Met deze uitleg is de werking van het programma hopelijk te volgen. 


Test u het programma 7 met willekeurige, maar continue, functies. 
In dit voorbeeldprogramma hebben we de functie y = e” } IX. cos(x) 
gekozen, hetgeen de grafiek van een gedempte trilling oplevert. 
Hieronder zijn enkele ideeën voor minder moeilijke functies. 


functie regel 1000 waarde voor a waarde voor b 
y = sin(x) y= SING) 0 271 (= 6,2832) 
yex? y= Xi 4 +4 
yet y =EXP(x) -3 +3 
yEXI-2KE- KX y= REKER DEKT -1 +3 











0, 1x 


y= e .cos(x) 


(zie regel 1000 van programma 7) 
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128 REN programma 7 rafiek 
een continue functie 


ëe3 
ELN ZNOUT “LINKER-INTERUAL GREN 
geen tNpur “RECHTER-INTERVAL GRE 
Nao re Ad THEN LET Caa: LET au 


ì LET 
152 Let “Áp=-190000: LET \pe1328 





sa 

158 LET dxaib-a) 7/64 

172 FOR xaa TO b STEP dx 

189 GO SUB 450 

1838 IF yohp THEN LET hpay 

290 IF ĳ<lp THEN LET Lpay 

818 NEXT x 

;Ere PRINT "GROOTSTE Y-WARRDE 

‚230 PRINT "KLEINSTE Y-WAARDE: " 
P 

‚#o INPUT "BOVENGRENS VOOR Y " 

„25e INPUT "ONDERGRENS VOOR Y " 


CL$ 
Ee LET Kyal7 
Pz 
LET ha@.5 

EOR Xzà TO b STEP dx 

GO SUB 1829 

LET XxaINT (KX#iX-a) bn 
LET HUE ete +h 
TETN T XLaxx: ler 


0 Ta se 
LET, cri LET yauyy 


y 
DRAW X@-XL,YR-yL 
LET Xluxä: dÉr yiey@ 


NEXT Xx 
LET x120: LET YIeINT (kyéí- 


OOPPROPPCR PPUGWGUGEEGOGGKDN DD. 
PONDUIRGE DX PE OB IROLPGOPOGE IGT AV 


GELBOS EÉUEFEOHESEELE U EELEESTES 


IS h) 

LET X2=255, EE 8 Kú 

Ir Varg Aje yic25S THEN PLO 
T xloyl: DR BXL, yd=yl 

geer Kie INT IKx&l-ad th): LET 

LET Xeuxl: LET van175 

IF X1>58 AND x142255 THEN P 
L XL,yl: DRAW XR-XL,YB=yl 

Eid INKEY&a"" THEN GO TO 458 

STOP 
1 LET ried (Blanes Oo) 
El RETURN 
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Nu volgen drie korte programma's. 


Programma 8 is de eerste van deze drie. Dit programma tekent tien 
in fase verschoven sinuskrommen, allemaal in hetzelfde coördinaten- 
stelsel. De vergelijking van het stelsel is 


y =sin(x*np), met als fase p = 5 (209), 


We vinden de tien vergelijkingen door voor n de waarden 0,1,2 t/m 
9 te kiezen. De lusvariabele j (regel 150) komt overeen met de 
beeldschermcoördinaat xx. Omdat we in dit programma de beeld- 
schermcoördinaat xx, dat wil zeggen j, in stappen van 5 naar 255 
laten lopen (regel 150), moeten we eerst de bij j behorende waarde 
voor x berekenen. Dan pas kunnen we de functiewaarde y = f(x) 
berekenen. We kunnen de evenredigheid 


Gera) : (bra) =j : 255 


natuurlijk ook gebruiken om x uit j te berekenen; hieruit volgt: 


(b-a) 


tk 


In het onderstaande programma zijn de gekozen intervalgrenzen 
a=0enb= 2m, zodat bovenstaande vergelijking wordt: 


Zn 
255 ” 





Xx 


Als we c = Ze nemen, wordt in het programma het argument van de 


sinusfunctie dus jxctnsp. 


In het programma wordt de waarde jxc als x berekend, dus zien we 
in regel 170 x+nxp als argument van de SINusfunctie. Natuurlijk 
moeten we de functiewaarde y = SIN(x+n«p) nog omzetten naar de 
beeldschermcoördinaat y,‚ zodat we krijgen 


y= INTCv*kaSIN (x+nep)+h). 


Door voor de schaaleonstante k de waarde 87 te kiezen (regel 120) 
krijgen we een mooie 'grote' tekening. 


Wie met kleuren wil werken kan dit programma uitbreiden door de 
diverse krommen andere kleuren te geven. Iets moeilijker zal het 
zijn de 'banen' tussen de sinuskrommen in te kleuren. Probeer het 
eens. 
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100 REM googranna 8 
Sinuskrömmen 
112 CLS 
180 LET ve87: LET Ke87: LET hed 
5: LET p=Pi/3 
132 LET c=È4P1/255 
142 FOR n= TO 
15@ FOR j=o TO 255 sTEP 5 
189 LET x=j4c 
‚ie CET VAINT WekeSin Gense 
+ 
182 IF jad THEN LET xiaj: LET 
gigs: 09 15 230 
80 LET x@ej: LET v2ey 
220 PLOT X1,176-v1 
210 DRAW xa-xl,vi-va 
228 LET xlexa: LET yieya 
230 NEXT J 
249 NEXT n 
25Q IF INKEY&="" THEN GO TO 280 
250 STOP 


De illustraties zoals hierboven en op de pagina hiernaast zijn door 
de auteur gemaakt op een plotter, waarbij de linkerbovenhoek de 
oorsprong is en niet zoals bij de ZX Spectrum de linkerbeneden- 
hoek. Om hetzelfde effect te krijgen als de illustraties aangeven 
moeten we in sommige programma's ook uitgaan van een oorsprong 
linksboven. Kijk welke figuur u krijgt als u in regel 170 niet 
v-k#SIN maar v+kxSIN neemt en in regel 200 niet PLOT x1,175-y1 
maar PLOT x1,ylen in regel 210 niet DRAW x2-x1,yl-y2 maar 
DRAW x2-x1,y2-y1. 


ZX Spectrum kleurcodes 














Code Kleur Code Kleur 
0 zwart 4 groen 
1 blauw 5 vaalblauw 
2 rood 6 geel 
3 paars ni wit 
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Programma 9 tekent een stelsel parabolen met de vergelijking 
y=-tx? +t 


Alle parabolen (bepaalde t-waarden) hebben dezelfde snijpunten 
met de x-as (x = len x=-1). 





198 REM progr dan, ta 8 
ParaboolstelseL 

118 CLS 

eer LET vel2ë: LET va87: LET he 

138 km-87 TO 87 STEP 18 

148 "RR xe-118 TÔ dn STEP 5 

158 LET XX=INT LU+X+h) 

160. LET varken 0400 ek: LET 

VeiNT Ate h) 

179 ket ee LET XLaxx: 

LET vac; 8 TO 

188 T err: atËr Yan 

hid PLOT xa,yi: DRAÙ xe-xi pe 

Zêo LET X1ax@: LET yiey@ 

21% NEXT x 

22Ô NEXT K 

230 IF INKEYSn"" THEN GO TO 238 

248 STOP 


Voeg toe 145 IF x <= O THEN INK 1 : GOTO 150 
146 INK 6 


of 145 IF INT((k/10)/2)*2=k/10 THEN INK 1 : GOTO 150 
146 INK 6 


om wat kleur in de tekening aan te brengen. 


In dit programma maakt het niet uit of de oorsprong linksboven of 
linksonder ligt. Het plaatje is symmetrisch rond de horizontale lijn 
midden in het scherm. We kunnen in regel 160 y=INT(v-y*h) nemen 
en in regel 190 PLOT x1,175-yl: DRAW x2-x1,yl-y?; dit geeft 
dezelfde figuur, die alleen van bovenaf getekend wordt in plaats 
van onderaf. 


Soms willen we de oppervlakte tussen de grafiek van een functie 
en de x-as berekenen (de integraal bepalen) of laten tekenen, 
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Programma 10 kleurt zo'n oppervlak tussen de x-as en de grafiek 
van de functie 





y = cos») — sos(3e) nf EE) her Sein . 


Ook nu is de lusvariabele j de grafische coördinaat xx. De waarde 


voor a is -mn en die voor b is tr, waaruit volgt dat (zie p.21) 


j.2n 
255 We 


Ook nu nemen we in het programma c = Ze 


x= 








SOEP pREEPERLEDE 
Pe SHPSOUIEGRODE 
S_SSEESELESSEEEE 


ror 





Probeer de regels 160, 170 en 180 zo te veranderen, dat u een 
figuur krijgt die gespiegeld is ten opzichte van de horizontale lijn 
in het midden van het scherm (kijk naar programma 8). 


Probeer dit programma zo te maken dat de drie oppervlakken in 
bijvoorbeeld rood, wit en blauw gekleurd worden. 
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Niet-overal-continue functies 


De functie y = Rn kan niet met behulp van programma 7 gete- 
kend worden. Zouden we bijvoorbeeld voor x het interval -5s X: 5 
kiezen, dan liggen hierin de punten x = 3 en x =-2. Dit zijn de 
twee punten waarvoor de noemer van bovenstaande functie nul is. 
Als de computer bij het tekenen bij één van deze twee punten 
belandt, zal op het scherm een foutmelding (DIVISION BY ZERO) 
‘delen door nul' verschijnen en zal het programma afgebroken wor- 


den, 





Dit zal ook gebeuren bij logaritmische functies met een niet-positief 
argument of bij wortelfuncties met een negatief argument. Zelfs bij 
het tekenen van continue (nette) functies kunnen problemen optre- 
den. We komen hierbij in de problemen als de functiewaarden heel 
groot of heel klein worden. Transformatie van dergelijke waarden 
naar beeldschermcoördinaten (0-255 en 0-175) heeft dan geen 
enkele zin meer, omdat de aard van het verloop van de functie in 
het geheel niet meer tot uitdrukking komt. Het is daarom beter om 
de functiewaarden van een continue functie naar boven en beneden 
te begrenzen. Dit heeft tot gevolg dat de grafiek van een continue 
functie er uit zou kunnen zien als de grafiek van een niet-continue 
functie. 


Het zal duidelijk zijn dat je een algemeen programma voor het teke- 
nen van een willekeurige continue of niet-continue functie niet 
zomaar even opschrijft. Dergelijke programma's kom je in de vak- 
literatuur dan ook niet of nauwelijks tegen, en mocht je er wel een 
tegenkomen, dan betreft het òf een heel groot programma òf een 
programma dat voor een bepaalde functie, waarvan men van te 
voren weet waar de discontinuiteiten zitten, geschreven is. 


We geven nu een verbazend kort programma voor het tekenen van 
de grafiek van een willekeurige functie y = f(x). Veel wiskunde- 
leraren, scholieren en studenten zullen nu hun ‘oren! spitsen. U 
begrijpt dat. gezien het bovenstaande, hierbij wel enkele beperkin- 
gen gelden: 


E: Degene die het programma gaat gebruiken moet een beetje kun- 
nen programmeren. De functie moet namelijk in gedeelten in 
een subroutine (vanaf regel 1000) beschreven worden. 

2, De programmagebruiker moet voldoende wiskundige kennis 
bezitten om te kunnen bepalen voor welke x-waarden functies 
moeilijkheden kunnen opleveren. 


We zullen zien dat, normaal gesproken, slechts een paar BASIC- 
regels nodig zijn om de functiebeschrijving te programmeren. Het 
berekenen van nulpunten met behulp van een of ander numeriek- 
wiskundig algoritme is in het geheel niet nodig. 
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In het onderstaande programma 11 gebruiken we twee variabelen 
fz en fa als zogeheten vlaggen (Flags). Een vlag is een variabele 
die slechts twee waarden (vaak 0 en 1) kan aannemen. 


Om de werking van de vlak fz duidelijk te maken geven we hieron- 
der de subroutine 1000 met de functiebeschrijving van de functie 


Kie 
x7-x-6 


1000 LET n=x*x=x-6 

1010 IF n=0 THEN LET fz=1: RETURN 

1020 LET y=Cx*xt3)/n 

1030 IF y<lp OR y>hp THEN LET fz=1: RETURN 
1040 LET fz=0: RETURN 


De vlak fz wordt alleen op 1 gezet als het punt (x,y) niet op het 
scherm getekend kan worden. Dit is het geval als de functie niet- 
continu is in het punt (x,y) (n = 0; x = 3 en x =-2) of als de 
functiewaarde buiten het opgegeven functiebereik (lp<y<hp) valt. 


Waarvoor dient nu de tweede vlag fa? Bekijk het volgende stukje 
programma eens: 


200 LET fa=1 

210 FOR x=a TO b STEP dx 

220 LET x2=INT(kx*(x-a)th): GO SUB 1000 
230 IF fz=1 THEN LET fa=1: GO TO 310 
240 IF fa=1 THEN GO TO 290 

250 LET y2=INT (ky* (y-Lp)+h) 

260 _ PLOT x1,y1 

270 DRAW x2-x1,y2-y1 

280 LET xl=x2: LET yl=y2: GO TO 310 
290 LET x1=: LET y1=INT (ky# (y-Lpdth) 
300 LET fa=0 

310 NEXT x 





Als fz=l dan wordt fa ook gelijk aan 1 gemaakt en wordt de volgen- 
de x-waarde bekeken (FOR-lus) zonder dat het nieuw berekende 
punt met het laatstgetekende punt verbonden wordt. Is fz echter 

0 (berekende punt ligt in het beeldvlak) dan wordt onderzocht of 
fa daarvoor soms op 1 gezet is. Is dit zo (tweede THEN) dan moet 
het nieuw berekende punt als nieuw beginpunt (x1,y1) gekozen 
worden en dit mag dan ook niet met het vorige getekende punt ver- 
bonden worden; fa wordt in dit geval nul gemaakt. 





Niet-overal continue functies 27 


Dus als zowel fz als fa nul zijn, wordt het nieuw berekende punt 
verbonden met het laatstgetekende punt. 


Om te zorgen dat het programma probleemloos werkt moet voor het 
berekenen van het eerste punt van de grafiek, dus voor de FOR- 
lus, de vlag fa op 1 gezet worden, waardoor we er zeker van zijn 
dat de waarden xl en yl berekend worden. Hier volgt het volledige 
programma: 


180, KE piogranan 11 grafiek 
villekeurige functie 


Ùs* 
128 INPUT "LINKERGRENS VOOR X 
138 INPUT "RECHTERGRENS VOOR X 
niee INPUT "BOVENGRENS VOOR Y " 
‚Ì5o INPUT "ONDERGRENS VOOR Y " 
16: ee THEN LET Caa: LET an 
EE Kx TBES/ CDE LET kye17 
P=-lp ha,6 
E Herth BE 
T fa 
IR xea tro b STEP dx 
Blas eeNt xe KAD +): G 
fzal THEN LET farsi: GO 
fasl THEN GO Tö 292 
T 
hi 


ME 


L 
LE 
CL 
LE 
Fo 


0 
Ld 


18 
YRRINT (Ky&(y=-lp) +h) 


zy 
XaeXl, vgev 
xlexa: der” Java: Go T 


dex@: LET VA=INT (kue 


GOSELSCHEEGETENE 
x 
be 


Li 
1 
IF 
IF 
LE 
PL 
DR 
Li 


DK -) ba 
GEN OWIDDORD DIEN te 


OPO: DEBIARP GEOPGOG 
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ET 
ET 
h) 
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maren 


@ 
WEEG ALGE E09 


“GR PSSOD IPEA 


20 
Ed 
PHP 

d 


er 


AW RX, 


® 
Fc nen 
4d de 


mi: LET vaal? 
eo Xi4a5E Tren PL 
RAU _x 53: 
Eine PEEN GS To oo 
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eres GEruud 
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pn 
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O9 GOOGOXEE EXEL: EE0-erG: 
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RX ete 

THEN LET fZel: RETUR 
X ax 
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WohPp THEN LET fz 
2: RETURN 
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ac 
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Dit programma tekent de grafiek van de functie 


Ein ee : kies voor a, b, hp en Ip- resp. 


„5, 5, 10en -10 








7 











Wilt u een andere functie tekenen, bijvoorbeeld y = In(x?-2), her- 
schrijf dan subroutine 1000 (regels 1000,,1010 en 1020): 


1000 LET u=x#x-2: IF u<=0 THEN LET fz=1: RETURN 
1010 LET y=LOG (u) 


U krijgt dan deze grafiek: 








Kies voor a, b, hp en Ip de waarden -25, 25, 8 en -5. 
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Hier volgt een hele fraaie: 











Verander de regels 1000 en 1010 in: 


1000 LET n=x-3: IF n=0 THEN LET fz2=1: RETURN 
1010 LET y=3-x+LOG (ABS ((x-1)/0)) 


Kies voor a, b, hp en Ip de waarden -5, 10, 8 en -4. 
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3 Krommen in poolcoördinaten 
en in parametervorm 


In het vorige hoofdstuk hebben we ons beziggehouden met het teke- 
nen van de grafiek van een continue of niet-continue functie, waar- 
van de vergelijking als y = f(x) geschreven kon worden; dus in 
cartesische vorm. Nu gaan we grafieken tekenen van functies waar- 
van de vergelijking in poolcoördinaten geschreven wordt of in de 
zogeheten parametervorm. Dit levert zeer fraaie ‘beelden! op. 


Krommen met poolcoördinaten 

Functies, waarvan de grafiek een gesloten kromme laat zien, zijn 

vaak eenvoudiger met poolcoördinaten te beschrijven dan in carte- 

sische vorm. Als voorbeeld noemen we de 'hartkromme' (kardioide), 

waarvan de vergelijking in poolcoördinaten eenvoudig ìs, namelijk: 
r=k(1 + eoso), 


maar waarvan de cartesische vorm nogal ingewikkeld is, namelijk 


Gleyt-ka)t _, 
KEER) 





Om de volgende programma's te doorgronden (niet om ze te draaien!) 
is een beetje theorie nodig. Het poolcoördinatenstelsel wordt in de 
wiskundelessen op school niet of nauwelijks behandeld. Eigenlijk is 
dit jammer, want hiermee (en met de parametervorm) kunnen juist 
de mooiste functies (en daarmee de fraaiste grafieken) beschreven 
en getekend worden. 


In een pooleoördinatenstelsel wordt elk punt in het platte vlak met 
twee coördinaten, te weten r en p, bepaald. r is de afstand tussen 
het punt en de oorsprong en p (phi, spreek uit: fie) is de hoek 
tussen de horizontale lijn door de oorsprong en de lijn door de oor- 
sprong en het punt (zie tekening op p.31). 
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—- nul-richting 


oorsprong 


De oorsprong van het poolcoördinatenstelsel leggen we in het 
middelpunt van het beeldscherm (u=128, v=88). Zoals we weten ligt 
de oorsprong van het beeldschermcoördinatenstelsel in de linker- 
benedenhoek van het beeldscherm (HOME-positie), De nul-richting 
in ons poolcoördinatenstelsel is horizontaal en wijst naar rechts. 
p=90° (1/2 radialen) is verticaal naar boven; p=180° (n radialen) 
is horizontaal naar links en p=270° (3m radialen) verticaal naar 
beneden. 


We gebruiken de volgende twee transformatieformules om het, uit 
de functievergelijking berekende, punt P(r‚p) om te zetten in het 
beeldschermpunt P(x,‚y): 


x 
y 


INT(u + r*COS(p) + 0.5) èn 
INT(v + rxSIN (p) + 0.5) 


In onderstaande figuur zien we hiervoor de verklaring. 











175 i 
PE: 
5 + a Pay) = 
E P(rcosp, rsinp) 
rsinpt 
t 
4 
EE 
(us) t reosw 
1 Je 
1 
' 
' 
1 
ii 
0 + El 
0 128 255 


De hoek wp, in de programma's voorgesteld door de variabele p, 
doorloopt steeds in positieve richting (tegen de klok in) het inter 
val van 0 tot 2m. Dit is een interval in radialen (0-360 in graden). 


In alle programma's is als lusvariabele de hoek w in graden gekozen, 
Met de formule p = w.(m/180) (in de programma's p=w*rd) wordt de 
hoek in radialen berekend. Het voordeel van een lusvariabele 
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die het interval 0-360° doorloopt is dat de stapgrootte waarmee de 
lusvariabele steeds wordt opgehoogd een geheel getal kan zijn (1 

bijvoorbeeld) en dat hierdoor de programma's beter leesbaar zijn. 
Een neveneffect is dat alle programma's bijna woordelijk in Pascal 
zijn over te zetten; Pascal kent immers geen gebroken stapgrootte 
in een lus. 


Meer theorie hebben we niet nodig. De volgende programma's spre- 
ken grotendeels voor zichzelf. 











Deze figuur ontstaat als we in programma 12 in regel 1000 
LET r=COS(14*p) nemen. 
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Programma 12 tekent de grafiek van de functie 
r = keos(nw) 


In het programma kiezen we k = 87 enn = 4. 


108 REM programma 12 Bor iek 
van de Functie r=CÔS (4&phi) 

118 CLS 

128 LET u=128: LET „288 

139 LET h=ê,S: EET ka87 


148 LET rd=PI/18 

158 FOR we@ TÔ 58e STEE To 

169 LET pzwtrd: GO SU 

179 LET XzINT woeker sense WE an 
188 LET g=INT (vekartSIN (Pp) +h 


) 
188 IF ER THEN LET Xlex: LET 
: GO TO 239 
LET x@sx: LET Wasy 
@ PLOT x1,yl: DRAW X2-xl,ya- 
F3 wek xi=X2: LET ylsy2 
Le dre Inker THEN 60 TÔ 248 
2 
2 


LET r2CO5 (dtp): LET rzABS 
2 RETURN 


© euwwmv” u 
PSTN + 


ent 


Het programma tekent een achtdelige draaisymmetrische figuur; een 
bloem met acht blaadjes. Als u met dít programma gaat experimente- 
ren zult u snel ontdekken dat voor even waarden van n een 2n- 
bladerige en voor oneven waarden van n een n-bladerige bloem ont- 
staat. Kunt u een tweekleurige bloem maken? 
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Programma 14 tekent sinuskrommen die cirkelvormig gekromd zijn. 








Pp 
168 LET pewérd: GO SUB 1990 
178 LET XzINT (u+r&COS (p) th) 
+ 
Li 


198 LET YeINT (ver SIN (p) 
208 IF pz@ THEN LET Xkex: 
Lay: 89 TO 248 


hb) 
ET 


18 ET Xaz=X: LET ya=y 
220 PLOT x1,g1: DRAU X2-x1,y2 
Zao LET Xisx&: LET yl=y2 
240 NEXT w 
259 NEXT k 
268 IF INKEY&s"" THEN GO TO 260 
276 STOP 
288 
1090 LET raB5O+kaSIN (44p) 
12822 RETURN 


Probeert u eens: 
1000 LET r = 24TAN(24p)+kASIN(2XCOS (SIN (6*PID) 


voor een schitterend spinneweb. De variaties zijn echt onbegrensd! 
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Programma 15 tekent een bloemvormige figuur met blaadjes die in 
het midden aan steeltjes vastzitten. In het programma is n= 4 geko- 
zen. U kunt gerust andere waarden voor n proberen. Het effect is 
hetzelfde als bij programma 12, 


IN 
A 


he REM programma 15 bloemen 
CLS 


129 LET vel28: LET ve88 
130 CET he8.5: LET rdePI/180 
135 LET naé: LET czë. 25: 

REM eerst de snie oemen 


142 FOR ke2® TO STEP 7 
159 FOR w=d TÔ Seo STE 
EL LET peerd: GO SUB 1298 


198 LET W=INT (ver SIN (Pp) +h) 
208 IF en THEN LET Xl=x: LET 


gia: Go TO 248 

18 LET xaax: Ed zy 

gee PLOT Rida Dr W Ke-xl,J2 
Sao EET Klete LET yi=yà 

249 NE 

258 NES 

268 LET c sao: LET Kantise/nsacd 
: REM dan de stele 

220 En Jek TON 

258 LET pajtPpl 

29D LET XL2INT (uer#COS (p) th) 
309 LET YA=INT (ver+SIN (Pp) éh) 
gie LET XBzINT (W+r#C0S (p+PI) 
+ 

ore LET yR2INT (versSIN (pePI) 
+ 

ze PLOT X4,yi: DRAW X2-X1,V2- 
y 

348 NEXT j 

358 IF INKEY$="" THEN GO TO 35@ 
Dee STOP 

aa LET rakeil+ceABS (SIN (nap) 
1823 RETURN 
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Bezit u een kleurenmonitor of een kleurenplotter dan kunt u na 
elke doorloop van de k-lus de afdrukkleur veranderen. Ook de 
stelen van de bloem kunnen een eigen kleur krijgen. Tekent u een 
aantal van dergelijke bloemen naast of onder elkaar dan hebt u een 
begin gemaakt met '‘computer-art'. 
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Spiralen zijn altijd geliefde figuren. Programma 16 tekent logarit- 
mische spiralen of spiralen van Archimedes. Deze laatste soort spi- 
ralen hebben als vergelijking: 

r=e.p 


In het programma op p.40 is voor c de waarde 3 gekozen. Kies 
gerust andere waarden, maar wel tussen 0,5 en 20. 


De logaritmische spiraal heeft (naar Bernoulli) de vergelijking 
zee), 


In het programma dat de onderste spiraal getekend heeft hebben 
we voor k de waarde 85 en voor c de waarde -0,2 gekozen, 


r=2p 


p=g5e 0,20 


Omdat de logaritmische spiraal zich oneindig vaak rond de oor- 
sprong zal winden moet ervoor gezorgd worden dat het programma 
na bepaalde 'tijd' wordt afgebroken. 


EN 
S 
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REM Programma 18 spiralen 


LET ua128: LET ve 
LET Be 5: LET Paeisieo 


„3 
FOR Met TO 19629 STEP 3 
LET pawkrd: GO SUB 1000 
LET XEINT (u+r&COS (pen 
CET VaINT (v+r4SIN (Pp) +h 
Eid Ee THEN LET xlax: LE 
TE 
if x<@ ER en OR y<@ OR 
75 THEN 60 TO 
LET x@ax: LETS 
PLOT X1,g4: bri %e- x1,ve 


REKE: Xlexê: LET yley2 
IF INKEYS&z"" THEN GO TO 268 
STOP 


OORD NE DIE PLC DEEP 
DSD IARGENDE SPS IDUP WES 
SSSS6S8 SSHAUGGISEESDHEE 


LET racap 
RETURN 


en 


Voeg toe regel 245 u = 132 : GO TO 150 en verander in regel 205 
de opdracht GO TO 260 in GO TO 245, en zie hoe dit de figuur ver- 
fraait. 


Voor het tekenen van een logaritmische spiraal maken we c gelijk 
aan -0,2 en herschrijven subroutine 1000 als volgt: 


1000 LET r=854EXP (c*p) 
1010 IF r<5 THEN STOP 
1020 RETURN 
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Als laatste programma met poolcoördinaten geven we het programma 
17 voor het tekenen van een willekeurige, in poolcoördinaten gefor- 
muleerde, continue of niet-continue functie. De programmastructuur 
komt overeen met de structuur van programma 11 (zie p.27). U kunt 
daar de werking van de vlaggen fz en fa nog eens bestuderen. 

Ook de andere variabelen hebben dezelfde betekenis als hun naam- 
genoten in programma 11. Als voorbeeld in het programma hebben 
we een vrij ingewikkelde functie, die niet overal continu is, geko- 
zen: 


De sin(1,5.w) 
1-2.cosp 





De figuur is getekend met a =-2, b =2, Ip =-2, hp =2, wo = 0° 
en wn = 7209, Om de grafiek er optisch mooi uit te laten zien wordt 
er niet getest op een negatieve waarde van r. 


Wilt u een andere functie proberen, verander dan alleen iets in 
de regels 1000 t/m 1090. Verander de regels 1100 t/m 1200 niet. 





Deze zeldzame vlinder heeft als vergelijking: 
___4,sin(1,5p+2) 
Sl ‚cos( Ip) 
cos(p). (1 + DEEL) 


Voor a, b, Ip, hp, wo en wn nemen we in programma 17 de waarden 
-4, +4, -4, +4, O en 720. 


42 Krommen in poolcoördinaten en in parametervorm 


198 REM programma bed grafiek va 
nde functie raf (phiù 

118 CLS 

„180 INPUT "LinKERGRENS VOOR x 
13e INPUT "RECHTERGRENS VOOR X 
‚4e INPUT "ONDERGRENS VOOR Y  " 
‚150 INPUT "BOVENGRENS VOOR Y " 
“ie, deld “STARTWAARDE VOOR PHI 


ze Îneur “EINDWAARDE VOR PHI 


oiRELET, kxea85, o- …à): LET Kyel7 
7 = 
123 EEE hed,5: LET rdaPI/180 


LE 

FOR, wewo TO wn 

LET paytrd: GO SUB zone 
IF fzai THEN LET fan GO 


IF fazi THEN GO TO 302 
LET X2aINT (kX&ix=a) th) 
LET S@eINT (ky&iy=lp) th) 
PLOT x1,yl: DRAW XB-XL,yE- 


LET xiex2: LET yiew@: GO T 


LET XAzINT (KX%0=8) +h) 
LET VEEN (kie Cy LP) th) 
â= 


OS SOSOGUESEGE 


NE uw 

BEEN THEN GO TO 34@ 
LET nal-24C08 (pl: IF ne=@ T 
LET fzel; TURN 

LET raSiN (Sap sm) / 

LET xer#COS (p): (Br yer&SI 
IF x<a OR Xòb RER aise OR yr 


EE far BETON 


SAGTSG ESESSOE, 


PTezZerIe 
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De parametervorm 


De meest interessante krommen krijgen we bij functies waarvan de 
vergelijking in parametervorm wordt opgesteld. Dit houdt in dat 
zowel de x- als de y-coördinaat als functie van dezelfde, derde, 
parameter t worden uitgedrukt. In de natuurkunde zien we vaak 
de tijd als parameter (vandaar de t). In de wiskunde is de parame- 
ter t bijna altijd een hoek in radialen. 


Zo is de parametervorm van een cirkel met straal r en het middel- 
punt in de oorsprong: 


Xx = reost en y = rsint 


Uit deze twee vergelijkingen kan gemakkelijk de cartesische vorm 
x?+y?=r? gedistilleerd worden. De vergelijking van een ellips met 
het middelpunt in de oorsprong en met een halve lange as a en een 
halve korte as b is: 


xXx = acost en y = bsint 


Het is eenvoudig programma's te ontwikkelen voor het tekenen van 
stelsels concentrische of excentrische cirkels. Ook stelsels ellipsen 
zijn, dankzij de parametervorm, eenvoudig te tekenen, We doen dit 
echter niet, omdat de figuren niet zo bijzonder zijn en vrij bekend, 
In hoofdstuk 1 hebben we trouwens al een ellips met behulp van de 
parametervorm geprogrammeerd. 


De volgende programma's tekenen figuren die u mogelijk nog niet 
kent en die zich voor computer-graphies bijzonder goed lenen. 
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Programma 18 tekent een Lissajousfiguur. In het algemeen is de 
parametervorm voor een dergelijke figuur: 


x= Kj-sin(f,.t+p) + kj-cos(f,.t) 
y= kg-sin(fg.t+pg) + k‚-cos(f,.t) 


Hierin is p de fase, f de frequentie en k de amplitude. Om een mooie 
figuur te krijgen hebben we de volgende waarden gekozen: 


kj=k3=87, Ê=16, p‚=0, Êgel7, p3=39 en kgekgefgefg=0 


102 REM programma 1 
Lissajousfiguren 


188 BRINT “TOETS K1,F1,PA,K2,FR 
Ì8o ENEUT kaslljPiske, fa 
INT ÍTOETS Ko,F3,PS,K4,F4 
13, pe,k4r ft 

Tv 


K3 
138: LET'veóe 
8.5: CET rdePI/180 


en 
e 
ee 
D' 
0: 
ej 
Z 


ao HD 
ourrrOrmmz D 
ONMRD AAT 


e zcd 
uu 


368 
d: GO SUB 1002 
+X th) 


(u 
Lv 4dth) 

EN LET xx: LET 
L 


A44 

HEX 
SunuS 
HE 

VAZZe A 


E 
erijso 


TO & 
LET x@asx: LET yeey 

PLOT xA,y1: DRAW X2-Xl,uR- 
LET Xlexa: LET pie 


EXT w 
A end THEN GO TO 29@ 


EC 
GOGO ADD DE DD DIDO et 


ET KekAsSIN (FARt+P1) +KRRC 
65 
EE ANKRASEN CESRtep 3) akar 


P SPOOON AAI PGNPGOPINZO 
-G-SSIOES SMEGOSEGSEE 


ETD GHZ 
se ha 


POROr 
© 

LI 

©. 

Et} 

Me 

5e 

C 

5 

2 





Experimenteer met het programma! Zoek zelf de juiste waarden voor 
de constanten zodat u mooie figuren krijgt. Natuurkundigen zullen 
u hierbij graag helpen; zij weten welke waarden u moet nemen! 
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Programma 19 presenteert een niet alledaagse figuur. Deze figuur 
wordt dikwijls de ‘vliegenkopfiguur' genoemd. Om de vliegenkop 
horizontaal op het beeldscherm te krijgen moet de parameter t het 
interval 90°-450° doorlopen. 


M programma 19 vliegenkop 


u=1i28: LET va87 
h23,5: LET Ke3@: LET rd 


@ 

OR we9D TO 458 STEP 3 

LET tawtrd: GO SUB 1000 
LET Xe INT (U+X+h) 

Li Y=INT (V+y th) 

I BU THEN LET xiex: LET 
b 
Pi 


aad 
2ux: LET Ed 
JT X1,y1: DRAW X2-KL,VR- 


lexe: LET yi=yâ 
KEY&e"" THEN GO TO 23 


€ 
X 
Zex 


S 
SaGEnED PRDEWECHEEEE DELE 


STOP 


„Her Xak4SIN (2#U) #(2,5+COS 
LET gekaerCOs (3et) 
RETURN 


DE*SGR ODP SOPOIDRPHEANEE 
S5-666 SOS GOLUSSELENGEGE 


eran 


Het tekenen van de x- en y-as is achterwege gelaten om de 'kop' 
beter te laten uitkomen. 


Voor een echte vliegenkop met ‘ogen’ moet u in de appendix kijken. 


Hoe geraffineerder u de functies x = f(t) en y = f(t) kiest, des te 
ongewoner worden de figuren. Probeer hier hoe creatief u kunt zijn. 


De 'huiswiskundigen!' bij computerfabrikanten hebben vaak hun 
eigen lievelingsfuncties. Die zie je dan ook vaak in een demonstra- 
tieprogramma optreden. 
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Programma 20 tekent een stelsel krommen. Hewlett-Packard publi- 
ceerde deze tekening enige tijd geleden. Onder het programma staat 
een tabel met waarden voor a en b, die heel mooie figuren opleve- 
ren. Probeer zelf andere combinaties. U zult verrukt zijn over wat 
u ziet! 


a=-6/b=1 





programma 20 vlinders 
Usl2S: LET va87 

„0.5: LET rdaePI/160 

k LET kya8 

UT "Toets a @n bin “;a 





AG HEEPOD 
mor zmmmt mn 


na=3 TO 

R wed TÔ See STEP 1 

Er" tewerd 

Den dk (U en DRC 
AT VARENS IN (Den tbes 


ET XXeINT (CU4KXEX +) 
GEBR CASE ve 
bi exx: 

GO TO 2e® 
LET X@axx: LET yaw 
PLOT X1,yi: Baat x ENT 
LET Xiex@: LET yiey@ 
EXT u hels 


N 
EXT _n 

F INKEYEs"'" THEN GO TO 300 
TOP 


fe) 


5 

++ 

Dororrone Diaz 
r 


d H 
BEOVDLEND DUNZANEERET eee 


u 
€ 
e 


PSOONEBHCREN P SODND UPGDPG 
G8EEE GOPEGS-E-SSESE LOGSGE 


GHZ 


Tabel voor mooie figuren 








b 1 2 1 2 1 1,5 edion 



































Kies in deze programma's de stapgrootte voor hoek w (FOR w=.....) 
gerust 1; dit duurt wat langer maar geeft mooiere tekeningen. 
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Het laatste programma uit dit hoofdstuk tekent ‘supersymmetrische! 
figuren. Ook nu geven we een tabel met waarden voor de parame- 
ters a, ben e. Het is ongelooflijk hoe de figuur totaal verandert 
als we andere waarden voor één of meer parameters kiezen. Het 
idee voor programma 21 is afkomstig uit het Amerikaanse tijd- 
schrift Creative Computing. 


a=2/b=1/e=? 





a=6/b=6/e=4 


102 REM programma 21 symmetrisc 
he, krommen 

110 CLS 

„i#O INPUT “Toets a,b,c in “sa, 
13e LET v=12s: LET vase 

149 LET he@.5: LET rdePI/100: L 
Er kaê? 

is CLS 

180 FOR wad TO 268 

de LET tuwerd: LET rakeSIN (Cc 
tan LET X@aINT (U+r#COS (att+ 
„19® LET VAsINT WerasiN (bete 
200 IF wed THEN LET xiex@: LET 
EE ee 60 TO 

„Ee PLOT x1,yi: DRAU xR-XL,VR- 
È22 LET Xlexa: LET vieva 

230 NEXT u Ee 

240 IF ENKEvan"" THEN GO TO ao 
aso sro 


Kijk in de appendix hoe u de waarden voor a, b en c samen met 
de figuur op het scherm kunt afdrukken. 
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Tabel voor mooie figuren 




















a # 6 4 1 3 Dei ets En 
b 7 6 | 6 1 3 2 EAA 
e 3 4 1 4 5 dk aan pede 

| 


























Het zou jammer zijn als u bij het zien van een mooie figuur de waar- 
den voor a, b en c niet hebt genoteerd. Verander het programma 
zo, dat na afloop van het tekenen de waarden van a, b en e op het 
beeldscherm of op de printer worden afgedrukt (zie de appendix), 
Ook heel leuk zijn de combinaties: 


a 20 -40 100 
b zl -40 -0,5 
e 3 10 3 


a=4/b=6/e=1 





























a=-1/b=-40/e=-0,5 

















4 Tekenen van driedimensionale 
figuren 


In dit en het volgende hoofdstuk gaan we ons bezighouden met het 
tekenen van driedimensionale liehamen (zo heet dat in de wiskunde) 
zoals kubussen, prisma's, pyramiden, kegels en bollen, en met het 
tekenen van driedimensionale grafieken van functies. Zo'n drie- 
dimensionale grafiek is een vlak in de ruimte met een vergelijking 
van de vorm z = f(x,y). Zo krijgen we de „Beken nde 'Mexicaanse 
hoed! (zie p.73) als we de functie z = e” 0) tekenen (e is het 
grondtal van de natuurlijke logaritme; e = 2,71828...). 


Er bestaan veel programma's waarmee driedimensionale figuren 
getekend kunnen worden. Waarom we toch ook in dit boek derge- 
lijke programma's laten zien komt, omdat er aan de meeste van deze 
programma's drie nadelen kleven: 


1. _De meeste programma's voor driedimensionaal-tekenen zijn 
voor een bepaald graphies-systeem ontworpen en zijn slechts 
met veel inspanning geschikt te maken voor andere systemen. 


2. De wiskundige theorie die aan het driedimensionaal-tekenen 
ten grondslag ligt wordt zelden of zeer summier behandeld. 


35, Veel programma's zijn uiterst ingenieus ontworpen en draaien 
op de kleine mierocomputers, in BASIC, heel langzaam. 


Met de volgende programma's en stukjes theorie hopen we deze drie 
nadelen uit de weg te ruimen. 


Er bestaan verschillende methoden om een driedimensionaal lichaam, 
bijvoorbeeld een kubus, in een plat vlak te projecteren. Wij hante- 
ren de projectiemethode, die u wellieht op school gebruikt hebt (of 
nog gebruikt). Stelt u zich een doorzichtige kubus voor met ribben 
van draadijzer. U staat voor deze kubus een beetje rechts van het 
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midden en kijkt er zo'n beetje schuin bovenop. Uw gezichtsstralen 
projecteren elke hoek, met de daaraan verbonden ribben, van de 
kubus in een, achter de kubus liggend, projectievlak. Zo ontstaat 
de bekende 'schuine' afbeelding van een kubus. 


Blame 











De verticale en horizontale ribben worden op ware grootte getekend. 
De ribben die naar achteren lopen worden verkort weergegeven. 


a!'=k.a k = verkleiningsfactor 


De rechte hoek tussen de ribben BA en AD in het grondvlak lijkt 
eveneens kleiner te zijn geworden (zie a in de bovenstaande 
figuur). U kunt gemakkelijk nagaan dat door het vastleggen van 
alpha (bijvoorbeeld 45°) en k (bijvoorbeeld 0,5) de projectierich- 
ting eenduidig bepaald is. 


Wat we nodig hebben zijn transformatievergelijkingen die voor een 
bepaalde a en k de coördinaten x,y,z van een punt op de ruimte- 
lijke kubus projecteren op de coördinaten x' en y' van het gepro- 
jecteerde punt in het platte (projectie-)vlak. Deze formules vormen 
de kern van alle volgende programma's. 
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Afleiding van de transformatievergelijkingen 


Het projectievlak is ons beeldscherm. We gebruiken ook nu een 
schermresolutie van 256 bij 176 punten. De oorsprong van het ruim- 
telijke coördinatenstelsel (xyz) moet precies in het midden van het 
beeldscherm geprojecteerd worden. De oorsprong van het beeld- 
schermeoördinatenstelsel ligt in de linkerbenedenhoek. De x'-as 
wijst naar rechts; de y'-as naar boven. 


De afbeelding van het ruimtelijke punt P(x,y,‚z) op het punt 
P'(x',y') in het platte vlak komt tot stand met behulp van de vol- 
gende transformatievergelijkingen : 


x'=u+x+tk,y.cosa èn 
y'=v+(k.y.sing + z) 


Hoe we aan deze vergelijkingen komen is uit onderstaande figuur 
direet (nou ja, direct, .) duidelijk. 














Nemen we c = k.cos(a), s = k.sin(a) en h = 0,5 dan zijn de trans- 
formatieformules in BASIC : 


INT (u*x+exy+h) voor de x-coördinaat en 
INT (v+s*y+z+h) voor de y-coördinaat. 


Meer wiskunde hebben we voor dit hoofdstuk niet nodig. 
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Programma 22 laat zien hoe je van elk, door rechte lijnen begrensd 
lichaam (polyeder) een projectie kan maken. Hiervoor is het nood- 
zakelijk dat de waarden van de coördinaten x,y,z in alle hoekpun- 
ten bekend zijn. Als voorbeeld van het gebruik van dit programma 
tekenen we een kubus met een ingeschreven achtvlak (oktaëder). 
De gevolgde programmeertechniek leidt ook bij een viervlak, prisma 
of pyramide tot het gewenste resultaat. 


We nemen aan dat het middelpunt van de kubus samenvalt met de 
oorsprong van het ruimtelijke coördinatensysteem. Om het program- 
ma ‘sneller’ te maken nemen we de coördinaten van de acht hoekpun- 
ten ABCDEFGH van de kubus en de coördinaten van de zes hock- 
punten IJKLMN van de achthoek op in DATAregels. Het tekenen 
van de ribben wordt door de letterstring z$ bestuurd. Zo betekent 
z$ = "ABBCCDDA" dat de computer van A naar B, van B naar C, 
van C naar D en van D naar A rechte lijnen moet trekken. Met de 
CODE-opdracht halen we steeds één letter uit de string z$ en maken 
er een getal van (A=1l, B=2, C=3, ...….). Hiermee moet de werking 
van het programma duidelijk zijn. Kies voor alpha 45° en voor k 
0,5; dat geeft de mooiste projectie. 
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198 REM programma 2ë 
KOESS het B-vran 


118 CLS 
189 INPUT "Alpha in graden (B) 
zee ANPUF vwerkteiningsfactor ú 


Pea LET umi28: LET v 


«8: 
154 LET h=0.5: LET vatBi/ieo 
168 LET want LET cak&COS (w) 
LET sak SIN 


E tw 
120 ga xild): bam yiia): DIM z 


2 


1 


2 
Ze IereGOEELLEGGEX HIN! 


z 


8 
jel TO 14 


Sn 

ae 
DO! 

0 


NEX 

ges 

FOR_n=1 2 : 
BERB'Ee U? LET velen (zet 
FOR mel TO \-1 STEP 2 
LET i=COD ED 


Ed 

E Zn 
LET jxCODE z&in) (m+1) Bd 
sytij: DRAW XU) 


BARD GOD 


A BORRE 
GGSSESHESESSS-SORE GSESSSG 


kelk 


n 
NKEY Sa" THEN 60 TO 920 


9 


-60,-50,-69,60,-62,-50 
62, 62,-bo,-68, do, - 
=66,-62,66,50, 5e 6 
3550 EI4SE, 


EE EEA 

8807026070 ordr0:5d 
ARSCEDORAEBFOGDHEFPGG 
“EU IK IL IMJKKLLMMJJNKNL 


VOUGID WHZ + 
PDDDDD Sumze 
444AAÁ 
EDDDDP Vn 


ZORDS OON S EN 
v U 
Lb D 
ERE 
» P 
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In programma 23 wordt in een kubus een zadelvlak getekend. Het 
geeft het idee van een gekromd vlak in een kubus. 





Bekijk de bovenstaande figuur. De diagonalen BG en ED zijn elk 

in acht (n) gelijke stukken verdeeld en de zo ontstane punten zijn 
paarsgewijs door een rechte lijn met elkaar verbonden. Kies in het 
programma voor n, bijvoorbeeld, de waarde 32. De bovenste kubus 
hebben wij getekend, de onderste de computer! 
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199 REM programma 23 
an EN RODE REL Pddeuviar 


188 Zeur “Alpha in graden (45) 
130 NPUT vverkleiningsfactor ( 


‚s) 

135 INPUT “Hoeveel Lijnen (92) 

2 

140 LET usi2g: LET ve 

188 LET VaaeS: LET Vleer 

169 CET vaatrd: „LET EER EES 0) 
: LET sak#SIN (w 
gij oi td: BIM veer: omm zet 
189 FOR je1 TOS 

199 READ K‚p:z 

Beo LET XLI BINT tuex+cs th) 
210 LET yIJBINT (V4sey+z+h) 
228 NEXT j 

252 CLS 

250 RERD ze: LET LaLEN (26) 

250 FOR mei TO L-1 STEP à 

270 LET isxCODE zein) -64 

252 LET JsCODE Zi (m+1)=54 
288 PLOT xlil,ytid: DRAY xj) 
KE ld 9 

300 NEXT m 

505 FOR Jee TO_n 

318 LET XAMINT (x(@) ee iX (7) =X 
(8) Zn th) 

315 LET VIMINT (v(@) eeu (7) =y 
B) zn th) 

329 LET XBZINT (XB) tje (4) =X 
(5) ) /n+h) 

325’ LET vasINT (y(B)4Jtty td) -y 
(5) zn th 
„39e PLOT xi,vs: DRAU x8-x1,v2- 
335 NEXT J 
340 PLOT Xt@),via): DRAW x (7) -x 
Dy CP) -y (2) 

46 PLOT X(5),ytB): ORR Xx (4) Xx 
B yi) =y (5) 

344 IF INKEYän"" THEN GO TO 347 
343 Srok 

58 DATA 5045043045, SL LED 
560 DATA 69,bo ‚6e, - 
578 DATA <36 EER EON ae se 
359 DATA ‚éo,6b,-da,éo ó 
420 DATA “aBscoDDAREBFCAONEFFGO 


I 


55 
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Met programma 24 kunnen we, naar keuze, cilinders, kegels of 
afgeknotte kegels tekenen. Als we voor rl en r2 dezelfde waarde 
invoeren, krijgen we een cilinder, Als r2 kleiner is dan rl krijgen 
we een afgeknotte kegel en voor r2 = 0 krijgen we een kegel. 


We plaatsen het lichaam zo dat het grondvlak bij z = -60 ligt en het 
bovenvlak (de top) bij z = +60, Het programma tekent zeven 
breedtecirkels met een onderlinge afstand van z = 20 en 16 lijnen 
op de kegel- (of eilinder-) mantel loodrecht op de breedtecirkels. 
Tot slot wordt verticaal de z-as getekend. 
198 REM programma 24 
cylinders, kegels 
118 CLS 
“Alpha in graden (45) 
NPUT "Verkleiningsfactor ( 
NPUT "straten r1 en ra 
Usi28: LET va88 
10,8: CÉT rdaPi/180 
drsirler@) 76 


weasrd: LET 
úw): 


zó, 
Ho 
zi 
u: 
4 
ä 


© 
RO FTE Hi 
Lelielied 
44 
> 
u 
a 


ror mm 


vermelde Pi 


0: 


OPNAM 


n=ë 
LET s=k4SI 


EOGSOErS 


z 


SCONBHAP ES 
De 4d 

a 

E 

Le 

* 

br] 

kel 

a 


za-69 TO 69 STEP 28 
Er ruf lenede 
FOR wad TO 360 STEP 39 
LET wizwerd 
LET x=r#COS (wi): LET yer 


wi 

IF wi<>® THEN GO TO 308 

LET XAsINT (U+X+C#Y+N) 

LET YIaINT (V+Styezeh) 
TÔ 348 


LET XRRINT (Uaxecuyth) 
LET YazINT (V+S&yeZ+h) 
PLOT x1,yi: DRAW XB-x1 


LET xiexà: LET yleya 
XT 


“ 
LET nan+1 
NEXT Z 
FOR wa@ TO 360 STEP 23 
LET wiewerd 
LET xar14COS (wi): LET ver 


E 
POR ROMIG TORD 


‚ya 


BRERPDERRERGIGOGEDE GOOR 
HOSGSEEEE GEGEOLEEZEECEEE 

Z 

hl 


EFO IGAR OENE 


1 N twi) 

Ö LET XLZINT (uex+CHy+h) 

12 LET Y1aINT (vesty-6Oth) 

23 LET XarR#COS (Wi): LET var 
BRSIN (wid 

3) LET XRMINT (U4X+CHY +) 

43 LET YRAINT (vasey+6O+h) 

58 PLOT XA,y1: DRAW XR-X1,VB- 
y 

58 DET Mi @ RAL 178 

U@: 

258 TRINE Et ie ENEN TE To «ee 
438 STOP 


Vindt u de breedtecirkels te ‘grof! neem dan in regel 230 als stap- 
grootte van de FOR-lus bijvoorbeeld 3. Het tekenen duurt dan wel 
een stuk langer! 








Tekenen van driedimensionale figuren 57 








Erg leuk is het tekenen van een bol met breedtelijnen en meridianen. 
Zoals bekend is kunnen we elk punt op het boloppervlak eenduidig 
vastleggen met de straal van de bol (r), de geografische breedte 
(@) en de geografische lengte () (zie figuur). Het omzetten van 
deze boleoördinaten r‚p,‚A in cartesische coördinaten (x,y,z) 
geschiedt met de volgende drie vergelijkingen: 


Xx = PCOSPCOSÀ 
y = rcospsin) 
z = rsinp 





De bolvergelijking in cartesische vorm is overigens x?+y?szt=r?. 
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In programma 25 wordt de hoek wp door w en de hoek } door p 
vertegenwoordigd. Kies voor alpha 90° en voor k 0,5. Andere 
waarden vervormen de bol, 


programma 25 bol 

PUT "Alpha in graden (99) 
TÊNPUT “verkleiningsfactor U 
@ INPUT "Straal, maximaal 78 


ET ual2S: LET va 
ET huô.5: LET EaSPirie0 


tcd 
cuk&C0S úw): LET SmkesI 


goe ee 


AROCDHESOONOA :e-G DPO 
a 





z 


BRPPRCHEG- XGEGG-GEORGEERENDDOKDEDLDNENDE ve 


GRLEDPPOG +OAIBEGLOPOOFONRS + 0 


We-93 TO 99 STEP 15 
ET wiswaerd: LET rier#COS 


FOR P=3 TO 368 STEP 3 
EEE Plupard: LET xel.15er 


pi) 
LET garikSIN (pi): LET ze 
(wad 


IF pzd THEN LET XIsINT (U 
LET YL&INT (v+Z+h): GO TO 


LET X@RINT (U4X+CHy th) 
LET va@zINT ie: Ee hrleed- HOO 
PLOT x1l,yi: U Ne-xl,ve 


LET Xiex@: LET viav@ 
XT p 


a 
ze 


+ 


AZ 
arOnzZ 
15 


EEROZHerE Fe 


2 B 182 STEP 15 


u 
en 

Poi td 

Gao 16 360 STEP 3 

T winwerd: LET rier#C0S 


T Xei,15er1#CO8 (P1) 
T WarieSiIN (Pi): LET ze 


) 

wed THEN LET x1sINT (u 
dant YIMINT (v+styerzeh 
ET XRINT jarber ete) 

T WAEINT (V4Sty+zeh) 
LOT X1,V4: DRAW XB-X1,a 


ta 
z 

€ 

nk 

Uem MD 


G6s65 ssogrsoorsssreRs SOPETSHEHE-GEGESSSS, 
E3 


N 

NEXT p 

är INKEYän"" THEN GO TO 458 
TOP 








Zie de appendix voor ‘mooiere! bollen. 
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Een geliefd demonstratieprogramma van computerleveranciers is een 
om zijn as draaiend driedimensionaal lichaam. Meestal wordt als 
‘draai-as' de z-as gekozen, die dan samenvalt met de lengte-as van 
het lichaam. 


Mierocomputers die alleen een BASIC-vertolker bezitten zijn te lang- 
zaam om zo'n draaïing real-time uit te voeren. Er zijn namelijk nogal 
ingewikkelde trigonometrische berekeningen voor nodig, die tame- 
lijk veel tijd in beslag nemen. Daarnaast duurt het tekenen van het 
liehaam in een bepaalde stand in BASIC zo'n 1 à 2 seeonden. Zouden 
we in een BASIC-programma steeds het lichaam tekenen, uitwissen, 
draaien, tekenen, uitwissen, draaien, ..… enz., dan zien we het 
draaiende lichaam in plaats van een vloeiende beweging een schok- 
kerige beweging maken. De beste oplossing voor dit probleem is de 
volgende. 


Deel de totale middelpuntshoek van 360° in n even grote hoeken. 
Voor elk van deze hoeken 
200 


u 


= je = 360° 
1 » @ Djs eeens & n.w. 360 


n 1 

berekenen we van te voren de coördinaten van de hoekpunten van 
het lichaam. Voor elke stand slaan we deze hoekpuntscoördinaten 
in een array op. Al deze berekeningen voeren we in BASIC uit bij 
een leeg beeldscherm. Als de berekeningen klaar zijn wordt een 
machinetaalprogramma uitgevoerd dat de eerste groep coördinaten 
uit de array haalt, het lichaam tekent, na een bepaalde tijd het 
beeldscherm wist, de volgende groep coördinaten ophaalt, 
enz. Deze oplossing heeft een redelijk vloeiende draaiing tot gevolg. 
Dit programma is geschreven op een CBM 3016 die een 6502-micro- 
processor bezit, Hier doen we echter anders (zie volgende pagina). 
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Programma 26 laat een driezijdig prisma om de z-as draaien. De z- 
as is zowel de lengte-as als de symmetric-as van het prisma. Het 
programma tekent een prisma en wacht vervolgens met het draaien 
en opnieuw tekenen tot we een toets indrukken. Elk volgend prisma 
wordt over zijn voorgangers heen getekend. We kunnen het draaien 
dus vertraagd gadeslaan. 


EN 
e 
Ss 


REM Eeten 2e 
draaiend prisma 


“Alpha in graden (46) 
INPUT “Verkl&iningsfactor « 
INROT "oraasnoek in graden 


De 
0-6 :G6 





LET ade LET vas 
EET he@,S: LET E8sBr/ieo: L 


LET eed LET cek#COS (u) 
K&SIN 


NK: 


m 
GEER RRODARDKDDEE PAP Er 7 


* 


XC HY th) 
Stg -6+h) 


Vestyr60r 


= ++ 


T 
DPOO ONGEHPONPSOONN OAARAG 


GSE SCEEZEESELEEEME- GE 


jl: DRAW %ijt 
DRAU X(UJ+ 


J+4): DRAW 
Wij+d) 


ze € 
Dd 


4 
TI 
m 
z 
e 
bed 
4 
[e] 
e 
@ 
e 


388 STOP 
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Programma 27 tekent de projectie van een regelmatig twintigvlak 
(ikosaëder). Het is bekend dat er slechts vijf 'echt-regelmatige'! 
veelvlakken (polyeders) zijn, namelijk cen tetraëder (regelmatig 
drievlak), een kubus (regelmatig zesvlak), een octaëder (regelmatig 
achtvlak), een dodekaëder (regelmatig twaalfvlak) en een ikosaëder 
(regelmatig twintigvlak). Een ikosaëder heeft 12 hoeken, 30 ribben 
en 20 vlakken. Een vlak is een gelijkzijdige driehoek. Als we pro- 
gramma 22 willen gebruiken om zo'n ikosaëder te tekenen dan heb- 
ben we de coördinaten van alle 12 hoekpunten nodig. Als we het 
ikosaëder in een speciale stand zetten kunnen we de hoekpuntcoör- 
dinaten met behulp van de techniek van de ‘gulden snede! relatief 
eenvoudig berekenen. Wie zich voor de hierachterliggende wiskun- 
dige theorie interesseert moet er maar eens cen meetkundeboek op 
naslaan. De 'gulden-snede'-techniek deelt een lijnstuk AB in twee 
stukken AT en TB op zo'n manier dat de volgende evenredigheid 
geldt: 


AB : AT =AT : TB 
| monmnnmanscmnsnesnteenn 
A T B 


Kiezen we voor AB de lengte 1, dan kunnen we de lengte van AT 
berekenen; immers 


AB _ AT 


AT TB 





de 5 
dus Tt alst = AT 


dus dan moet gelden t? = 1-t ofwel t?+t-1=0, 
Met de alombekende abc-formule berekenen we nu 





EESTE 
102 2 
De positieve wortel geeft: t-d 5 {5 


Dit is de lengte van het grootste van de twee stukken die we krij- 
gen als we een lijnstuk, ter lengte 1, volgens de ‘gulden snede! in 
twee stukken verdelen (t = 0,618: 1-t = 0,382). 

Deze waarde wordt in regel 180 berekend en in de regels 220-330 
gebruikt om de hoekpuntcoördinaten van het ikosaêder te bereke 
nen. Om te zorgen dat het twintigvlak voldoende groot getekend 
wordt zijn alle coördinaten met 70 (f) vermenigvuldigd. Kies voor 
alpha (a) 90° en voor k de waarde 0,4. Andere waarden verteke- 
nen het beeld. Alpha = 180° en k=0,6 geeft trouwens iets onver- 
wachts! 


62 
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192 REN programma 27 ikosaeder 
118 CLS 
122 INPUT "Alpha in graden (90) 
jä 
ee INPUT "verkleiningsfactor « 
43 LET uz128: LET va88 
152 LET h 5: LET rduPI/180 
152 LET waatrd 
nee eer Cak&COS (w): LET sek&SI 
u 
vee LET ta{S9R (5)-1)/2: LET fz 
189 CLS 
208 DIM x (12 
218 OIM zilg 
283 LET Xx (1) L 
ET zildmef 
239 LET x(2)a®: LET yl@)s-ftt: 
LET za) mf 
243 LET x (31 LET y(9)a@: LET 
Ees afer 
58 LET X(4)m-f: LET yié)n®: LE 
T zie) sefet 
260 LET x(B)efat: LET wief: L 
ET_zi5)a0 
270 LET x(B)u-fet: LET yibef: 
LET z(6) 22 
280 LET x(7)afat: LET y 07) fb: 
LET z(7) u@ 
290 LET X(S)z-ftt: LET y(O mf 
LET z{8)2@ 
300 LET xi9)laf: LET v(9)e2: LET 
zie ret: (20) f: LET y (18) 8 
x mt: v „0: 
LET z(i®) ef st 
329 LET X(11)2®: LET yiiefet: 
LET z(iiaf 
330 LET x(l2)sd: LET yilgie-f et 
: LET Zil) ef 
348 FOR nal TO 3 
En READ Z&(n): LET LuLEN (z&í 
n 
362 FOR tul TO l-1 STEP 2 
370 LET imCODE Zin) (m) -84 
388 LET JsCODE z&in) (m+1) -64 
zer LET XIRINT (uexti)scay tin 
Kd LET yIsINT (veseyli)sz tij 
bd LET XBRINT (Uexij) teng tjd 
Nn 
+a0 LET VAsINT (Vestyij)+z (jd 
+ 
430 PLOT x1,va: DRAW xe-x1,vg 
die NEXT _m 
<52 NEXT n 
4568 IF INKEYGe"" THEN GO TO 460 
473 STOP 
EE 
Saa 
512 
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In het vorige hoofdstuk hebben we ons uitvoerig beziggehouden met 
het tekenen van driedimensionale lichamen. We herhalen hiervan nog 
eens de belangrijkste punten: 


Te We gebruiken de parallelprojectie om een driedimensionaal 
liehaam met n hoekpunten P(x,y,z) in een plat vlak te tekenen. 
Voor de projectiehoek alpha (a) gebruiken we bij voorkeur de 
waarde 45° of 60°. De schuin-naar-achteren lopende ribben 
worden korter getekend dan ze in werkelijkheid zijn. Als ver- 
kleiningsfactor kiezen we vaak 1/2 of 1/3. 


2. Elk punt P(x,‚y,z) van het ruimtelijke lichaam wordt op een 
punt P'(x',y') in het projectievlak geprojecteerd. De hierbij 
gebruikte projectieformules zijn: 

x'=u+txXtk.y.cosa 

y'=v *+k.y.sina + z 


Zoals gebruikelijk zijn u en v de coördinaten van het midden 
van het beeldscherm, Tot nu toe hebben we u = 128 en v = 88 
verondersteld. 


In de volgende BASIC-programma's zien we de bovenstaande trans- 
formatievergelijkingen in de vorm: 


xg = INT (utxxtexyy+h) 
yg = INT (v+tsxyytzth) 


De betekenis van de gebruikte variabelen is: 


x8,y8 : beeldschermcoördinaten x',y' 
XX,YY de lopende coördinaten x,y 
e k*COS (ward) 

s : k*SIN (werd) 

w : projectiehoek in graden 
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rd 1 factor 7/180 voor omrekenen van graden in radialen 
k : verkleiningsfactor 


Tekenen van driedimensionale functies 


Als paradepaardje van menige demonstratie van Hoge-Resolutie- 
Graphies wordt vaak een programma gebruikt dat een of andere 
fraaie grafiek van een driedimensionale functie tekent. Als leek 
vraag je je af hoe ze weten welke functies ze moeten nemen, hoe het 
tekenen van de grafiek van zo'n functie geprogrammeerd wordt, en 
hoe ze het programmatechnisch voor elkaar krijgen dat de ‘onzicht- 
bare lijnen’ ook inderdaad niet getekend worden. Bij dergelijke 
demonstraties wordt doorgaans geen programmalisting getoond. 
Mocht dit wel het geval zijn dan is het programma vaak zo machine- 
afhankelijk dat het omzetten van het programma naar een andere 
machine lastig, zo niet ondoenlijk is. 


In dit hoofdstuk hopen we u antwoord te kunnen geven op de vol- 
gende vragen: 


1. Hoe vind ik 'mooie! driedimensionale functies? 


2 Hoe ziet een algemeen programma voor het tekenen van een 
dergelijke functie eruit? 


3, Hoe onderdruk je het tekenen van de onzichtbare lijnen 
(hidden lines)? 


We beginnen met de definitie van een driedimensionale functie. 


Elke functie van de vorm z = f(x,y) heet een driedimensionale func- 
tie en vormt een, meestal, gekromd vlak in een driedimensionaal 
(ruimtelijk) coördinatenstelsel. Speciaal voor niet-wiskundigen volgt 
nu een voorbeeld van een functie z = f(x,y) en het gekromde vlak 
dat als 'grafiek' bij de functie hoort. 


Stel dat we voor z = f(x,y) de volgende vergelijking nemen: 


(xt +y? 
ze 7195) (e is het grondtal van de 
natuurlijke logaritme; 


e = 2,718284183...) 


Voor elk punt (x,y) in het (platte) X-Y-vlak kunnen we nu de 
daarbij horende waarde van z uitrekenen. Als x = 0,5 en y = 0,1 
dan is 
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z g0,25:0,0D 
z 0,26 
= 8 


> 


z = 0,771, want 2,71828... tot de macht -0,26 is 0,771 


Zet nu (in gedachten) in het voetpunt P(0,5;0,1) in het X-Y-vlak 
een staafje met een lengte van 0,771 rechtop. Doe ditzelfde voor 
nog een groot aantal punten (x,y). 





Leg nu over al deze staafjes een elastische folie. Dit (golvende) 
stuk folie vormt een goed model van het vlak met vergelijking 
OE Ct di] 
EE . 
Hoe dit eruit zou zien kunt u zien op p.73. 
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Zoek met behulp van programma 7 uit hoofdstuk 2 een willekeurige 
eontine functie die symmetrisch ten opzichte van de y-as is. De 
grafiek van de functie moet ‘bergen en dalen' vertonen (de functie 
moet maxima en minima hebben). Erg geschikt zijn trigonometrische 
functies (sin, cos) en combinaties van deze functies. Ook geschikt 
zijn functies waarin alleen termen voorkomen met teven-machten! 
van x en exponentiële functies. 


Enkele voorbeelden: 





: y= sin(x) 


y= cos(x) - sos) á sene) A goat) 
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Bekijk nu de grafiek van zo'n functie in het interval -a*xza. Dit 
‘stuk grafiek gaan we draaien rond de y-as. Zo ontstaat een, ten 
opzichte van de y-as, draaisymmetrisch driedimensionaal vlak. Als 
we nu de y-as als z-as kiezen, dan wordt de vergelijking van een 
dergelijk draaisymmetrisch ruimtelijk vlak (rond de z-as) van de 
vorm 


ze f(r) met r = Vx?+y? 


We zullen dit aan een eenvoudig voorbeeld proberen te verklaren. 
We kiezen als voorbeeld de eenvoudige paraboolvergelijking y=x?. 


var? 





Stel we bekijken de grafiek voor -2sx7 2: 











Als we nu de x-as rond de y-as laten draaien (de x-as komt als het 
ware loodrecht het papier uit), dan ontstaat een draaisymmetrisch 
vlak rond de y-as. Dit is een kegel: 
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/ 





Elk punt P op de kegelmantel heeft de eigenschap y=r?, waarbij r 
de afstand is van het geprojecteerde punt P' tot het punt 0. We 
kunnen het vlak V dat door de ronddraaiende x-as wordt gevormd 
als volgt tekenen: 











Welnu, als we nu de z-as als y-as nemen en we nemen de y-as in 
vlak V loodrecht op de x-as, dan zien we dat de afstand r geschre 
ven kan worden als 


r=Vx?y? (stelling van Pythagoras, zie rechts in de 
volgende tekening) 
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en nu wordt y=r? geschreven als z=r? met r? = x?+y?, dus de ver- 
gelijking van de paraboloïde (kegel) wordt dan 


z= xt? 











Dus y=x? wordt bij draaiing om de y-as en bij een z-as die de plaats 
van de y-as inneemt z = x?+y?. 
Voor de bovengenoemde functies geldt het volgende: 








bd er wordt zee O7) 

y= NE wordt ze sm) 
vcd cos(3x) , IE) 5 eene) wordt 
y =cos(r) - gaste) ES cosi) En sos{le) 


met r = Vx?+y? 


Voor deze klasse van draaisymmetrische figuren ontwikkelen we een 
algemeen tekenprogramma. 
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Antwoord op de tweede vraag 


Bekijk de onderstaande tekening. U kijkt recht bovenop het ruimte- 
lijke vlak van een driedimensionale functie. We willen dat het op het 
grondvlak geprojecteerde vlak een vierkant is met zijden van 2a. 
Straks zullen we zien dat, als we een mooie tekening van zo'n ruim- 
telijk vlak op het beeldscherm willen maken, a=93 een heel goede 
keus is. 





(0,0) 





a 


ann 
Dy 

Dx 
je en y=-93 


a a 











Als ‘doorgewinterde! programmeur ziet u natuurlijk direct dat het 
programma voor het tekenen van het vlak een geneste FOR-NEXT 
(dx,‚dy) zal bevatten! Als we een vlak als grafiek van zef(x,y) wil- 
len tekenen, moeten we namelijk voor een aantal combinaties van x en 
y de functiewaarde z = f(x,y) berekenen, We beginnen met y =-93 
en laten x met stapjes dx (bijvoorbeeld dx=3) van -93 maar +93 
lopen. Voor elke combinatie x,y met y = -93 berekenen we de func- 
tiewaarde z = f(x,y). Zo ontstaan de punten (x,y,z) van het ruim- 
telijke vlak waarvan de geprojecteerde punten (x,y) in de boven- 
staande tekening op de onderste zijde van het vierkant als bolletjes 
getekend zijn. Nu verhogen we y met dy en laten x weer van -93 
tot +93 lopen. Dit geeft de tweede reep van het vlak. Op deze wijze 
ontstaat het hele vlak, waarvan de punten (x,y,z) natuurlijk nog 
getransformeerd moeten worden naar beeldschermcoördinaten (x',y'). 


Een eerste globale opzet voor het tekenen van een 'vierkant'-stuk 
vlak is: 
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FOR y=-93 TO 93 STEP dy 
FOR x==93 TO 93 STEP dx 
GO SUB 1000 


NEXT x 
NEXT y 


In de subroutine 1000 wordt voor een punt (x,y) de waarde van 
z = f(x,y) berekend. Met de bekende transformaties wordt vervol 
gens het punt P(x,y,z) overgebracht naar een punt P'(x',y') op 
het beeldscherm. Nu kan het berekende punt P(x,y‚z), dat in het 
vlak ligt, als het punt P'(x',y') op het scherm worden getekend. 


Jammer genoeg lenen niet alle driedimensionale functies zich voor 
de intervallen -93sx593 en -93sy<93; vandaar dat de waarde voor 
a met een INPUT-opdracht ingelezen wordt. Als we in het program- 
ma toch -93sxxs93 en -93Syys93 kiezen, moet de ingetoetste 
waarde a als volgt gebruikt worden: 


x= En = a 
= XX.95 en Y=YY-55 
dat wil zeggen er geldt: -asxsa en -asysa. 
De z-waarden van met name de trigonometrische driedimensionale 
functies zal tussen -1 en 1 liggen. Om deze waarden wat ‘op te bla- 
zen! kan een vermenigvuldigingsfactor (k1) worden ingetoetst. 
Goede waarden voor kl liggen tussen 30 en 80. 


We kunnen ons programma nu als volgt opschrijven: 





begin tekenprogramma voor de functie z = f(x,y) 





maak beeldscherm schoon 





lees w‚ k, a en kl in 





geef u, v,‚, h, rd, c‚ enz. een waarde 





maak beeldscherm schoon 





voor y is -93 tot +93 in stapjes dy doe 
voor x = -93 tot +93 in stapjes dx doe 








bereken de functiewaarde z = f(x,y) 





bereken de beeldschermcoördinaten xg en yg 








teken het punt (xg,‚yg) op het beeldscherm 





wacht tot een toets wordt ingedrukt 





einde tekenprogramma voor de functie z = f(x,y) 
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Zoals beloofd laten we nu zien waarom de intervallen -93sxxs 93 en 
-93syys93 zo geschikt zijn om driedimensionale functies te tekenen 
met een hoog oplossend vermogen. We willen graag de fraaie pro- 
jectie met q=45° en k=0,5 gebruiken. Wat worden dan de beeld- 
schermcoördinaten xg en yg voor de hoekpunten 'linksonder' en 
‘rechtsboven! van het vierkant op p.69? Deze punten bepalen 
immers of de hele grafiek op ons beeldscherm van 256 bij 176 punt- 
jes getekend kan worden. Voor kl nemen we 50, dat ligt zo'n beetje 
tussen 30 en 80! Nemen we voor de z-waarde ook 50, dan gaat het 
punt (x,y,z) met coördinaten (-93,-93,50) over in het beeldscherm- 
punt (xg,yg) met 


xg = INT(128-93-0,5.93.cos(45)+0,5) > 2 
yg = INT(88-0,5.93.sin(45)+504+0,5) > 105 


Deze waarden liggen inderdaad binnen ons 'graphic-scherm' 
(Osxgs255 en Osygs176) en de ‘breedte! (256) van het scherm 
wordt heel goed benut, kijk maar naar de coördinaten (xg,yg) voor 
het punt (93,93,50) rechtsboven: 


xg = INT(128+93+0,5.93,cos(45)+0,5) > 254 
yg = INT(88+0,5.93.sin(45)+50+0,5) => 171 


We benutten dus haast de hele schermbreedte (2£xgs254) voor het 
tekenen van de driedimensionale figuur. In programma 28 (regel 
270) en in programma 29 (regel 310) is rekening gehouden met even- 
tuele negatieve xg-waarden en xg-waarden groter dan 255. Als dit 
voorkomt, worden deze waarden op respectievelijk 0 en 255 gezet, 
zodat het programma niet door een foutmelding afbreekt, 


Op de onder- en bovenrand (yy=-93 en yy=+93) gelden vaak heel 

kleine z-waarden, zodat bij a=45°, k=0,5 en z=0 de grafiek op het 

scherm zal liggen tussen 
n z-waarde 

yg = INT(88+0,5.93.sin(45)+0+0,5) > 121 

en yg =INT(88-0,5.93.sin(45)+0+0,5) > 55 


Ons tekenvlak is dus ongeveer 2sxgs254 èn 
55sygsl2l 


In het bovenstaande programmavoorstel wordt gebruik gemaakt van 
‘puntgraphies'. Hierbij worden de punten puntje voor puntje gete- 
kend. Om een enigszins acceptabele tekening te krijgen moeten erg 
veel ‘puntjes! berekend worden (dx=1,dy=1). Dit duurt in BASIC 

(geneste lussen en vele trigonometrische berekeningen) erg lang. 

Voor de onderstaande tekeningen betekent dit al snel een tekentijd 
van meer dan 30 minuten, soms wel een uur. Dit is weinig bevredi- 
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gend! Veel sneller gaat het als we de 'vectorgraphies'-methode 
gebruiken. Deze techniek hebben we in alle voorgaande programma's 
gebruikt; we verbinden steeds twee naast elkaar liggende punten 
door een recht lijntje. Hierbij hoeven we lang niet zoveel punten 

te berekenen als bij de 'puntgraphics'-techniek. Bovendien ziet 

de tekening er beter uit. 


Als antwoord op de tweede vraag komen we dan met het volgende 
algemene programma voor het tekenen van draaisymmetrische ruim- 
telijke vlakken. In programma 28, en ook in programma 29, zien 
we voor xXx de variabele q en voor yy de variabele p*. 


108 REN genoaranee 2e grafiek va 





‘Alpha in graden (45 
Verkleiningsfactor £ 
ur “REChterorens voor x 
‚“Vaergrotingstactor (3 
ueikê: LET veen 

haô.5: LET rduPI/280 
Sake09s werd): LET sak 
É 

dxn3: LET dya6: LET afa 








« @ …— 
@ 
JP HeRAhez Zi Zee 
neef 2 ek la 
„Jed d 


ren TO 93 STEP dy 


TO 93 STEP dx 
f: GO SUB 1980 
TE (u+q+cep+h) 
Ml 


acof 
Á 
c 
Hi 


HHHHHTT TORD 


(V+SRP+Z+h) 
N LET X 


75 Ti et jans 
3 xiuxg: 
y 


pe 
« 

rc 
F4 
led 


T Xt: er ve 
PLOT x1,yl: DRAW at. 


LET Xlex@: LET viewe 
En 
IF INKEYGa"" THEN GO TO 95% 
STOP 

LET ZakIkEXP (axe ey ey) 
RETURN 


Met dit programma maakt u de tekening op p.73. Kies op uw ZX 
Spectrum de volgende waarden: w=45°, k=0,5, a=3 en k1=70. U 
krijgt dan het mooie 45°-perspectief. 





*In een FOR-NEXT-lus mogen alleen éénlettergrepige lusvariabelen 
voorkomen, vandaar q en p in plaats van respectievelijk xx en yy. 
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= bte? 
zee 9D 




















U zult over de bovenstaande tekening niet tevreden zijn. De voor- 
grond is goed, maar in het achterste stuk zijn ook alle lijnen die 
voor ons onzichtbaar zijn getrokken. In onderstaande tekening zien 
we het resultaat met dezelfde waarden voor w‚ k, a en kl, maar 
waarin de niet-zichtbare lijnen ook niet getekend zijn. Hoe krijgen 
we dit voor elkaar? 
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Bekijk de onderstaande tekening. Een punt P(x,y,‚z) van het te 
tekenen vlak is alleen dan zichtbaar wanneer: 


a. het punt 'hoger' ligt dan alle voorgaande punten met dezelfde 
x-coördinaat of 

b. het punt 'lager' ligt dan alle voorgaande punten met dezelfde 
x-coördinaat. 


Vanzelfsprekend geldt dit ook voor de op het beeldscherm gepro- 
jecteerde punten P(xg,yg). 
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Plxryrz) 
« 




















We gaan als volgt te werk. We gebruiken twee arrays (hl(xg) en 
h2(xg)) waarin we voor alle xg-coördinaten in ons vlak de respec- 
tievelijk kleinste en grootste yg-waarde bewaren. Vinden we voor 
een punt met een bepaalde xg-waarde een bijbehorende yg-waarde 
die óf kleiner is dan hl(xg) óf groter is dan h2(xg), dan is het 
punt zichtbaar en wordt hl(xg) of h2(xg) gelijk gemaakt aan deze 
nieuwe kleinste of grootste yg-waarde. Geldt voor een bepaalde xg 
dat hl(xg)<yg<h2(xg), dan is het punt (xg,yg) niet zichtbaar. 
hi(xg) en h2(xg) zijn te zien als twee horizonnen. 


Als voor een bepaald punt P1(x,y,‚z) uit het vlak het beeldpunt 
Pl(xg,yg) berekend is, gaan we kijken of yg kleiner dan of gelijk 
aan hl(xg) is (yg<=hl(xg)). Is dit zo dan is het punt Pl op het 
beeldscherm zichtbaar. Nu zetten we de vlag fl op 1 en we maken 
hixg) gelijk aan yg (hl(xg)=yg). Is yg groter dan hi(xg) 
(yg>hl(xg)) dan is Pl onzichtbaar en blijft de vlag fl op 0 staan. 


Nu gaan we naar het volgende punt Pz(x,y,z) uit het vlak (we ver- 
hogen x met dx). Opnieuw berekenen we het beeldpunt P2(xg,yg). 
Weer wordt gekeken of yg<=hl(xg) en zonodig wordt de vlag f2 op 

1 gezet. De punten Pl en P2 worden alleen dan door een rechte lijn 
(lijntje) verbonden als beide vlaggen fl en f2 de waarde 1 hebben, 

dus als fixf2=1. 


Ditzelfde doen we voor de tweede ‘horizon’ h2(xg). We kijken of 
yg>=h2(xg), ………. ‚ enzovoorts. 


Omdat we bij een vaste y-waarde de x-waarde steeds met dx opho- 
gen (bijvoorbeeld dx=3) en hiermee steeds twee buurpunten 
Pl(xg,‚yg) en P2(xg,yg) berekenen slaan we als het ware een aan- 
tal punten over waarvoor geen hl(xg)- en h2(xg)-waarden berekend 
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worden. Door lineaire interpolatie tussen de horizonwaarden (hl(xg) 
en h2(xg)) van twee buurpunten zouden we de horizonwaarden 
voor de hiertussenliggende punten kunnen berekenen. We hebben 
deze waarden wel nodig, omdat het kan voorkomen dat, als we y 

met dy verhogen, we xg-waarden vinden waarvoor nog geen 
hl(xg)- en h@(xg)-waarden berekend zijn. 


Een programma dat de hierboven geschetste werkwijze zou volgen 
zou, in BASIC geschreven, veel te langzaam zijn. Om de tekensnel- 
heid acceptabel te houden, dat wil zeggen niet langer dan 10 minu- 
ten tekenen voor één figuur, passen we de volgende wereenvoudi- 
gingen toe: 


1. We berekenen alleen de bovenste horizon, die we h(xg) noemen. 


2, We passen geen lineaire interpolatie toe bij het berekenen van 
de array-waarden in h(xg). Alle beeldpunten in een interval 
ter lengte dx krijgen dezelfde h(xg)-waarde. 


Als we de stapgrootte dx maar klein genoeg kiezen (een goede waar- 
de is dx=3), krijgen we ondanks deze twee simplificaties toch accep- 
tabele tekeningen. Hieronder volgt een korte uitleg van het pro- 
gramma. 


De regels 


200 DIM h(256) 

210 FOR L=1 TO 256 
220 LET h(l)=-1000 
230 NEXT Ll 


zorgen ervoor dat, voor het tekenen begint, het hele scherm zicht 
baar is, Dit is zo als de onderrand van het scherm de horizon is. 
De waarde -1000 geeft in feite een horizon die ver onder het beeld- 
scherm ligt. De array h bevat in principe voor elke pixel op de 
horizontale beeldschermlijnen (x-waarden) een horizonwaarde. Als 
we dx=1 nemen, hebben we inderdaad 256 van deze horizonwaarden 
nodig. In programma 29 gebruiken we dx=3, hetgeen betekent dat 
we op de horizontale beeldlijnen steeds 2 pixels overslaan, waar 
voor dus ook eigenlijk geen horizonwaarden bekend hoeven te zijn. 
In feite hebben steeds drie naast elkaar liggende x-pixels dezelfde 
horizonwaarde. We laten dit aan de hand van een voorbeeld zien. 


Stel dat in programma 29 voor een bepaald punt P(x,y,z) berekend 
is dat dit punt geprojecteerd wordt op het punt met beeldscherm 
coördinaten xg=40 en yg=126 en laten we aannemen dat dit berekend 
is voor een waarde xx>-93 (in het programma voor q>-93). Na het 
berekenen van xg en yg (regels 290 en 300) komt het programma 
bij regel 330. Omdat we even aannemen dat xx(q)>-93 gaat het pro- 
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gramma verder met regel 370 (door GO TO 370 in regel 330). De 
regels 370 en 380 luiden: 


370 LET f2=0: LET L=INT(xg/dx)+1 
380 IF yg>h(l) THEN LET f2=1: LET h(lD=yg 


De vlag f2 wordt eerst op nul gezet. Vervolgens wordt de index 1 
voor de horizonarray bepaald. Voor xg=40 en dx=3 vinden we 1 = 
INT(40/3)+1). Dit geeft 1=14, In feite zouden xg=39 en xg=41 ook 
l=14 opleveren. Steeds leveren drie punten op de x-as dus dezelfde 
horizonwaarde op. Dit komt omdat dx=3 gekozen is. Nemen we dx=1 
dan levert elke xg-waarde tussen 0 en 255 een andere l-waarde op. 
Voor xg=0 en dx=l krijgen we I=INT(0/1)+1). Dit geeft l=1; het 
eerste horizonarray-element. Voor xg=255 en dx=1 vinden we l= 
INT(255/1)+1, dus 1=256. Dit is het grootste horizonarray-element. 
Voor dx=l gebruiken we dus alle horizonarray-elementen. Terug 
naar ons voorbeeld. De waarde voor l ís dus 14 als xg=40, yg=126 
en dx=3. In regel 380 wordt vervolgens bekenen of de berekende 
yg-waarde groter is dan de op dat moment geldende horizonwaarde 
bij xg=40. Is yg inderdaad groter dan did horizonwaarde (yg>h(l) 
is waar), dan is het punt (xg‚yg) zichtbaar en wordt de nieuwe 
horizonwaarde bij xg=40 gelijk aan de waarde yg. Stel dat de hori- 
zonwaarde h(14) gelijk was aan 130. Regel 380 test dan: 


IF 1262130 THEN LET f2=1: h(14)=126 


maar 126 is kleiner dan 130, dus blijft 130 de horizonwaarde bij 
xg=40 en de vlag f2 wordt niet op 1 gezet maar blijft 0. Dit betekent 
dat het punt (xg=40, yg=126) een onzichtbaar punt is en dus wordt 
er geen lijn getrokken uit het laatst getekende punt naar dit zojuist 
berekende punt. In regel 410 wordt hierdoor fl ook weer nul (fl= 
f2), zodat we eerst twee keer een f2-waarde van 1 moeten hebben 
om weer een zichtbaar lijntje te kunnen trekken. Op deze manier 
ontstaat een tekening met verborgen lijnen (hidden lines), waarvan 
een aantal na het programma is afgedrukt. In programma 29 zien 

we de variabelen q en p de rol van xx en yy vervullen. 


— 2 2 
Met w=45°, k=0,5, a=3 en k1=70 krijgt u de grafiek van z=e C*°*Y5) 
zoals die een paar pagina's eerder is getekend; de niet-zichtbare 
punten die achter de hoed liggen zijn inderdaad niet te zien! 


Hoe kleiner we dx en dy maken, des te gedetailleerder zal de figuur 
worden. Maar ook 'des te langer zal het tekenen duren'. De combi- 
natie dx=3 en dy=5 is een goed compromis tussen de tekensnelheid 
en de mate van gedetailleerdheid. 
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100 REM programma 29 
grafiek van zefix,y) 

& hidden Lines 

BUT “Alpha in graden (45 

eur “yerkieiningsfactor ‘ 

PUT “Rèchtergrens voor x 

"a 

ur „rverarotingsfactor Kt-] 

u=iëë: LET ve68 

0.5: LET rdaP1/180 

=k4COS (ward): LET ser 


"3: LET dy=5: LET afm 


gi 


ZNZ ZI 


S-61:90S8 
7 
G 


e 
@ 


ROZ HO Mara He H-HO 


OFMFOM 


* 
OHNE HOP Wee 
me MEI - 2 

Ae HH 


* 
PDIOIUN EA ee 


Ó: 


AMD 
3 


R_pu-93 TO 93 STEP dy 

LET gaptaf 

FOR q=2-93 TO 93 STEP dx 

LET xeataf: GO SUB 1008 

LET X9eINT (U+q+cepth) 

CET YOZINT (v+sep+Zth) 

F X9<@ THEN LET x9=@ 

FR X9>258 THEN LET, had 
Pr yo OR +175 THEN PR 

‘royte ki ': STUP 

IF q>-83 THEN GO TO 378 

LET faz@:; LET LeINT (X9/d 


Er A THEN LET flul: 
ey 
LET Xiexg: LET yiayg: GO 
LET f@=3: LET LaINT (x9/d 
IF ygrhit) THEN LET Fagel: 
hCUmyg 
LET x2eXg; LET vert 
@ IF fi+fâai THEN PL xL,y 
DRAW XE-XL,YR-YL 
LET Xiexë: LET yieya: LET 
- 
2 NEXT a 
2 NEXT Ë 
2 Pd INKEY gu" THEN GO TO 440 
1089 LET zekisEXR C-Xax-y ty) 
RETURN 


HH 


APP POO NINA PG 


E 
GAZEEELELKDDD 


x 


h 


mn 


br PUTG-GOUF GT 
eN NEES 
EN 
& 


x 
OIOPGEGIGFEG SAOSEEELLESPEELEZEEG 
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Dit programma is een algemeen programma voor het tekenen van 
draaisymmetrische ruimtelijke vlakken. Als u andere figuren wilt 
maken, verander dan alleen de functie in de subroutine 1000 en 
kies mogelijk andere waarden voor de variabelen. 
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Nu volgt een aantal voorbeelden. 


1000 LET r=SQROaxtyky derd 
1010 LET z=k1#COS(r)-COS(Bkr)/34COS (5ér)/5-COSC7kr)/7 
1100 RETURN 


U kiest: w=45, k=0,5, a=180°, k1=35. 





cos(3r) , cos(5r) _ cos(7r) 


z = cos(r) — Lj 5 7 


1000 LET r=SQROkxtyay derd 

1010 IF r=0 THEN LET z=k1: RETURN 
1020 LET z=k1#SINCrd/r 

1100 RETURN 





U kiest: w=45, k=0,5, a=1080, k1=50. 
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1000 LET r=SQRxixtyty) 
1010 LET z=k1kEXP(-COS(r/16)) 
1100 RETURN 


U kiest: w=45, k=0,5, 2=90, k1=30. 





z= exp(-cos(j) 


Voor de volgende twee tekeningen geldt: 


1000 LET r=SAR(xtxtyiy) 
1010 LET z=k1*COS(r/16) respect. LET zek1SIN(r/16) 
1100 RETURN 





U kiest: w=45, k=0,5, a=90, k1=50. 


Ox 


” Tr nn 
z= cos(75) z= sin( 





r 
Te) 


Met dit programma hebben we talrijke andere grafieken gemaakt. 
Elke keer als u een ‘leuke functie! tegenkomt, kunt u meteen de 
daarbij horende draaisymmetrische figuur maken. 
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Tot slot volgt nog een programma voor het tekenen van een fraaie 
viertoppige grafiek. U bent hem misschien wel eens in een compu- 
tertijdschrift tegengekomen. Om de snelheid op te voeren is de sub- 
routine 1000 in het hoofdprogramma opgenomen. U kunt dit natuur- 
lijk ook in alle andere programma's doen. 


108 REM R52ICSAMA 30 
oie functie 


148 CLS 
„520 INPUT “alpha in graden (45- 


3 
138 INPUT “verkleiningsfactor £ 
„3,55 tk 

142 LET Uwi2S: LET va 

158 LET hs@,5; LET KasPi/ieo 
152 LET cek#COS (werd): LET sek 


SIN (ward) 

zee LET dxa3: LET dyz5: LET Kia 
188 DIM h (285) 

198 FOR Lul TO 256 

208 LET hile-1290 

218 NEXT L 

228 CLS 

239 FOR y=-83 TO 93 STEP du 

ad LET mietos (WrRePI/9I-PIJ + 
250 FOR x2-93 TO 93 STEP ax 
288 LET meaC08 (xx2rPI/93-PI) 

+ 

20 LET zeklemiema 

280 LET X9eINT (Uax+ctyth) 
298 LET YAzINT svn) 
308 IF _X>=83 TH Go 348 

„ee LET f1iz8: tr Cant ard 
+ 

328 IF yg>hil) THEN LET flai: 
LET hil)myg 

TO ae LET Xlaxg: LET yimyg: GO 

Sed LET fF@x=@: LET LeINT ixg/d 
+ 

352 IF yg>htl) THEN LET Faal: 
LET hieyg 

360 LET Kasxo: LET v2u 

378 IF flxfasl THEN Bret Xl 

1: DRAW XR-X1,yR-yd 

38e LET XlaX@: LET yizy2: LET 
leta 

390 NEXT Xx 

408 NEXT E 

418 IF INKEYSa"" THEN GO TO 410 
«28 STOP 








ze (eos Eme) 
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6 Turtle-graphics en 
LOGO-simulatie 


LOGO is vooral door de turtle-graphics beroemd geworden. Er is 
geen andere programmeertaal waarmee zo eenvoudig zeer moeilijke 
grafische figuren gemaakt kunnen worden dan met LOGO. In dit 
hoofdstuk zullen we vijf LOGO-programma's in BASIC vertalen. Deze 
BASIC-programma's bevatten wel 25 tot 30 regels, terwijl LOGO 
hiervoor slechts 3 of 4 regels nodig heeft. 


“Eerst LOGO, daarna andere programmeertalen" 


is het motto van vele informatici en pedagogen die zich met de 
invoering van de informatica op scholen bezighouden. 


Wanneer men deze stelling wil begrijpen, moet men niet alleen de 
taal LOGO en haar mogelijkheden, maar ook de omgeving waarin 
LOGO ontwikkeld werd goed kennen. Een korte historische terug- 
blik lijkt daarom op zijn plaats. 


Het schrijven van computerprogramma's is een intellectuele presta- 
tie en een creatief proces. Aan de wijze waarop wij programma's 
ontwikkelen, herkennen wij direct onze manier van denken. Aan 
het Massachusetts Institute of Technology (MIT) heeft men al in het 
begin van de zestiger jaren een speciale programmeertaal ontwik- 
keld, met behulp waarvan men kunstmatige intelligentie ging bestu- 
deren. Men noemde deze taal LISP, een afkorting van LISt- 
Processing. 


In LISP werden programma's geschreven die een met kunstmatige 
zintuigen uitgeruste elektromechanische muis in staat stelden een 
uitweg uit ieder willekeurig doolhof te vinden. In LISP worden pro- 
gramma's ontwikkeld die de computer tot een medespeler met leer- 
vermogen maakt. Hoe meer spelletjes de computer tegen een mense- 
lijke tegenspeler speelt, des te beter wordt hij. Goede zetten van 
de tegenstander onthoudt hij en legt hij vast in zijn geheugen, ter- 
wijl hij de eigen slechte zetten voor toekomstige spelletjes uit zijn 
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geheugen wegveegt. Een ‘afvalprodukt' van deze studies aan het 
MIT zijn de moderne schaakprogramma's. 


LOGO is een hoog ontwikkeld dialect van de LISP-taal. Seymour 
Papert, leerling van de bekende onderzoeker Jean Piaget en mede- 
werker aan het MIT, heeft in twaalf jaar tijd LOGO ontwikkeld. Jean 
Piaget heeft in zijn bekende boek Hoe kinderen leren de kinderlijke 
denkstructuren laten zien. Hieruit blijkt dat tekenen en knutselen 
tot de eerste creatieve handelingen van kinderen behoren. Wij zul- 
len zien hoe LOGO deze bezigheden ondersteunt. 


Amerikanen gaan door voor een volk dat graag experimenteert. Het 
was hun echter duidelijk, dat jonge kinderen, wij denken dan aan 
zes- tot dertienjarigen, niet in staat zijn een programmeertaal 
zoals BASIC, laat staan Pascal, te leren en algorithmen voor nume- 
rieke, niet-numerieke of grafische problemen te schrijven. Daarom 
heeft men een ‘kinderlijke! taal gemaakt (waarachter echter een 
imponerende software schuilgaat), waarmee het kind met gemak 
tekeningen kan maken. 


Wordt het LOGO-systeem ingeschakeld, tegenwoordig meestal een 
64K-mierocomputer met bijbehorende software, dan verschijnt in 
het midden van het beeldscherm een kleine driekhoek, waarvan de 
top naar boven wijst. De kinderen noemen die driehoek 'turtle', het 
Engelse woord voor schildpad. Met behulp van eenvoudige bevelen 
kan het kind de schildpad willekeurig over het beeldscherm laten 
rondlopen. En al near gelang het wenselijk is laat de turtle wel of 
geen spoor na. 


De volgende LOGO-opdrachten zijn beschikbaar: 


FORWARD 100 100 passen vooruit 


BACK 50 = 50 passen achteruit 

RIGHT 90 = draaiing van 90° met de klok mee 
LEFT 45 = draaiing van 45° tegen de klok in 
PENUP = haal pen van papier 

PENDOWN = zet pen op papier 

HIDETURTLE = haal Turtle van het beeldscherm 


Normaal gesproken geldt de toestand "PENDOWN", Als we bijvoor- 
beeld de hoofdletter F op het scherm willen tekenen, kan dat in 
LOGO als volgt: 


FORWARD 100 RIGHT 90 FORWARD 50 
RIGHT 90 PENUP FORWARD 50 PENDOWN 
RIGHT 90 FORWARD 50 

HIDETURTLE 


Voor 100, 50, 90 en 45 kunnen ook andere waarden gekozen worden. 
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Omdat LOGO een vertolkend systeem is, wordt elke opdracht (na 
het geven van RETURN) direct uitgevoerd. Zo op het oog lijkt het 
slechts leuk speelgoed! Laten we nu eens een LOGO-programma 
bekijken; liever spreken we van een procedure: 


TO VIERKANT 
REPEAT 4 (FORWARD 75 RIGHT 90) 
END 


Het LOGO-systeem weet dat tussen TO en END een procedure (een 
stuk programma) gedefinieerd wordt. De procedure tussen TO en 
END (in ons voorbeeld VIERKANT genoemd) maken we zelf, Zouden 
we bovenstaande procedure intoetsen, dan zal op het beeldscherm 
een vierkant met zijden ter lengte 75 getekend worden. De program- 
meurs onder u kennen vast de opdracht REPEAT waarmee een her- 
halingsstructuur geprogrammeerd kan worden. 


Als we de procedure VIERKANT in het LOGO-systeem ingevoerd 
hebben, kunnen we hier ook gebruik van maken. LOGO onthoudt 
alle procedures die wij zelf invoeren. We kunnen dan ook later 
nieuwe procedures ontwikkelen waarin we gebruik maken van eer- 
der ingevoerde procedures (zoals VIERKANT). Nu volgt een proce- 
dure met een parameter: 


TO VIERKANT :ZIJDE 
REPEAT 4 (FORWARD :ZIJDE RIGHT 90) 
END 


Toetsen we nu VIERKANT 150 <RETURN?> in, dan zal LOGO cen 
vierkant met zijden van 150 tekenen. Niets verhindert u een proce- 
dure in te bedden in een volgende procedure en die weer in een 
volgende en die weer ………. ‚ enzovoorts. Alleen de beschikbare 
geheugenruimte is hierbij een remmende factor. Bekijk het volgende 
LOGO-programma. 


LOGO-programma nr. 1 


TO VIERKANTPATROON 
REPEAT 8 (FORWARD 20 LEFT 45 VIERKANT 75) 
END 


Dit eenvoudige drieregelige LOGO-programma tekent een patroon 
van acht vierkanten (zie p.87). We zullen dit programma straks in 
BASIC vertalen. Het zal blijken dat hiervoor enige wiskundige en 
programmeertechnische kennis nodig is. Het LOGO-programma kan 
door een leerling van de lagere school gemaakt worden, terwijl het 
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BASIC-programma dat deze acht vierkanten tekent slechts door 
scholieren van de hoogste klassen van het voortgezet onderwijs kan 
worden gemaakt. Het wordt nog moeilijker als we intikken: 


TO STROOK:AANTAL 
REPEAT:AANTAL (FORWARD 100 VIERKANTPATROON) 
END 


Tikken we vervolgens in STROOK 5 <RETURN> dan maken we 

heel eenvoudig een fraai ornament. Probeer dit maar eens in BASIC 
of Pascal. LOGO wordt pas echt interessant als we van de mogelijk- 
heid gebruik maken dat een procedure zichzelf aanroept (recursivi- 
teit). De turtle-graphics berusten op het principe dat we een pro- 
cedure parameters geven en dat deze procedure zichzelf steeds met 
andere parameterwaarden aanroept. Bekijk de volgende procedure: 


TO VEELHOEK:ZIJDE:HOEK 

FORWARD:ZIJDE LEFT :HOEK 
VEELHOEK:ZIJDE:HOEK — hier roept de procedure 
END zichzelf aan 


Als we nu VEELHOEK 200 144 intikken, zal het LOGO-systeem een 
regelmatige vijfpuntige ster tekenen. We zullen echter op de STOP- 
toets moeten drukken om het tekenen te stoppen. Brengen we in 
deze procedure een stopmechanisme aan en laten we de procedure 
steeds de waarde van ZIJDE veranderen, dan krijgen we de wel- 
haast bekendste LOGO-procedure: 


LOGO-programma nr.2 


TO TURTLE:ZIJDE: GROEI: HOEK 
FORWARD:ZIJDE LEFT :HOEK 
MAKE:ZIJDE:ZIJDE+:GROEI 
IF:ZIJDE>200 THEN STOP 
TURTLE:ZIJDE:GROEI: HOEK 
END 





Als er een dubbele-punt voor een LOGO-woord staat, weet het 
LOGO-ssyteem dat het om een naam van een variabele gaat en niet 
om de naam van een procedure of een LOGO-opdracht. In de boven- 
staande LOGO-procedure zien we hoe de procedure zichzelf steeds 
met een nieuwe waarde voor ZIJDE aanroept. Ook dit programma 
vertalen we in BASIC. Hiermee kunt u elke denkbare rechtlijnige 
turtle-grafiek maken. U hoeft alleen de waarden van de invoerpara- 
meters (ZIJDE, GROEI en HOEK) te veranderen. 
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Het derde LOGO-programma tekent een reeks vierkanten in spiraal 
vorm. Het zijn bekende figuren in het land van de computer- 
graphics’. 


LOGO-programma nr. 3 


TO VIERKANTSPIRAAL:ZIJDE:HOEK 
IF:ZIJDE>150 THEN STOP 
VIERKANT:ZIJDE LEFT:HOEK 
VIERKANTSPIRAAL:ZIJDE+4: HOEK 
END 


Ook dit LOGO-programma vertalen we in BASIC. Met deze drie pro- 
gramma's hebt u niet alleen met LOGO maar ook met de ‘turtle- 
graphies' kennis gemaakt. Wilt u meer weten van deze graphies- 
wereld, dan wijzen wij op het boek Turtle Geometry: The Computer 
as a Medium for Exploring Mathematics van Harold Abelson en 
Andrea dì Sessa, uitgegeven door Cambridge, MA:MIT Press, 1981, 
Liefhebbers van wiskunde, informatica en LOGO vinden in dit boek 
een schat aan informatie. 


Natuurlijk zijn we in LOGO niet gebonden aan ‘rechte lijnen', In 
LOGO hebben we de beschikking over alle wiskundige functies en 
bewerkingen. Alle programma's uit dit boek kunnen in LOGO 
geschreven worden. Zij zouden dan vast korter, overzichtelijker 
en begrijpelijker geweest zijn. In LOGO kunnen we heel eenvoudig 
met lijsten en tabellen werken. In LOGO hoeven we de variabelen 
niet te declareren. Achter een naam in LOGO kan gewoon een getal, 
een string, een n-dimensionaal veld en nog veel meer schuil gaan. 
Ook kan het de naam van een procedure zijn. We hoeven ons niet 
te bekommeren om het reserveren van geheugenruimte (DIM- 
opdracht in BASIC), noch om het aangeven van het type (real, 
integer, character, boolean, enz.); het LOGO-systeem regelt dit 
zelf. 


Wie in LOGO programmeert leert automatisch netjes programmeren. 
Begrippen als procedure, recursie, herhaling, toekenning worden 
op een natuurlijke manier aangeleerd. Met een paar opdrachten is 
een kind al in staat ingewikkelde figuren te tekenen. Met LOGO 
degraderen we de computer niet tot een programmeerbare zakreken- 
machine, maar halen we er alles uit wat erin zit. 
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Vertaling LOGO-programma nr.1 in BASIC 


De kern van elk LOGO-graphics-programma wordt gevormd door 
twee opdrachten: 


FORWARD (resp. BACK) :ZIJDE 
LEFT (resp. RIGHT) :HOEK 


Hoe vertalen we deze opdrachten in BASIC? Bekijk hiertoe de vol- 
gende figuur. 


y 


dx=s. cos (w) 
äy=s.sin(w) 








HOME x 


De PEN staat in punt P(x1,y1) en moet in de richting van hoek w 
over een afstand s verplaatst worden. In Ps(x2,y2) aangekomen 
moet de PEN (eigenlijk de turtle) een hoek van dw° linksom maken. 
De volgende BASIC-opdrachten voeren dit uit: 


INPUT x1,y1,w,s,dw 

LET h=0.5: LET rd=PI/180: LET wi=wird 
LET x2=Cx1ts*COS (w1) thd 

LET y2=Cy1tS*SIN (W1)thd 

PLOT x1,y1: DRAW x2-x1 „y2-y1 

LET x1=x2: LET y1=y2 

LET wewtdw: IF w>360 THEN LET wew-360 
LET wi=wird 





Deze opdrachten komt u in alle volgende BASIC-programma's tegen. 
Meer theorie is niet nodig! De programma's moeten met bovenstaan- 
de informatie gelezen kunnen worden. 











Vertaling LOGO-programma nr. 1 in BASIC 


Experimenteer met programma 31 LOGO-1 vierkantpatroon. Met 
w=90°, s51=20, dw=45, 52=75 en n=8 Krijgen we de onderstaande 


tekening. 

















1 
nen (LOGO-4) 
ten startpun 


hting (grade 
sing "is 


am 
kantsp 
F 


ee 55 
a sou 


140 INPUT wareber ‚Linksom 
“jd 

Mel INPUT “zijde vierkant “is 
162 INPUT “aantal vierkanten " 





a 
179 LET ha3,5: LET rdaPI/188: L 
EI Ag 


188 FOR Jel 

Ps a LET Kink iss 1kCOB (wi) 
Er LET yRaINT (y1+sieSIN (wi) 
+ 


228 PLOT X1l,yl: DRAW XB-xl,VE- 


y 

23) LET xlexa: LET vient 

243 LET axsmx1i: LET ay 

LET grt: IF üje F3 THEN 


36 
LET wieward 
838 FOR km@ TO 2 
LET X@eINT (X1+324CO8 (ul 


88 
Roder +h) 
LET WRZINT (y1+s@#SIN (Wi 
+K/RRPIJ éh) 
jee BLT K1,yA: DRAW X2-XL,yB- 
y 
318 LET x1ex@: LET vley@ 
323 NI k 
330 PLOT X1,yl: DRAW axxl‚av- 
y 
340 LET xX1räx: LET yizay 
358 NEXT j 
3588 IF INKEY$="" THEN GO TO 368 
378 STOP 
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w=90, 51=30, dw=48, 52=70 en n=15 geeft de onderstaande tekening. 
Als dw een deler is van 360, is n gelijk aan 360/dw. Kiest u voor 
dw een willekeurige waarde, neem dan voor n een groot getal, bij- 
voorbeeld 100, en breek het tekenen met de stoptoets af. De 'kinde- 
ren! doen dat in LOGO ook op deze manier. 





Vertaling LOGO-programma nr.2 in BASIC 
Het programma lijkt sterk op het vorige en behoeft daarom geen 


commentaar. 


188 Ren proacam 32 


ma 

rtle grafiek (LOGO-2 
ENT “coordinaten ter tpun 
vraa, Ensor” "peginrichting (grade 
138 iNeuT “yvêrplaatsi “is 
Ae INPUT “draaihoek Nlinksòn 
wrd 
52 INPUT “toenane zijde "ds 
162 LET ha0.5: LET rdal 181/186: kh 
Er siawdrd 


„san LET X2=INT (X4+$#C08 (wi) th 
190 LET Y2sINT (yl+saSIN Wi) th 
Re22 IF X@c2 OR X2>255 OR y2O O 

Sales THEN GO TO &6 
PLOT x21,yl: DRAU x2-x1,ya- 


LET Xlax2: LET 
ies EET beurd: IF vake THEN 


240 LET wisward: LET sustds 
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Het is handig de ingevoerde waarden na het tekenen op het scherm 
of op een printer te laten afdrukken. In de volgende tabel zien we 
de diverse waarden voor de daaronderstaande tekeningen. 











Tekening x1 yl 1 w T s dw ds 
A 128 88 90 5 144 3 
B 128 88 90 5 123 2 
C 128 88 90 7 92 2 
D 128 88 90 3 72 1 
E Ï 128 140 90 130 145 0 














7 
Sa 





hed 
NN 
SN, 


LS 4 
/ 
ijf 
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Vertaling LOGO-programma nr.3 in BASIC 


182 REM proprennn 33 
vierkantspiraal (LOGO-3) 
118, INPUT "coordinaten startpun 


‘ kut? 

a ee “beginfichting (grade 
120, INPUT “zijde beginvierkant 
„48 INPUT "draaihoek ‚Linksom 
‚dw 


189 INPUT “toename zijde “;ds 
160 LET he0.5: LET rdaPI/180: L 
ET wlewerd: LET vlag=® 

179 CL$ 


188 LET axaxl: LET ayzyl 
190 FOR fu@ TO PI STEP PI/2 


200 1802 

219 IF vlagel TMEN GO TO 396 
zee PLOT xl,y1: DRAW x@-x1,ya- 
230 


HEE xl=x@: LET viva 
RSD PLOT Xl,y1: DRAW AX-Xl,ay-y 


260 LET x1max: LET wima 

270 LET wewtdw: IF wen3b2 THEN 
LET waw 384 

89 LET Uimwtrd: LET Sus$+ds 
298 T 


GO TÔ 19 
302 IF INKEY&a"" THEN GO TO 308 
318 STOP 
ae LET X@aINT (xl+s#COS (wief) 
+ 
ze LET Y@INT (yA+S#SIN (wier) 
+ 
1829 IF X2<@ OR XR)2S5 OR VRO O 
R ya2>175 THEN LET vlagal 
1238 RETURN 
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Voor de boven het programma staande tekening kozen we x1=100, 
y1=120, w=90, 5=5, dw=15 en de=3,5. 


Vergelijk deze drie BASIC-programma's eens met de overeenkom- 
stige LOGO-programma's. In BASIC moeten we zelf alle graphic- 
software schrijven, terwijl deze in LOGO als machine-routines aan- 
wezig is. 


Het 4e en 5e LOGO-programma 


Tot slot van dit hoofdstuk geven we nog twee BASIC-programma's 
waarmee de turtle behalve rechte wegen ook kromme wegen kan 
bewandelen. Dit voegt een nieuw element aan de turtle-graphics 
toe. De LOGO-structuur zullen we stapsgewijs verfijnen. Als uit— 
gangspunt nemen we de onderstaande figuur. Deze bestaat uit een 
lijnstuk L, met daaraan vast een cirkel met straal R. We noemen 
deze figuur voor het gemak even 'de steelpan'. 


Als de turtle de steelpan tekent, begint hij in A. Hij legt vervol- 
gens een afstand L in positieve x-richting af en komt in B. Hier 
draait hij 90° met de klok mee en legt daarna de omtrek van de cir- 
kel af tot hij weer in B uitkomt. Daar draait hij zieh 90° met de 
klok mee en kijkt dan rechtvooruit naar het punt A. Met deze 
steelpan kunnen we twee soorten tekeningen maken. Laten we 

deze mogelijkheden bekijken: 


LOGO-programma 4 


TO CIRKELFIGUUR-1:LIJNSTUK:STRAAL:HOEK 
STEELPAN:LIJNSTUK:STRAAL 

LEFT : HOEK 
CIRKELFIGUUR-1:LIJNSTUK:STRAAL:HOEK 
END 
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LOGO-programma 5 


TO CIRKELFIGUUR-2:LIJNSTUK:STRAAL:HOEK 
REPEAT 4 (STEELPAN:LIJNSTUK:STRAAL LEFT 90) 
LEFT :HOEK 
CIRKELFIGUUR-2:LIJNSTUK: STRAAL: HOEK 

END 


Deze twee LOGO-programma's kunnen nog niet uitgevoerd worden, 
want de procedure STEELPAN is nog niet in het LOGO-systeem 
gedefinieerd. We gaan de structuur verfijnen: 


TO STEELPAN:LIJNSTUK:STRAAL 
FORWARD:LIJNSTUK RIGHT 90 
CIRKEL:STRAAL 

RIGHT 90 

END 


Ook in deze procedure zien we een nog niet gedefinieerde proce- 
dure, CIRKEL, opduiken. Om de turtle een cirkel te kunnen laten 
maken zullen we deze procedure CIRKEL moeten maken. Veel LOGO- 
versies bevatten reeds zo'n voorgedefinieerde CIRKEL-procedure. 


Wij gebruiken in het navolgende programma de CIRCLE-functie van 
de ZX Spectrum. In de regels 230 t/m 260 zien we dan ook: 


230 LET ax=x2: LET ay=y2 

240 LET U=INT (axtr*COS (w1) th) 
250 LET v=INT (aytr&SINCW1)+h) 
260 CIRCLE u‚v‚r 





Vertaling LOGO-programma nr. 4 in BASIC 





Vertaling LOGO-programma nr.4 in BASIC 


198 REM Edet 34 
cirkelfiguur-l (LOGO-4) 
118 CLS 


128 INPUT "coordinaten startpun 
zaal ENBUT “peginrichting (grade 
Eer ENPUT “Lengte van de steel 
158 INPUT "straal van de cirkelt 
266, INPUT “draaihoek , Linksom 


178 LET he8,5: LET rdaPI/18@: L 
ET uiaward 

182 CLS 

192 LET xxexl: LET Wymyd 

223 LET X@RINT (x14l#COS (WW) +h 


212 LET YRRINT (y2+L&SIN (wi) th 
228 PLOT x1l,yi: DRAW X2-XL,WR-y 


1 

230 LET axumxe: LET ayuyaä 
240 LET UeINT (ax+r&G03 (with) 
85 LET VaINT (ay+rasIN (wiek) 
258 CIRCLE u‚v‚r 

273 LET wavedu: IF w>e36O THEN 
LET waw-380 

260 LET wiawerd: LET Xizax: LET 


EEE 
83 IF Xlaxx AND yimyy TMEN GO 
To ® 

GO TO & 


308 ee 
312 IF INKEY&e"" THEN GO TO 918 
328 STOP 


Ga 
e 


De bovenstaande tekening krijgen we door x1=128, y1=120, w=90, 
1=80, r=20 en dw=170 te kiezen. 
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Vertaling LOGO-programma nr.5 in BASIC 


AOOTREN RIRREt IIP 2 (LOGO-5) 


128 SREUT “coordinaten startpun 
130 INPUT “"beginrichting (grade 
er 
140 INPUT "Lengte van de steel 
168, INPUT “straal van de cirkelt 
166, INPUT “draaihoek , Linksom 
158 LET ha0.5: LET rdsP1/180: L 
ET wizwird 
189 CLS 
199 FOR Jai 
06 "LET REnINT ocielacos (wie 
RLO LET VRRINT (WA+LKSIN WIJ+ 
22O PLOT X1,yi: DRAW KR-X1, ya 


232 LET axsxa: LET a 
248 LET UzINT (axers 58e twi) +h 


259 LET vaINT (ayer#SIN (Wi) th 


) 

did CIRCLE u,v 
27 LET presse: Car W°368 THEN 
Er Waw 3! 

258 ger” Samwerd: LET X1isax: LE 


Tad rk 


XT j 
340 LET wvawrdu: IF W>n36D THEN 
LET waw -36 
319 LET, et 
328 GO TO 
338 STOP 





Voor deze tekening geldt: x1=128, y1=88, w=0, 1=40, r=10 en dw=45. 
Met andere waarden voor deze variabelen kunt u de mooiste turtle- 
plaatjes maken. 











7 Educatieve toepassings- 
programma’s 


De voorgaande 35 grafische programma's zijn geen toepassingspro- 
gramma's. Ze waren bedoeld om een overzicht te geven van de 
mogelijkheden van graphics met hoog oplossend vermogen. In dit 
hoofdstuk zullen we vijf praktijkgerichte grafische programma's 
presenteren. Deze programma's zijn: 


Tekenen van een landkaart 

Maken van een histogram (een stavengrafiek) 
Demonstratieprogramma voor de breking van lichtstralen 
Demonstratieprogramma voor de 'speldenworp' van Buffon 
Prooi-roofdierpopulaties 


ar co ee 


1. Tekenen van een landkaart 


Educatieve programma's over topografie gebruiken dikwijls land- 
kaarten van een land of van een werelddeel die door de computer 
getekend worden. We zullen laten zien hoe de computer zo'n land- 
kaart kan tekenen. We gaan de landkaart (althans de grensomtrek) 
van Zwitserland tekenen. Waarom Zwitserland? Wel, in de eerste 
plaats omdat de auteur dit gekozen heeft en in de tweede plaats 
omdat Zwitserland geen eilanden heeft. Het tekenen van de contou- 
ren van Nederland gaat in principe net zo, alleen kosten al die 
eilanden een hoop extra werk, vandaar! 


Hoe tekenen we nu de omtrek van een bepaald land? Heel eenvoudig! 
We pakken gewoon een atlas, een stukje overtrekpapier, een pot- 
lood, een lineaal en millimeterpapier. We kiezen in de atlas een land 
of werelddeel uit en trekken het op overtrekpapier over. Daarna 
leggen we het overgetrokken plaatje op millimeterpapier, waar we 
van te voren een coördinatenstelsel op getekend hebben (een x- 

en een y-as). We bepalen nu van een (groot) aantal punten op de 
omtrek van het land (of werelddeel) de coördinaten (x,y) en schrij- 
ven deze op. Als we ons op het beeldscherm ook een coördinaten- 
stelsel voorstellen en als we hierin de punten, waarvan we de coör- 
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dinaten in een programma opnemen, met elkaar laten verbinden, dan 
ontstaat een 'landkaart' op het scherm. 


In het onderstaande programma, dat 'de kaart! van Zwitserland 
tekent, is een ‘echte! kaart gebruikt met een schaal van 

1 : 2.000.000. Om een enigszins natuurgetrouwe weergave te ver- 
krijgen zijn de coördinaten van 90 grenspunten berekend. De coör- 
dinaten van deze punten, die opgegeven zijn in millimeters ten 
opzichte van de oorsprong van het gekozen coördinatenstelsel, zijn 
in DATA-regels opgenomen. Als het programma het coördinaten 
paar (0,0) leest, weet het dat de 'kaart' af ìs. 


De x-coördinaten liggen tussen 6 en 177. Om deze naar het hoge- 
resolutiebereik 0-255 te transformeren vermenigvuldigen wij zowel 
de x- als de y-coördinaten met de factor k = 1,3. Hierdoor blijft 
nog wat ruimte over om een kader rond de kaart te tekenen. De 
werking van het programma zal hiermee duidelijk zijn. 


We hebben op deze manier ook een kaart van Europa en zelfs een 
wereldkaart getekend. Het probleem met de eilanden hebben we als 
volgt opgelost. Als het programma in een DATA-regel negatieve x- 
en y-coördinaten leest, weet het programma dat dit punt niet met 
het vorige verbonden moet worden en dat dit punt dus het begin 
is van een apart stukje 'land'. 






ZWITSERLAND 











128 REM 
118 
128 

=i5: 
198 
142 

12 IN 

5 


158 
@zINT ( 


Tekenen van een landkaart 
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aart van Zwitserland 
„1,3: LET he8.5: LET u 
z175 

Eainr (uekaxsh): LET u 
EN] 

ZEN tuekaxth): LET U 


+ 


„iee PLOT en Deur: DRAW XBRL, 
178 BER, LET yiny@ 


200 IN 
219 STOP 
220: 

323 DATA 
5,104, 75, 1 
$19 bATA 
1,108, 94, 1 
420 GATA 
8,106, 118, 
439 SATA 
‚112, 112,1 
340 DATA 
2,145, 123, 
EEEN 

’ 
560 GATA 
74,59, 177, 
358 ATA 
e1,Se, 165, 
360 ÓATA 
43,38, 199, 4 
384 _ÓRTA 
22 29,125,1 
469 SATA 
16,25, 108, 
26180 
á&a ÓÂTa 
248,13 
433’ DATA 
3,29,45 
448 DATA 
‚6,3a 
463 DATA 
3,25,63 
469'DATA 
«46, 102 
470 /oATA 
102,89, 108 
486 DATA 


138: K0 THEN oo TO 16 


KE Barn THEN GO Po zoe 


328, 7j 107,70, 104,7 
1185781, 104,86, 105,9 
1,196,109,196, 195,10 
il8,102,117,108,118 
je, 126, 110,199, 10 
:98r144 70,14, 77,1 
,68,199,76,177,79,1 
7,52,171,51,167,56,1 
86,94, 162,36,187,42,1 
ige, 45,193,48, 132,40, 1 
122,11,119,12,114,20,1 
102,45, 94, AR, 05,5, DO, 
75,15,66,19,60,14,52,1 
37,29,99,36,37,40, 99,4 
16,38,18,33,13,29,6,28 
11,34,13,40,19,48, 12,5 
26,73,39,79,48,94, 42,9 
54,102,53,99,51,98,89, 
2, 


oe 
a 


Ë 
] 
6 
8 
C) 


BRbRGeNeer GES 0S0 
GO IRBOP-P POOIER O 
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In dit programma leggen we de oorsprong (u,v) linksbovenaan het 


scherm; u=15 en v=l75. De coördinaatparen x,y zijn namelijk door 
de auteur ten opzichte van een oorsprong in de linkerbovenhoek 


berekend. 


Kijk in de appendix voor een uitbreiding van dit programma. 
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2. Maken van een histogram 


Vaak willen we een aantal waarnemingen grafisch in een stavendia- 
gram of histogram weergeven. Hiertoe dient het volgende programma. 
Dit programma tekent een horizontale en een verticale as. De verti- 
cale as wordt in stukjes van acht schermpuntjes verdeeld. Dit komt 
overeen met 5% van de hele verticale as. Elk tweede streepje op 

deze as geeft de volgende 10% aan en is iets breder getekend. Hier- 
mee kan de lengte van de staven redelijk geschat worden. 


Het programma kan maximaal 100 waarnemingen verwerken, Uit este- 
tische overwegingen moet u echter niet veel meer dan 40 waarnemin- 
gen invoeren, omdat de staafjes anders meer op lijntjes dan op balk- 
jes gaan lijken, 


De waarnemingen worden zo 'geschaald! dat de grootste waarneming 
overeenkomt met de lengte van de verticale as. 


Het zou mooi zijn als we bij de staven, de assen en de schaalverde- 
ling tekst en getallen konden zetten, maar dat gaat heel lastig, dus 
doen we het nu maar niet, Denkt u erom dat de getallen die u 
intoetst de lengte van de staven aangeven. Als elke staaf een 
bepaalde klasse met waarnemingen voorstelt, dan voert u dus steeds 
het aantal waarnemingen (de frequentie) van een klasse in. Het pro- 
gramma kan uitgebreid worden door er een stuk voor te zetten dat 
ruwe gegevens inleest, die vervolgens netjes in klassen verdeeld 
worden, De frequentie van de waarnemingen in de klassen wordt 
vervolgens gebruikt om de stavengrafiek te tekenen. 
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Maken van een histogram 


198 REM programma 37 histogram 
118 CLS 
128 PRINT SHISTOGRAN TEKENEN, 
A39 PRINT "e=meememeeveven _ 
148 PRINT ; PRINT 
pe INPUT “Hoeveel gegevens (4 
EN 
î 
158 DIM an) 
178 LET mxs-1e1@: PRINT 
188 FOR Jel TO n 
198 PRINT "Waarde “;j;TAB (12) 
208 INPUT aij) 
205 PRINT aij) 
zere IF a(j)>mx THEN LET mxmalj 
228 NEXT j 
238 CLS 
240 REM horizontale as 
859 PLOT 2,18: DRAW 255,9 
280 REM verticale as 
278 PLOT 19,10: DRAW 2,165 
889 REM schaalverdeling 
298 FOR jal TO 1 
300 ET Xis4: LET yimj#16 
318 LET x&=1û: LET 


yeeyl 
gie PLOT XA,yl: DRAW x@-xl,ya- 
zeeer X1la7: LET ylay@+8: LET 
yv 
gen PLOT X1,yl: DRAW XR-XL,YR- 


y 
358 NEX 
Pae REN staven tekenen bubreedt 


57e LET beINT (36/0): LET he@, 


380 FOR jet T 
382 LET ziatg. =D eva15: LET yle 


EN 

403 LET xasXl: LET V@sINT (10+ 
165xa lj) ZmXeh) 

418 LOT X1,yl: DRAW xR-X1,V2- 


yi 
420, LET zlex8: GET 9iev2: LET 
Xänxie+b: LET vas 

85e "DAR vari gur 

„iso, LET vlet: Er Saave: Ler 
= 

<50 DRAU x2-x1,va-v1 


482 
478 bid INKEY Se" THEN GO TO 478 
489 STOP 


99 
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3. Demonstratieprogramma voor de breking van lichtstralen 


Met dit programma willen we laten zien hoe we de ZX Spectrum bij 

natuurkundelessen zouden kunnen gebruiken, Voordat we het pro- 
gramma geven leggen we nog iets uit van de natuurkundige begin- 
selen van de breking van licht. 


Als een lichtstraal vanuit de ene stof (bijvoorbeeld lucht) een ande- 
re, optisch dichtere, stof (bijvoorbeeld water) binnenkomt, wordt 
de straal op het scheidingsvlak van beide stoffen naar de verticale 
as toe gebogen. Zie onderstaande figuur. 


[Verticale Werticale 
jas as 

Medium 1 | 
(lucht) ij 









Medium II 
|__ water) 


Hierbij geldt de brekingswet van Snellius: 





ej en ez zijn de lichtsnelheden in respectievelijk de eerste en de 
tweede stof. De constante n heet de brekingsindex. Bij de overgang 
van lucht naar water geldt een brekingsindex van 1,33. Voor de 
overgang van lucht naar glas gelden andere waarden, enzovoorts. 


Als het licht vanuit een 'dichter' medium overgaat in een 'dunner' 
medium geldt het omgekeerde: de lichtstralen worden nu van de 
verticale as, die loodrecht op het scheidingsvlak van beide stoffen 
staat, afgebogen. De brekingshoek a kan nooit groter dan 90° zijn. 


Voor dit grensgeval (a = 90°) geldt: 


sin90° 


8 in90e = 
sing (want sin 90 1) 


=n > singt= 





dE 
n 
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Bij de overgang van lucht naar water geldt voor B*: 
ere AL 5 En Rr, el 
sing = 155 >= sinp*=0,7519 > @ 48,75 


Groter dan 48,75° kan @ dus niet worden, Dit betekent dat als het 
licht van water overgaat in lucht met een invalshoek a die groter is 
dan gf (zie rechter figuur), het licht niet meer het water ‘uitkomt’. 
Op het scheidingsvlak van water en lucht wordt het licht geheel 
teruggekaatst. Op dit principe berusten de moderne glasvezelkabels, 
waarin informatie in de vorm van licht wordt getransporteerd. 


In het onderstaande demonstratieprogramma kan de brekingsindex 
n ingetoetst worden. De lichtbron bevindt zich in het punt met 
schermcoördinaten (0,10). Het scheidingsvlak ligt horizontaal en 

is de lijn v=88. De invalshoek van een lichtstraal die van medium 2 
(de dichtere stof, onderste helft) in medium 1 (de lichtere stof, 
bovenste helft) overgaat wordt van 0° steeds met stapjes van 3° 
opgehoogd. Op het moment dat deze invalshoek groter wordt dan @* 
(deze is afhankelijk van de ingetoetste brekingsindex) treedt totale 
reflectie (terugkaatsing) op. Rond de figuur wordt een kader gete- 
kend. De eindpunten van de diverse lichtstralen worden met trigo- 
nometrische functies berekend. De commentaaropdrachten in het 
programma leggen nog eens uit ‘wat waar! gebeurt. 


De inverse functie van de sinusfunctie (de boogsinus of arcsinus) 
bestaat niet op de ZX Spectrum. We lossen dit op door de inverse 
tangensfunetie (ATN in BASIC) te gebruiken, 





Ge 
1 





Het voordeel van een dergelijke computersimulatie, boven het uit— 
voeren van het natuurkundige experiment, is dat we heel gemakke- 
lijk verschillende brekingsindexen kunnen invoeren zonder steeds 
andere apparatuur op te hoeven stellen en dat de lichtstralen als 
dunne lijnen zichtbaar gemaakt kunnen worden. 
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190 REM Brofane 38 
king van Licht 








aL 

128 SREur brekingsindex n ";n 
130 LET vw88: LET he@,5: LET rd 
=PI/180: LET bed 

148 PLOT B,v: DRAW her 

158 PLOT @,1/ DRAW 255,8 

188 DRAU 2,185: DRAW -266,2: DR 
AW_@,-15 

173 REM stralen in medium 2 en 
2 tekenen 

188 REM kende in graden; blabe 
ta in Cadi 

199 REM tlas lpha in radialen; s 
azsIN (a 

220 LET brb+3: LET bleberd: LET 
xls: LET iz 

21e ed zatint UZIETAN (bi) +h): 


bi 
220 EE xa255 JAN GO TO 3978 
238 PLOT X1,y4: DRAW xRexi,yR-y 


1: LET Xlsx2: LET yi=y2 

240 REM SIN (alphâ) en alpha be 

rekenen 

eend REM controle op totale reft 

stie 

863 LET san&SIN (bi): IF $>1 TM 

EN GO TO 338 

270 LET AlxATN [S/SOR (les4s)) 

880 LET XBEINT (X14+774TÄN (al)+ 

hì: LET yaa165 

292 IF XA<255 THEN ee zee yi: 

DRAW XR-XL,yR-yl: GO TO 

300 LET Xa=256: LET MEM ive t 

255-x1) TAN (ai) +hl 

318 PLOT xX1,y1: DRAW XR-xL,yBey 

1: 60 TO 298 

320 REM totale EEN SBC ETE 

338 LET Kandeel ve 

340 IF Xac255 TMEN PLOT x1,ui 

PRAIJ Xa -X1 Dh Go TO 

350 LET xt=255: LET YBSINT wel 

255-x1) /TAN (bi) +h} 

1058 PLOT X1,Wi: DRAW Hi=Xl,yaey 
60 TO 288 

‘520 IF INKEY&a"" THEN GO TO 378 


388 STOP 
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4, De speldenworp van Buffon (1773) 


Stelt u zich voor dat we in een plat vlak een aantal evenwijdige 
lijnen met een onderlinge afstand a tekenen. 


We werpen nu, zonder echt te ‘mikken’, een speld met een lengte c 
die kleiner dan of gelijk aan a is (c<a) op het vlak met de evenwij- 
dige lijnen. Hoe groot is nu de kans dat een speld een van de lijnen 
treft, dat wil zeggen snijdt of raakt? 


/ 


speld snijdt een lijn speld raakt een lijn 





Dit probleem kwam in 1773 op bij de Franse wiskundige Buffon na 
het zien van de Amerikaanse vlag met de ‘stars en stripes'. Buffon 
heeft berekend dat deze kans gelijk is aan 

2e 


an 


Omdat deze formule de constante n bevat, heeft men deze 'spelden- 
worp' vaak gebruikt om de waarde van n door simulatie te bepalen. 
We werpen hiertoe vaak, bijvoorbeeld een lucifer, op een papier 
waarop een aantal evenwijdige lijnen (met een onderlinge afstand 
die groter dan of gelijk aan de lengte van de lucifer is). Als de 
lucifer in n worpen k keer een lijn treft, dan geldt dat 


2e 


ongeveer gelijk is aan k 
an n” 


Beide uitdrukkingen geven namelijk een indruk van de kans dat 
de lucifer een lijn treft. 
We kunnen als volgt de waarde van n benaderen: 


ze _k zen 
Tl n 


> ne 
ak 
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Op deze manier heeft de astronoom Rudolf Wolf in 1850 met 5000 luci- 
ferworpen voor n de waarde 3,1596 gevonden (n = 3,141592654....). 


Met een computer kunnen we gemakkelijk duizenden worpen sìmule- 
ren. Deze programma's vinden we in bijna elk informaticaleerboek. 
Dit zijn echter programma's die vrij lang draaien, maar waaraan 
niets te 'zien' is. Als we bijvoorbeeld in zo'n programma voor n de 
waarde 10.000 zouden intikken, zien we eerst een tijd niets en ver- 
volgens verschijnt de mededeling dat van de 10.000 keer 6.349 keer 
een lijn is geraakt, hetgeen voor n de waarde 3,1501024 oplevert. 


We gaan een programma maken dat ook dergelijke berekeningen uit- 
voert, maar dat bovendien elke speld die geworpen wordt laat zien. 
We zien het experiment als een film aan ons voorbijgaan. Bekijk 
hiertoe de onderstaande tekening. 














o speld AB Is 55 lang 


De manier waarop een speld valt kan met twee toevalsgetallen 
bepaald worden. De twee toevalsgetallen bepalen respectievelijk 
de afstand ym van de speld tot de daarboven liggende lijn en de 
hoek w die de speld met de lijnen maakt: 


Osymza ; ym is de afstand van het midden van de speld 
tot de daarboven liggende horizontale lijn 


Osws180° ; wis de hoek (in positieve zin, dat wil zeggen 
bij draaiing tegen de klok in) die de speld met 
de positieve x-richting maakt 


Hierna volgt een illustratie van hoe het eruit kan zien en het pro- 
gramma. 
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109 REM programa A9 
Bhetdenworp van Bufion 


CLS 
120 PRINT "SPELDENWORP VAN BUFF 
130 PRINT "e=emmmemenmmmenenveenn 


149 INPUT "Hoeveel worpen “in 
150 LET med: LET h2@.6 


CLS 
170 REM Lijnen trekk 


en 
188 FOR 4-8 TO 178 STEP 56 
192 Plof B,y: DRAW 255,0 
239 NEXT vw 
219 REM n naal gooien en tekent 
n 
223 FOR Jel TO_n 
232 LET XMaINT (195&RND+3@+) 
B4@ LET YM=INT (SB+RND+63+h) 
250 LET waPExRND 
8250 LET dxun?.54COS (w): LET d 
YER? .S4SIN Cw) 
270 LET XIRINT (Xm-dXth): LET 
VÄnINT (jm dyn 
880 LET XBRINT Oxmedxeh): LET 
YÄINT Ek BE 
vee PLÔT x1,yi: DRAW XR-XL,VB- 
30e IF y1icsb@ OR yRr=115 THEN 
LET Matt 
312 NEXT J 
380 IF INKEY&e"" THEN GO TO 320 
338 CLS 
ziee ERINN “aantal worpen";TAB ( 
556 PRINT "aantal keer snijden 
“;TAB (21); m 
A60 PRINT “benadering voor Pi“; 
TAB (24); aen /M 
370 STOP 





Uit de waarden voor ym en w berekenen we de coördinaten van de 
uiteinden A en B van de speld. Het programma tekent hiermee de 
speld op het beeldscherm. 


Een speld treft een lijn als de y-coördinaat van A groter dan of 
gelijk aan 115 is of als de y-coördinaat van B kleiner dan of 
gelijk aan 60 is. De x-coördinaten van A en B doen er in het 
geheel niet toe. 


Het programma tekent als Ispeelveld' vier evenwijdige lijnen met 
een onderlinge afstand van 55. De spelden zijn trouwens ook 55 
lang. Als alle spelden geworpen zijn kunnen door een toets in te 
drukken de waarden voor n, k en n worden afgelezen. 
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Verander de regels 320 t/m 370 als volgt: 


320 PRINT AT 18,1;"aantal worpen";TAB(21);n 

330 PRINT AT 19,1;'"'aantal keer snijden'';TAB(21) ;m 
340 PRINT AT 20,1;"benadering voor pi';TAB(21); 2*n/m 
350 IF INKEYS="""" THEN GO TO 350 

360 STOP 


Nu worden de gegevens over het aantal worpen en het aantal keer 
snijden en over de benadering van pi onder de tekening afgedrukt. 


5. Prooi-roofdierpopulaties 


Het laatste educatieve programma is een (deterministische) simulatie 
van een ecologisch systeem. Het is een demonstratieprogramma voor 
een biologieles. 


Het ecologische systeem bevat gras, hazen en vossen. Tussen deze 
drie ecologische componenten gelden de volgende betrekkingen: 


1. De hazen eten gras en de vossen eten hazen. 

2. Als er meer gras groeit, neemt ook het aantal hazen toe. 
Deze hazen eten echter van het gras en verminderen zo hun 
eigen groei. 

3. Als er meer hazen komen, neemt ook het aantal vossen toe. 
Omdat vossen hazen eten, verminderen zij zelf hun groei, net 
zoals bij de hazen en het gras. 


Als we het aantal hazen op een bepaald tijdstip t aangeven met h(t) 
en het aantal vossen met v(t), kunnen we, volgens Lotka en 
Volterra (1920), voor het aantal hazen en vossen op tijdstip t+1 de 
volgende vergelijkingen opstellen: 





h(t+1) = h(t) +a.h(t) - b.h(D).v(t) vergelijking (1) 
v(trl) = v(t) +e.v(t).h(t) — dev(t) vergelijking (2) 


De toename van het aantal hazen tussen de tijdstippen t en t+1 is 
evenredig met het aantal hazen op tijdstip t‚ dus 


h(t+1)-h(t) =a.n(t) 


toename hazen aenei hazen op tijdstip t 


groeifactor 
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De afname van het aantal hazen is evenredig met het aantal aanwe- 
zige hazen (natuurlijk verloop) en met het aantal vossen, want die 


eten hazen, dus 


h(t+i)-h(t) =a.h(t) - b.h(t).v(t) 





groei van \ afname 
het aantal hazen jj aantal hazen 
toename 
aantal hazen 


Dit is vergelijking (1). We nemen in het programma aan dat er altijd 
genoeg gras voor de hazen voorhanden is. 


De toename van het aantal vossen is evenredig met het aantal aan- 
wezige vossen en met het aantal hazen (prooi). De afname van het 
aantal vossen is alleen evenredig met het aantal, omdat in dit 
systeem de vossen zelf geen prooidieren zijn. Voor de vossen krij- 
gen we dus: 


v(te v(t) =e.v(t).h(t) - d. v(t) 





groei van id afname 
het aantal vossen i aantal vossen 
toename 
aantal vossen 


De vergelijkingen zijn als zogeheten differentievergelijkingen opge- 
steld. Het gras speelt, zoals gezegd, eigenlijk geen rol; er is altijd 
genoeg om alle hazen te voeden. 


Bij het draaien van het simulatieprogramma hebben we de volgende 
waarden gekozen: 

x (= h(0)) = 200; y (= v(0)) = 20; a=0,3; b=0,01; e=0,002 en 
d=0,5. 


De grafiek die het programma tekent laat heel mooi de groei en de 
afname van de hazen- en vossenpopulaties zien. Het aantal hazen 
groeit eerst, bereikt een maximum en neemt vervolgens af. De vos- 
sen groeien ook, maar later, bereiken later een maximum en nemen 
dan ook af, waarna de eyclus zich herhaalt. In de biologie noemen 
we dit een dynamisch evenwicht. 


U hoeft de waarden voor de coëfficiënten a, b, c en d maar iets te 
veranderen of het systeem kan ontregeld worden. Hierbij neemt òf 
het aantal hazen enorm toe òf alle hazen en vossen sterven snel uit. 
Deze eenvoudige simulatie toont aan hoe desastreus een kleine 
ingreep in een bestaand ecologisch systeem dat in een dynamisch 
evenwicht is kan zijn. 
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hazen 


vossen 








182 nma 4@ 

OOidiër populatie 
119 
120 inpopuletie prooi 
dier ix 
„422 inpopulatie roofd 
148 tifactor prooidie 
ren 
152 amefactor prooidi 
eren 
zie tifactor roofdier 
n_{ 
178 amefactor roofdie 
ren 
189 
ad @.175: DRAW @,-175: DR 
Ee ET KuO,3: LET he@.5: LET v 
= 
298 FOR Used Ae 247 STEP 8 
210 LET xiw 
220 LET ke Ke Caere hdd. 
230 LET eidipe ep =d ty) 
„2%0 REM pop.prooidieren tekene 
250 LET YIaINT (vekexe+h) 
En nld Xaex1+B: LET Y&BINT (v 
2® De W@B OR V2175 THEN GO 
TO 368 


vie PLOT x,y: DRA xRexl,ya- 
299 REM pop, roofdieren tekene 


38 LET Närklrer LET GBeinT ov 
a8 ‘Ie Y2<@ OR „2175 THEN GO 

38 PLOT X1,y1: DRAU YR-X1,Va- 
LET Xaxp: LET veur 


NEXT u 
If INKEYüz"" THEN GO TO 380 
STOP 


2 








Appendix 


In deze appendix laten we zien hoe een aantal programma's uitge 
breid kunnen worden om leuke effecten te bereiken. Ook geven 
we nog een aantal tekeningen die u met de programma's kunt 
maken. Tot slot volgt dan nog een betere versie van het bolpro- 
gramma 25. 


Programma 4 


In de tekening, die programma 4 oplevert, zien we dat tussen een 
zeshoek en de daarbinnenliggende zeshoek zes driehoeken ontstaan. 
We gaan deze driehoeken om en om ‘inkleuren’. We doen dit niet 
alleen voor de driehoeken tussen de buitenste zeshoek en de eerste 
ingeschreven zeshoek, maar ook bij de driehoeken tussen de eerste 
ingeschreven zeshoek en de tweede ingeschreven zeshoek, voor 

de driehoeken tussen de tweede ingeschreven zeshoek en de derde 
ingeschreven zeshoek, enzovoorts. Er ontstaan dan drie ‘zwarte! 
spiraalvlakken. 


Als we een vlak, in dit geval een driehoek, dat geheel omsloten 
wordt door lijnen willen inkleuren, moeten we in dat vlak een punt 
lokaliseren. De schermcoördinaten van dat punt gebruiken we dan 
in de ‘paint-routine' om het vlak waarin dat punt ligt te kleuren 
(zwart of een andere kleur). Omdat er een groot aantal driehoeken 
gekleurd moeten worden, moeten we ook een groot aantal coördina 
ten van punten in die driehoeken berekenen. Deze punten moeten 
op een systematische manier in de FOR-lus (regels 220-350, p.115) 
bepaald worden. Als we in het X-Y-vlak een driehoek ABC tekenen 
waarvan de coördinaten van de hoekpunten (Xa,Ya);(XB:Yg) en 
Xe,Yo) zijn (zie de tekening op p.110), dan weten we uit de 
meetkunde dat de coördinaten van het zwaartepunt (Z) van de 
driehoek gelijk zijn aan 


Xa * XB *Xe Ya *YB* Ye 
{ 3 k 3 





)« 
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Het zwaartepunt van een driehoek is het snijpunt van de drie lijnen 
die elk een hoekpunt met het midden van de daartegenoverliggende 
zijde verbinden. In de tekening is Z het zwaartepunt en er geldt 
dus dat 
ASK ek Se 
Z 3 z 3 M 








XA Xz XC X5 





Kz 
Yz 


1/3 Xa+XB+Xo) 
1/3 (Ya+YBtYC) 


nu 


We kiezen in ons programma steeds het zwaartepunt van de driehoek 
als inwendig punt. Om de coördinaten hiervan te kunnen berekenen 
moeten we dus de drie x- en de drie y-coördinaten van de hoekpun- 
ten weten. We hebben in de figuur op pagina 111, die door pro- 
gramma 4 getekend wordt, één driehoekje als voorbeeld genomen. 
Van dit driehoekje zijn de hoekpunten A en B twee hoekpunten van 
de ingeschreven zeshoek, terwijl hoekpunt C een hoekpunt is van 
de daarvoor getekende (in ons voorbeeld de buitenste) zeshoek. 

We moeten dus steeds van twee zeshoeken de coördinaten van alle 
zes hoekpunten ter beschikking hebben, In het onderstaande pro- 
gramma hebben we, om de coördinaten van de hoekpunten van de 
zojuist getekende zeshoek te kunnen onthouden voordat de nieuwe 
berekend worden, twee arrays extra opgenomen, en wel p(7) en 
a(7). In regel 325 maken we p(j) en q(j) gelijk aan de coördinaten 
van de zojuist getekende zeshoek (x(j) en y(j)), vlak voordat x(j) 
en y(j) gelijk worden gemaakt aan de coördinaten van de hoekpun- 
ten van de volgende zeshoek (a(j) en b(j)) in regel 330). Zie p.115. 
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y 
Ys VC 
Yz F 
Eg 
EN 
= x(t) 
= xk) 
à = p(ke) 











(0,0) 


We zien vervolgens in de regels 263 en 264 hoe we de schermcoördi- 
naten van het zwaartepunt berekenen uit de twee hoekpunten van 
de nieuwe zeshoek x(k), x(k+1) en y(k) en y(k+1) en één hoek- 
punt van de daarvoor getekende zeshoek (p(k+1) en q(k+1)). In 
de FOR-NEXT-lus van regel 262 nemen we de stapgrootte 2 om 
steeds één driehoek over te slaan. Natuurlijk hoeven we dit inkleu- 
ren pas te doen als de tweede zeshoek getekend is, vandaar de 

IF n=1 THEN GO TO 270 in regel 261. Het inkleuren gebeurt ten- 
slotte in regel 265 en 266 met de opdrachten PLOT x,y en 
RANDOMIZE USR 64800. We leggen dit hierna uit. Na deze uitleg 
drukken we af hoe de ‘gekleurde! zeshoeken eruitzien, gevolgd door 
het hiervoor benodigde, aangepaste, programma 4. 


We hebben in het blad Your computer van oktober 1983 een routine 
gevonden voor het kleuren van een vlak, dat geheel door lijnen 
omsloten wordt. Willen we zo'n vlak kleuren dan moeten we de. gra- 
fische cursor met een PLOT-opdracht naar de coördinaten van een 
punt ergens binnenin dat vlak brengen. In programma 4 hebben we 
steeds als ‘inwendig! punt van een te kleuren (op te vullen) drie- 
hoek het zwaartepunt van zo'n driehoek genomen. Met PLOT x,y 
(regel 265 in het aangepaste programma 4) brengen we de cursor 
dus naar een punt binnenin zo'n driehoek. Met RANDOMIZE USR 
64800 roepen we een machinetaalprogramma aan dat we hiervoor 
speciaal in het geheugen gePOKEd hebben. Deze machinetaalroutine 
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punt (x,y) ligt in de dan geldende 


… In ons programma is dit zwart, maar wellicht kunt 


u deze kleur veranderen of zelfs drie spiralen van verschillende 


kleurt het vlak waarin het 
kleur maken, 


afdrukkleur 


Wilt u in uw grafische programma's ook vlakken vullen, dan zult u 


eerst de machinecode uit het onderstaande BASIC-programma in 


de BASIC-programma in te lezen en te draaien. Als alles goed gaat 


het geheugen moeten 'laden', U kunt dit doen door het onderstaan- 


vernietigt! het programma zichzelf en kunt u de, door dit program- 
ma in het geheugen, gePOKEte machinecode SAVEn. Als u vlakken 
wilt kleuren moet u eerst deze machinecode LOADen (zie p.114). U 


kunt dan altijd, als u dat wilt, aan een programma een PLOT- en 


een RANDOMIZE USR-opdracht toevoegen als u iets wilt kleuren. 








LEE 
paYafarsef 
HOSPES IEORELLESSKHEOPESEEEOSOEELLLELELLOEGOLSEELLEELEEESESELSLEELE 
AAD LDD ORO TUD HOON ORD IAD SOON DAD AAD HDO DAS AUD HDD DRO 
Acts NANANEBAAODDON DOOS td tt d tt 





Appendix 
ET_F=0:_DIM A L59) 
Bag DATA 7377,8807,8913,5421,63 
4D Bet EEEB, 4471,6418, 4636, 4563 
„9558, 8493 8634 8595, 3Â0I IDEE, 
1384,1284 423, Sa18 4938, 7058, 3 
225694 5489, 6622 
540 DATA 4770,4772,4771,S723,77 
52 9o27 10313, 1811,5509,6277, 634 
3e4d61, SE68 BES BTR, GIOS, 7O12, 
1ó55j5óba 1658 4da1 Sad 1180, 77 
840 RESTORE 529: LET Cr@: FOR I 
=ì To 50; RERD Air): LET GaC+ACI 
Dei: NEXT Is READ U; IF GU THE 
NOPRINT pERROR IN LINES 522 oR 5 
Eat FOR 1=19 TO 52 STEP 10: RE 
T Cs@: READ 


STORE I: LE 
625 PRINT I 
612 FOR Nei TO 19 STEP @: LET L 





628 ziL-48-(L257) #7) + (H-48 
Zn ee LET CeC+Uet (Ne) / 
638 POKE A,U: LET Azä+i: NEXT _N 
: IF C4>A(I/18) THEN PRINT "ERRO 
R_IN LINE ";I: LET Fal 

643 NEXT I 

658 IF Fed THEN PRINT "AL OK: B 
YE“: BEEP 1,0: POKE 23527,PEEK 2 
3635: POKE 23628 PEEK 23636: CLE 
AR 25828 

652 PRINT ''"CORRECT MISTAKES A 
ND RERUN. 
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Dit programma doet het zowel op een 16K als op een 48K Spectrum. 


Hebt u een 48K Spectrum, dan kunt u de machinecode naar een 
hoger adres verplaatsen, Hiervoor moet het volgende programma 


gedraaid worden. Denk erom: dit kan alléén voor de 48K Spectrum, 


maar het is niet noodzakelijk om de programma's samen met de 
paint-routine in dit boek te kunnen draaien. 


1 REM relocation pr am 

2 REM to conv B&R to 4Bk 
3 REM 

9 LET as32828 

19 FOR is32020 TO 32502 

Dn LET VaPEEK (i+i)+@56&PEEK U 


8 IF [PEEK iz4á@ OR PEEK i=194 
PEEK is=185 OR PEEK al OR 


OR 
PEEK in2d@ OR PEEK in204 OR PEEK 
i=295 OR PEEK im210 OR PEEK is 
12 OR PEEK im@18 OR PEEK ima2@) 
AND v>32800 AND vc32580 THEN POK 
E i+a,PEEK i: POKE i+i+a,lv+a)-2 
S64INT ((v+a)/256); POKE i+a+8,I 
AT GO) Z08B): LET isi+t@: GO tö 
30 POKE ita,PEEK 1 
43 NEXT í: POKE 64894,3: POKE 
64895,220: CLEAR 58600 
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Als u het programma met de paint-routine ingevoerd hebt en draait, 
zal het zichzelf vernietigen als het geen fouten in de DATAregels 
10 t/m 500 gevonden heeft. Het programma staat dan niet meer in 
het geheugen maar het machinetaalprogramma, dat gecodeerd in 

de DATAregels 10 t/m 500 staat, zit nog wel in het geheugen. U 
kunt vervolgens dit machinetaalprogramma SAVEn met de opdracht: 


SAVE "FILL" CODE 32000,500 


Hebt u op uw 48K Spectrum het ‘relocation program' gedraaid, dan 
moet u de machinecode SAVEn met: 


SAVE "FILL" CODE 64800,500 


Wilt u de paint-routine gebruiken dan hoeft u alleen maar het 
machinetaalprogramma in te lezen. U doet dit met: 


CLEAR 25800 : LOAD "FILL" CODE 
Op uw 48K Spectrum kunt u gebruiken: 
CLEAR 58600 : LOAD "FILL! CODE 


Mocht u de paint-routine in uw programma's gaan gebruiken dan 
moet u aan het begin van zo'n programma opnemen: 


CLEAR 25800 en bij 48K: CLEAR 58600 


We zien dit in regel 115 van het onderstaande programma 4, waarin 
we (op een 48K Spectrum) van de paint-routine gebruik maken. 
Het aanroepen van de paint-routine gebeurt dus door eerst een 
inwendig punt (x,y) van het te kleuren vlak te bepalen en met 
PLOT x,y de cursor ernaartoe te brengen; daarna kan de paint- 
routine worden aangeroepen. Dit kan met 


RANDOMIZE USR 32000 in een 16K Spectrum en 
RANDOMIZE USR 64800 als u een 48K Spectrum hebt. 


U kunt de routine ook aanroepen met LET A=USR 32000 respectieve- 
lijk LET A=USR 64800. 


De paint-routine volgt de eventueel gegeven INK-, BRIGHT- en 
FLASH-opdrachten. Zorg dat het te kleuren vlak geheel omsloten is 
door lijnen (er mag geen pixel in ontbreken) anders zal de routine 
het hele scherm vullen! 

In het onderstaande programma 4 ziet u hoe de routine wordt aange- 
roepen. Eerst geven we echter het resultaat van het programma. 
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100 REM programma & 
ingeschreven zeshoeken 
119 CLS 
115 CLEAR 58600 
188 DIM x(7) : DIM wi7) 
190 DIM aí7) : DIM Bi7) 
195 DEM PLA, DIM A7) 
140 LET Uwi2è : LET veeg 
150 LET r=88 : LET huô.5 
189 LET wa6OtPIJ1e0 
178 FOR _Jei TO 7 
189 LET wisjeu 
„298 LET x CENT twerkcos Wan 
“po LET y(JDSINT (vereSEN (Wi) 
+ 
210 Td 
B2O FOR nei TO 
230 FOR I= TO 6 
248 PLOT xj) ve 
‚25e DRAU XG kt ye -y 
89 NET} 
Sl EEneioTHEN GDI 270 
2 FOR kua TÔ 6 
263 LET XZINT (CX CKD +X (KID +P 
(kel 75) 
264 LET VeINT C(y(k)ew (ken) +a 
tk+19) 75) 
ZSE PLOT MAY 
25E RANDOMÍZE USR s4800 
RE7 NEXT k 
270 FOR kus TOG, 
289 LET a(KDeINT COX ex ther 
22+) 
298 LET BIKDGINT ((y thd ey tkel 
Beh) 
393 NEXT k 
‚gie LET alma : LET br eb 
320 FOR jat TO 7 
S25 LET POD: LET aj ey 
„5e LET xiJvzatj) : LET voDe 
J 
349 NEXT J 
359 NEXED 
68 IF INKEYGn"" THEN GO TO 360 
570 SToP 
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Programma 19 


We kunnen van de vliegekop, die programma 19 tekent, een ‘echte! 
vliegekop met ogen maken. Hiertoe tekenen we twee ogen (cirkels) 
en kleuren het vlak dat door de cirkels en de daarbuitenliggende 
lijnen begrensd wordt. Het resultaat en het gewijzigde programma 
19 staan hieronder. 

In de regels 221 en 222 tekenen we de twee ogen met de CIRCLE- 
opdracht. In de regels 223, 224 en 225, 226 kleuren we de vlakjes 
rond de ogen met de machinetaalroutine die we in het geheugen 
gePOKEd hebben, net zoals we dat bij het vorige programma gedaan 
hebben. 

Probeer de ogen en het vlak rond de ogen een kleur te geven door 
de INK-opdracht te gebruiken. Wellicht kunt u de ogen ook laten 
schitteren met de FLASH-opdracht? 





100 REM programma 19 vliegenko 
195 CLEAR 88550 8 i 
ìig CLS 

128 LET u=i28: LET vnaz 

139 LET heò,5: LET ke30: LET ra 
=PI/188 

140 FOR we92 TO 452 STEP 3 

180 LET tuwsrd: GO SUB 1800 
189 LET xm INT (U+xéh) 

170 LET YaINT (v+y*hì 

160 IF wa9@ THEN LET xiex: LET 
yiey: GO TO 229 

398 LEE THEN 

28 PLOT x1,y1l: DRAW xe-x1,ya- 
W 

210 LET xlax2: LET vieya 

280 NEXT w 

E21 CIRCLE 150,90,30 

222 CIRCLE 88, do, 10 

223 PLOT 1e@ óa 

224 RAÄNDOMIZÉ UsR eae00 

228 PLOT 78 92 

225 RANDOMIZE Ur 64820 
220 IF INKEYga"" THEN Go TO 29 
242 STOP 

250 
1809 LET MaktBIN (240 +(2.B+C08 
1818 LET yakteacos (oet 
1920 RETURN 








Appendix 117 


Programma 21 


Allereerst geven we een paar voorbeelden van wat je zoal met pro- 
gramma 21 kunt tekenen. Je zou dit wellicht teomputerkunst!' mogen 


noemen. 


a=1/b=-100/c=-200 
























a=3/b=100/c=200 
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Tekst op het grafische scherm 


Het zou leuk zijn op het grafische scherm de waarden voor a, b 

en c af te drukken naast de tekening. Bevalt zo'n tekening, dan 
kunnen we direet de waarden van a, b en c noteren. Op de ZX 
Spectrum is dit helemaal niets bijzonders. We kunnen zonder veel 
moeite tekst en graphics op het scherm afbeelden; dat hebben we 
in programma 39 immers al gezien. We volstaan daarom met het 
geven van programma 21, waarin PRINT-opdrachten zijn opgenomen 
voor het afdrukken van de waarde van a, b en c op het ‘grafische! 
scherm. 


122 REM programma @i 
symmetrische krommen 


LS 
128 SRBur “Toets a,b,c in " 
30 LET u=ië8: LET va88 








49 LET haô.5: LET rduPI/180: L 
ET ka8? 
158 CLS 
189 FOR wad TO 360 
zie LET tawerd: LET rak#SIN (c 
pies LET X&2INT (U+r#&COS (att) + 
„JSO LET WESINT WarksiN (beu + 
209 IF wzQ THEN LET XieX@: LET 
wizya; GO TO 23 
„21e PLOT X1,i: DRAU wax gan 
820 LET Xisxe: LET 
Er et a 
hat bre 
238 PRINT AT 80,1;"b ie 
Beo TP INKEv ant Hnên: Gbro aso 
9 2: 
250 SToP 
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Programma 25 


De auteur Marcel Sutter heeft onlangs in het tijdschrift Mikro+Klein- 
computer een verbeterde versie van het 'bolprogramma' 25 gepubli- 
ceerd. We laten dit nu zien. Het programma tekent bollen in allerlei 
standen en met verborgen lijnen. De draaihnoeken om de x-, y- en 
z-as alsmede de 'afstand' tussen de breedte- en lengtelijnen kunnen 
worden ingevoerd. Eerst volgt het programma, daarna enkele voor- 


beelden. 


5@ REM bol met hidden Lines en 
dräaing 
182 CLS 
185 
118 PRINT MET LENGTE EN BR 
EEDTESIEKELS 






PRINT 

dao, ENEUT “ORAAINOEK OM X-A5: 
flo INPUT “DRAAIKOEK OM Y-as: ” 
8ao INPUT “DRAAIHOEK DM Z-RS: " 
16: PRINT 

êRs INPUT "AFSTAND LIJNEN (19-4 


5: “id 

230 LET us128; LET va8 LET fm 
87: LET bmaPI/188: LET ha®, 
240 LET s1eSIN (akbh): LET S2sS 
IN (bébm): LET s3 teh) 

253 LET c1aCO8 (atb): LET c@aC 
9% (bbm): LET c92C0S (c#bm) 


5 
260 REM ROTATIEMATRIX BEREKENEN 
279 LET axmc2rr3: LET aym-cRHEI 
: LET azsse 
280 LET bxacles3e, 





sinsaecs 
290 LET byacitc3-sitsans?: LET 
bim-s lc 
308 CET cxuslasdeclesRecd 
310 LET Cyesincd+ciesaesd: LET 
EANCIMER 
438 REM OMTREK CIRKEL TEKENEN 
405 CL$ 
PEI CIRCLE u‚v‚r 
888 REM LENGTECIRKELS TEKENEN 
12 FOR (w@ TO 188-d STEP d 
515 LET f1is® 
523 FOR pa@ TO 360 STEP 5 
530 GO SUB 1908: REM XX,YY,ZZ 
Sao 
540 T' vy2® THEN LET f2n8: LE 
T_fizs®: GO TO 589 
Eee LET XbeINT (U4xXth): LET 
WBRINT (v+zZ+h): LET fiel 


see IF fimd THEN LET Araxb: L 
680 


er gtevbr EET bei: BOTS 

57 PLOT x1,yl: DRAU Xb-XL,vb 
zuk LET Klmeh: CET yiayb t 

H 


582 NEXT Pp 
588 NEXT L 


IN 
bo 
S 


Appendix 


RKELS TEKENEN 
Ö 9@-d STEP d 


60 STEP 5 
REM XX,YY,ZZ 


EN LET f2=0: LE 


68 
INT (U+XX+h): LET 
): LET fas 


Se An 


EA 


PEOEKENE -E-OHONINAGR UC ÁGER PUREE ARA 
SET ODD FP GADPPSODUWEN Oy Ar MGUEPEN 


Sessa 5 


er 


IF kalk THEN GO TO 798 


STO 
EM BOLCOORDINATEN => CARTE 
COORDINATEN 

ET xer#C05 (p&bm)#CO8 (Leb 
LET garxG0S (pabm)+SIN (Leb 
LET zar+SIN (p&bm) 

REN ON 2) DRAAIEN TOT P( 


LET nxsaxanvayeyazer 
LEL WjabX Xb bz Hz 
RE 


SODAESALEN- SC SHOPEDESASSH 
2 
6 
— 
T 
me 
zi 
Tad 
5) 
á 
x 
es 
u 
x 
T 
e 


R 
HE 
L| 


T ZEECXHKACH EN ACH HE 
TURN 


$ 
EN 
u 
EN 
L 
1 
1 
EN 
Xx 
1 
1 
1 
1 





draaihoek om x-as: 90 
draaihoek om y-as: 30 
draaihoek om z-as: 40 
afstand lijnen : 30 
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draaihoek om x-as: 0 
draaihoek om y-as: 40 
draaïhoek om z-as: 60 
afstand lijnen 2 30 


draaihoek om x-as: -50 
draaihoek om y-as: 20 
draaihoek om z-as: 60 
afstand lijnen 45 
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Programma 36 


We gaan de kaart van Zwitserland wat 
het gewijzigde programma. 


122 


195 
112 
128 
215: 
138 


na 98 
AA Mitsertand 


LET hed,8: LET U 
CET ed 

READ 

LET Le NT tuekexth): 


LET y 
EN LET 
LOT x1,175-yl: DRAU xR-x1 
T KAnGE LET yisyà 

X) 

*® ZIEN Go TO 159 
THEN 9 TO 198 
“THEN Bor 188,88: 
4800: LET aKaINK 


6 
THEN LET Url28: L 


mm 
xD 


RD <01 
De 
GZ 
U 
“oo 
402 
OH Fel 
VIND siNuw 
Oer Omas 


se OBTN, 2 
419 da ATA8 


Zj: 107,70, 104,7 
81,104,86, 105,9 
‚109,198,195,10 
102,117,108,118 
128,110,139,10 
44, 78,154,77,1 
93,76, 177,70, 1 
„52,171,61,167,56,1 
6,94,162,94, 157,48, 1 
'1g6,45, 199,48, 192, c0,1 
LaR,11,119,12,114,20,1 
'1oa,sB,94,4n,88,35,90, 
75,15,66,19,69,14,52,1 
37,29,39,96,97,40,30,4 
16,38,18,39,19,29,6,28 
11,94,13,40,10,45,12,5 
26,79,30,73,48,94, 42,9 


O+ 0 


„ 
@ 


a 
80 
@, 


tene 


PERPPGENPEPEOEPS 


LJ 
@ 


EJ 


vl 
BEREHARPS VD 


BAbOPEEEP GAG 


@- OPE pe 
Te 


: 
1 
67 
1 


st 
Û 


@ 
POGO OPGHO- Te T- 0 Oe U TO URG. Oe 


SDuD 
im À 
P 


N 


@ 
© 


BROIRA EP BSO 


Ln 


re 


a, 
DOLAPODALNPEDOD- DOD- 


L 
© 


Ed 
Falsh ig de ANdade doede Ja 


Ka 
pel 
» 


DADEPE ID De Dr De De De De De DOD 


Li 
4 A 
» B 


PN 
en POR A PDRORNAGGHGRGRG GE d- 


SONG PROPN- Ge OG de On Reeder 
BEAGLES OSE NERENEVSAGAENELONEGSEEE 


"GQ 


verfraaien. Eerst geven we 








Appendix 
&20 DATA _54,102,53,99,61,98,59, 
vigebidee, 
€30 DATA -4O,19,-10,10,-19,40,1 
28 DATA 19,19,40,10,40,-10,10, 
B12 DATA 10,-40,-10,-40,-10,-10 
"826 DATA -19,-40,10,0,0 


Het verschil met het programma op p.97 is dat na regel 190 een 
aantal opdrachten is toegevoegd. Door het intikken van een 'p', 
‘yv! of 'k' gebeuren er nu verschillende dingen. Draai het pro- 
gramma en kijk of u ook de hierbovenstaande illustratie op uw 




















scherm kunt krijgen, 
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Het zal duidelijk zijn dat de in dit boek gegeven programma's onbe- 
perkt uitgebreid kunnen worden. Ze kunnen vast ook verbeterd 
worden; de invoerwaarden zouden bijvoorbeeld gecontroleerd kun- 
nen worden, zodat het programma niet ergens halverwege afbreekt 
door een foutieve schermcoördinaat. We hebben dat allemaal niet 
gedaan, omdat we de programma's klein wilden houden, waardoor 
alle aandacht op de tekentechniek gericht blijft. 
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Over de inhoud van dit boek 


Dit boek bevat 40 grafische 
programma’s voor het programmeren 
met hoge resolutie in ZX Spectrum BASIC. 
Dit houdt in dat het scherm van uw 
ZX Spectrum computer verdeeld wordt in 
256 puntjes horizontaal en 176 puntjes 
verticaal. Op dit grafische scherm kunnen 
fraaie tekeningen gemaakt worden. 

Bij bijna elk programma wordt de 
nodige tekst en uitleg gegeven. Daarnaast 
bevat het boek veel illustraties van dat- 
gene wat u op het scherm kunt 
verwachten. Alle programma's zijn getest 
op een ZX Spectrum. De programmatekst 
is direct vanuit het geheugen van de 
ZX Spectrum op een printer afgedrukt en 
deze afdrukken zijn rechtstreeks in het 
boek opgenomen. 

Laat u niet afschrikken door de wat 
wiskundige uitleg bij de programma's. 
Ook voor niet-wiskundige Spectrum- 
bezitters is dit boek bedoeld. Het alleen 
maar intikken en draaien van de 
programma’s geeft al zulke fraaie 
beelden, dat het bestuderen van de 
erachterliggende theorie geenszins 


De illustratie op de voorkant is 
gemaakt met het onderstaande 
programma 


Mú 
L 
CLS 
123 LET vaag, : LET y2=175 
@ TO 25 TEP 36 
148 FOR bz@ TO 255 STEP 36 


NEXT à 
193 IF INKEY&z="" THEN GO TO 198 
STOP 


noodzakelijk is, alhoewel de ‘kijker’ hier- 
toe wel geprikkeld wordt. 

Naast het tekenen van allerlei 
eenvoudige en meer ingewikkelde 
functies (voor leraren en scholieren) 
worden de beginselen van het drie- 
dimensionaal tekenen uitgelegd. Ook de 
techniek van de verborgen lijnen (hidden 
lines) wordt gebruikt. Vijf programma's in 
BASIC laten zien wat voor figuren je met 
de taal LOGO kunt maken. Hieronder zijn 
ook de beroemde turtle-graphics. 

Ook bevat het boek vijf educatieve toepas- 
singsprogramma’s, waarin met graphics 
gewerkt wordt. 

In een appendix wordt een aantal 
programma’s verfraaid en uitgebreid om 
de lezer de vele mogelijkheden van ZX 
Spectrum graphics te laten zien. 

Deze appendix bevat een programma 
voor het inlezen van een paint routine 

(in machine code), die niet standaard in 
ZX Spectrum BASIC voorhanden is. 

Met deze routine kunnen vlakken worden 
gevuld of ingekleurd. 
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Over de inhoud van dit boek 


Dit boek bevat 40 grafische 
programma’s voor het programmeren 
met hoge resolutie in ZX Spectrum BASIC. 
Dit houdt in dat het scherm van uw 
ZX Spectrum computer verdeeld wordt in 
256 puntjes horizontaal en 176 puntjes 
verticaal. Op dit grafische scherm kunnen 
fraaie tekeningen gemaakt worden. 

Bij bijna elk programma wordt de 
nodige tekst en uitleg gegeven. Daarnaast 
bevat het boek veel illustraties van dat- 
gene wat u op het scherm kunt 
verwachten. Alle programma's zijn getest 
op een ZX Spectrum. De programmatekst 
is direct vanuit het geheugen van de 
ZX Spectrum op een printer afgedrukt en 
deze afdrukken zijn rechtstreeks in het 
boek opgenomen. 

Laat u niet afschrikken door de wat 
wiskundige uitleg bij de programma's. 
Ook voor niet-wiskundige Spectrum- 
bezitters is dit boek bedoeld. Het alleen 
maar intikken en draaien van de 
programma's geeft al zulke fraaie 
beelden, dat het bestuderen van de 
erachterliggende theorie geenszins 


De illustratie op de voorkant is 
gemaakt met het onderstaande 
programma 


CLS 
123 LET y1=@ : LET y2=175 
130 FOR az@ TO 255 STEP 36 
148 FOR bz@ TO 255 STEP 36 
159 PLOT a,y1l 


noodzakelijk is, alhoewel de ‘kijker’ hier- 
toe wel geprikkeld wordt. 

Naast het tekenen van allerlei 
eenvoudige en meer ingewikkelde 
functies (voor leraren en scholieren) 
worden de beginselen van het drie- 
dimensionaal tekenen uitgelegd. Ook de 
techniek van de verborgen lijnen (hidden 
lines) wordt gebruikt. Vijf programma's in 
BASIC laten zien wat voor figuren je met 
de taal LOGO kunt maken. Hieronder zijn 
ook de beroemde turtle- graphics. 

Ook bevat het boek vijf educatieve toepas- 
singsprogramma's, waarin met graphics 
gewerkt wordt. 

In een appendix wordt een aantal 
programma's verfraaid en uitgebreid om 
de lezer de vele mogelijkheden van ZX 
Spectrum graphics te laten zien. 

Deze appendix bevat een programma 
voor het inlezen van een paint routine 

(in machine code), die niet standaard in 
ZX Spectrum BASIC voorhanden is. 

Met deze routine kunnen vlakken worden 
gevuld of ingekleurd. 


173 NEX 

189 NEXT a 

199 IF INKEY&z="" THEN GO TO 199 
22 STOP 
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